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Problem in Self-Heating of a Spherical Body 


S. M. Genensky* 


An analytic steady-state solution is developed for a spherical body in which heat is 
generated according to a first-order unimolecular-reaction law, and lost at the surface in 


accordance with Newton’s law of cooling. 


The temperature within the sphere depends 


chiefly on the radial distance from the center, but also on the ambient temperature, the 
surface heat-loss coefficient, and the material properties, which are assumed constant. 


l. Introduction 


The analysis described was pursued in connection 
with an investigation into the self-heating of fibrous 
materials carried on at the National Bureau of 
Standards.’ 

Comparison of the results of this analysis with the 
experimental results indicated that the assumption of 
a single first-order unimolecular reaction was an 
for the material investigated 
experimentally. However, the analysis serves to 
indicate the relative importance of the various 
parameters entering into the problem, and may be 
helpful in indicating the kind of data that would be 
of greatest use in further investigations in which the 
assumptions may be applicable. 

A homogeneous sphere of radius B is generating 
heat, under steady-state conditions, in accordance 
with a first-order unimolecular-reaction law, Ae~ ‘“/*”, 
and losing heat from its surface to an atmosphere at 
constant temperature, 7,, according to Newton’s 
law of cooling. The temperature of the sphere, 7(r), 
is independent of time and depends upon the radial 
distance, 7, from the center of the sphere. Within 
the sphere, transfer of heat by means other than 
conduction is considered negligible. The thermal 
conductivity, *; constant, /; rate of reaction- 
frequency product, A; and activation energy, £, 
of the material of the sphere are assumed to be known 
constants, as are the temperature, 7), and gradient 
dT/dr=0, at the center of the sphere. The problem 
is to find the surface temperature of the sphere, 
Ty, at r=B, and the heat-transfer coefficient, h, 
between the surface of the sphere and a surrounding 
atmosphere. 


oversimplification 


us 


2. Analysis 


Mathematically the problem becomes 


eT 2d’ / 
| Sat r 1 ta kT—0 (O<) B), 


*Present address: Brown University, Providence, R. I 
D. Mitchell, New light on self-ignition, Natl. Fire 
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Then eq (1) to (4) become 


’V 2dV 


= ae (5) 


dn? 
(6) 


dV 


QO at 7 
dn " 4) 


dV BhR 
i dy ) 1 7 I 


It is further assumed that (a) the temperature 
gradient within the sphere exists and is continuous 
throughout the interval O<»<1, and (b) that the 
absolute value of the difference between Tn), 
0<n<1,and 7, is small in comparison to 7». -Under 
these restrictions, a technique presented by Cham- 
bre * proves useful to observe that 


[T's af? T,|- (3) 


] V i as V\' 
AV | u-4 )| ve (a)? (9) 


where AV=V,—V. Now, because |AV/V,!< 1, the 
series converges, and for a sufficiently large positive 
N 


(1/V,)>3(AV/V,)' is a very good approx- 


i=0 


integer, N, 


imation to 1/V. 
Therefore, (5) may be written approximately as 


N 
1 ~fAV\i 
@V .2dV. .. v2) 

+(e ‘i=0 0 O<n<1) 
n dn 


0 
dy (10) 


?P. L. Chambré, On the solution of the Poisson-Boeltzman equation with ap- 
plication to the theory of thermal explosions, J. Chem. Phys. 20, 1795 (1952). 








Now 
R/AV\' s V\ : ; ; 
Ss . _ ( 11) | : 
2 (+) 2 (1-7) = +1)—B,V+ B21 
l i > y-. 
where 
B < n! | y 
‘N : ; ) A 
ae — 1 >) | ae 
and in the special case V= V, 
(N+1)—B,V+B,V? 1)*ByV*=1 
V=0,1,2 
Therefore, (10) may be written 
av 2dV 4, Ve, & (Bi V—B,V24 1) *ByV 
dy’ ' n dn 
(11) 
Because (11) is analytic and regular except at »=0, 


and further as dV/dy=0 at »=0, the solution of (11 


may be written in the form 


V=V(n) 2ain' (O<n<l). (12 
Recalling (6), eq (12) at n=0 yields 
TR 
V,.=V(0)=——=a (13 
hk 


and further differentiating (12) once and considering 
(7), it is found that 

a,=0. (14) 

Differentiating (12) twice, substituting these first 

two derivatives into eq (11) and rearranging terms: 


i ee ae | BV —By\ 1I)’ByVS . 
3 (72+ )a;n' . De' , (15) 
V+1 
where D=( e and the a, are given by 
| (BV —BsN IYBy VN 
1 d 1 pe | 
=e 
m (7 r) (a—2)! dy‘ ~* 
(16 
for 7=0, and 1=2,3,4, 
Using the coefficients obtained from (13), (14), 
and (16) and substituting them into eq (12), V can 


be computed for 0<n<1. 

In particular, for 7»=1, eq (12) gives the dimension- 
less surface temperature V(1), and thus recalling 
that 7(n) =(2£/R)V(»), the surface temperature, 7(1), 
which equals 7',, is easily found. Differentiating 
eq (12) once and evaluating this derivative at »=1, 
h can then be found by solving eq (8), for all the 
other factors involved in this equation are either 
given or have been computed. 
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The first six nonzero coefficients have been evaly. 
ated and are 


do= Vo 
(’ l 
a — ih ' } 
6 
’ 
( 2 
Qs 120172‘ j 
hy = : 
" eee aes 
“o1008V2115 Vv, 34‘ 
: > ’ { } 
Cc" 122¢" 122¢’, 10 
as=>-: 4 ott? 6~ CO rT = Of 
51840Vil 63V2 21, 3 | 
( 503207 S428C" 1486¢€' 560 
Gun - i eee aca —ye 
" 4435200V3 |. 405V3 135V2 " 27V, 9 |' 
Observe that a.,_; ((=1,2,3, .) are zero. This 


follows from a,;=0 and the spherical symmetry of 
the problem. 

Figure 1 is a plot of V(1), as a function of C and 
1/V,. where C was allowed to vary over the range 
10-'<('<10" and 1/V, took on selected values in 
the range 5<1/V,<100. The points were obtained 
by using | 


<a a 
(17) 


over its range of applicability for the values of C 
and 1/V, within the above limits. 

In the foregoing, it has been assumed that the | 
temperature at the center of the sphere 7) is known, 
and also that the physical condition indicated by 
eq (8) is satisfied. ‘ 

The more common steady-state problem of finding 
the temperature at the center of a sphere that gener- 
ates heat according to a first-order unimolecular 
reaction law and loses heat according to Newton’s 
law of cooling may be solved by using the analysis 
described above. In this case, h, PR, k, A, h, T,, 
T,, and B are assumed known. 

If the sphere being considered really satisfies the 
assumptions of the problem, then its temperature 
must satisfy eq (1). Because eq (1) is of the second 
order, only two boundary conditions can be pre- 
assigned. Further, because the temperature at the 
center of the sphere must remain finite, one of these 
preassigned boundary conditions must be d7/dr=0 


*- 


at r=0. This leaves but one free boundary con- 
dition. However, at the surface of the sphere both 
=) Di ili ave 
——|7 2 (18) 
( ad = }- | B 
and 
T=7T, at r=B (19) * 


must be satisfied; and because only one of these two 
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boundary conditions may be formally assigned, it 
must be assumed that the other is consistent with 
the problem. Therefore, figure 1, which represents 
eq (17), may be used to solve this problem. Because 
R, Ek, A, and B are known, both C and V(1) (recall 
that 7,=FEV(1)/R) can be computed. Thus, by 
using figure 1, 1/V) may be found by interpolation, 
provided the point (C, V(1)) under consideration 
lies within the region of applicability of eq (17), as 
indicated above. Because the center temperature 
of the sphere is 7), that is LV,/P, the problem is 
solved. 

Here no use was made of condition (18), and as 
mentioned earlier, this condition must be satisfied 
independently, if the assumptions made are fulfilled. 
Thus the degree of consistency between conditions 
(18) and (19) serves as an indication of the applica- 
bility of the assumptions made in this analysis. 


WASHINGTON, March 19, 1954. 
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Properties of Arsenic Sulfide Glass 
Francis W. Glaze, Douglas H. Blackburn, Jerome S. Osmalov,? Donald Hubbard, and 


Mason H. 


glass were 
propert ies 
expa 


sulfide 
Among the 
electrode funetion, 


Samples of arsenic 


physical properties 
durabilitv, hvdrogen 
hardness, 

The 
the infrared 
its physical properties are inferior to those of « 


main mterest 


investigated to 


Black * 


determine their chemical and 
are the hygroscopicity, chemical 
elastic moduli, modulus of rupture, 


reported 
nsivity, 


transmittance, and indices of refraction. 
in arsenic sulfide glass, 
It is transparent to infrared energy out to approximately 13 microns, 


for the present, is for its transmittance in 
3 Although 


rdinary silicate glasses, they are superior to 


those of some of the available materials used in infrared opties. 


1. Introduction 


Arsenic sulfide (As.S;) glass, although reported in 
1870 [1],! has received consideration only in recent 
years for its spectral transparency to 
infrared radiation. It has a usable transmittance 
of radiant energy out to 12.0 or 12.5 uw, depending 


desirable 


on thickness [2,3] 

The present investigation undertaken — to 
determine the chemical and physical properties of 
this glass and also to obtain pertinent information 
as to its suitability for an element of an 
optical system when exposed to the atmosphere. 


was 


use as 


Much of the material contained in this paper has 
been the subject of a report to the U.S. Air Force {4}. 


2. Preparation of the Glass 


The arsenic sulfide glass for most of the tests was 
made in the apparatus shown in figure 1. The 
arsenic and sulfur, in stoichiometric ratio of 2 to 3 
plus a 0.25 percent excess of sulfur, were allowed to 
react in the reaction tube in an atmosphere’ of 
nitrogen, and then the arsenic sulfide was distilled 
at approximately 735° C into the accumulator. <A 


typical time-temperature schedule for a 5,000-g 
batch of arsenic sulfide is given in table 1. A slug 


of glass, chilled by means of a stream of nitrogen 
impinging on the outlet tube, retained the distillate 
in the accumulator. Nitrogen was used throughout 
this operation to prevent oxidation. When the 
distillation was completed, the glass in the accumu- 
lator was stirred to homogenize it. Heat was then 
applied to the outlet tube, the stirrer was stopped, 
and the first slug of glass was caught in the tube 


This work w DI fir t Dep ent the Navy, Bureau 4 . ps 
of Ships, and later by the United States Air Force, Wright Air Development | suspended from the bottom of the molding oven. 
Center. Materials Laboratory a hl . & . 

The mold was then quickly moved into place and 

Present adadre Eastman Kodak ¢ Rocl N. ¥ . 

Dect 1 { patent has been applied for on this one-step procedure for making arsenic 

Figure n bracket | ! nees at the end of this paper sulfide glass 

TABLE 1 Arsen sulfide production run 
Operation sheet 
Mi Accumulator Collecting oven 
Tir ( Outlet Remarks 
tubs 
| B Top Bottom lop Bottom 
( ( ( ( ( ( ( ( 
S00 Current on, 
a-) 14 im 2 430) 270 6 
i. iM is $4 1) $55 100 
10:00 | au 12? 40) 425 160 
wu) 2 44 40 is) 245 145 140) 
11:00) ; i 140) 570 20 300 225 22) Stirrer on (approximately 240 rpm). 
11:30 til 2 120) wh 5 320 255 25) Distillation completed. 
2-00 5 40 325 D5 20) 
12: 3 wi) 40 325 20) 255 
1:00 Hit) 540 320 265 260) Cut stirring speed to approximately 180 rpm 
0) M5 25 $25 25 20) 
2-00 525 WO 325 pie) 255 
2-30) 51S 100 320 260 D5 Cut stirring speed to approximately 120 rpm 
On ALS 100) 320 4M) 2h) 
Ww a) 10) 325 iM) 7%) 
1-00) 17 1) 320) iM) mu) Cast 
Ro 1 
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the slab cast. After flow practically ceased (about 
3 min), a glass plate was placed over the mold and 
the mold and contents moved to an annealing oven 
at 185° C. It was held at this temperature for 
12 to 16 hr and then slowly cooled to room temper- 


ature. For pieces weighing 2.0 to 2.5 kg, the whole 
annealing operation took 34 to 3% days. No 


attempt was made 
the annealing oven. 

Although the apparatus shown in figure 1 per- 
mitted the use of either technical-grade or distilled 
arsenic without appreciably altering the quality of 
the glass, the large amount of residue left behind by 
the technical grade attacked the 


to control the atmosphere in 


















The transmittances of arsenic sulfide glasses made 
at the Bureau by similar procedures and of materials | 
of comparable purity will agree within +2 percen 
between 2 and 8 yu, except possibly in the region of | 
the water absorption band at about 2.8 x. . 

Nothing is known about the details of preparation | 
of the commercial specimens. The commercial and 
experimental glasses are appreciably different in ap. } 
pearance and spectral transmittance in the visible 
and very near infrared as shown by data obtained | 
with a Beckman quartz spectrometer, model DU | 
(fig. 2). 











reaction tube, TTY 
thereby shortening its life. J 
Some of the experiments were made on samples of 4 
glass obtained from a commercial source ® and : } 
others on experimental melts prepared at the Bureau é 
by the one-step process described above. 2 
The large thermal expansivity of the As.S, glass | 2 
presented a distinct hazard to molding blanks suc- | & 
cessfully. The blanks, while being formed in Pyrex | 4 MMERCIA 
7740 or Vycor, often stuck to the walls and bottom | 2 
of the molds. Upon cooling, the differential expan- | © , ' 
sion between the mold and specimen and the adhesion 2 
of the arsenic sulfide glass to the walls of the mold | @ 
frequently resulted in severe fracture of both mold | ° 
and specimen. This difficulty was overcome by 
lining the mold with aluminum foil, which eliminated ls 
sticking and allowed complete freedom of adjustment f 
to the large dimensional changes that occur during oF 2.1 te | oo | 
cooling. WAVE LENGTH, MILLIMICRONS 
Fictre 2. Comparison of the spectral transmittance of rail 
* The commercial samples were projuce1 by the American Optical Co., Sout! sam ple s of arsenic sulfide glass from SCO to 1,000 Mig 
bridge, Mass Specimens 6 mm thick 
5 
; 
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u $4 eat Joint STIRRER 
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( 64mm 0.4. 24 in.ig) Tr 
. / N, OUT 
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(40m mM 0.4.) Yj Ww 
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Fiaure | Schematic drau ing ol the apparatus 
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FIGURE 3 


3. Chemical Properties | 


3.1. Hygroscopicity 

The hygroscopicity and chemical-durability data 
were obtained by the powder and interferometer 
procedures described in earlier publications [5, 6]. 
In crushing and screening the powdered samples 
(approximately 1.5 g that passed a No. 150 Tyler 
standard sieve) for hygroscopicity determinations, 
the “fines’’ adhered conspicuously to the mortar, 
sieve, pan, walls of the hood, and the skin and 
clothing of the operator. Because the material 
toxic, suitable precautions should be taken during 


IS 


this operation. 


Figure 3, plotted from data given in table 2, shows 
the hygroscopicity curves (water sorbed versus 
hours of exposure to 98% relative humidity) for 
specimens of arsenic sulfide glass from two different 
sources, compared with those for fused SiQ,, Vycor, 
Pyrex 7740, a window glass, and Corning 015. 
These silicate furnish a reference index 
covering a broad cross section of commercial usage. 
Although the results obtained on the two samples 
differ from each other, both are considerably less 
hygroscopic than Pyrex 7740. One of the samples 
even compares favorably with fused SiO, and Vycor, 
although this should not be interpreted as an indica- 
tion of good chemical durability, as shown in figure 4 


glasses 


urability of arsenic sulfide glasses compared with similar prop- 


erties of typical commercial glasses. 














TABLE 2 Hygroscopicity, voltage characteristics, and chemical d 
W ater sorbed Voltage departure at pH 
Gl ‘ 
1 hr 2 hy 2 1.1 6.0 8.2 10 2 
é mM n" ne mi mer mi 
Fused SiO 6.3 x ‘ ‘ ‘ 
V veor s. 1 1 1 “ * 
| t 0 0 105 2 342 447 
AsSs(N BS). 0 G5 1US 303 438 
0 120 141 332 421 
A{seS;(commercial 13 19 
Pyrex 7740 16 2s 0 77 If 28S 
W indow 23 7% 
Corning 01 61 12 0 0 0 0 0 
* No defir YH response ND. not detectable. A. attack iS. owe 
70 T 
9 CORNING 5 
« 60 7 os 
- / 
j 
WwW / 4 
= WINDOW 
— 
2 
wy 20 
oO j 
Lv 
a / 
= / 
U j 
~ 
”) 40 
2 
a / f 
wo) / 
ould 
+ wD / / 
= / F + PYREX 7740 
’ / , Pe 
ra) / / i 
WwW a 
@ / y 
S y “ 
Oo 20 j / ” ii < 
a j A j ener 
y, Yi id COMMERCIAL 
« | | Sf + 
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om / 4 / 
a ff a P NBS 
F1or-/ AZ vysoR — O° 
Pa - ——eEE 
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a ae 
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Figure 4. Chemical durability of arsenic sulfide glass com- 


pared with similar data for Corning 015 glass and fused 
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where fused silica and arsenic sulfide glasses are com- 
pared. The differences between these two As,S, 
glasses, as shown by the hygroscopicity data, and 
which were also evident in the transmittance curves 
of figure 2, may be the result of impurities or varia- 
tions in the arsenic-sulfur ratio 


3.2. Chemical Durability 


a. Swelling and Attack 


The chemical durability of the sulfide 
glass was reported as the surface alteration in inter- 
ference fringes brought about by exposure of optically 
flat surfaces to various solutions under controlled 
conditions of temperature, time, and pH. Test 
solutions chosen were the Britton-Robinson universal 
buffers [7] over the range pH 2 to pH 11.8, hydro- 
fluorie acid, and sulfuric acid. 

There is general belief that the lower the hygro- 
scopicity of a glass, the greater the ability of the 
glass surface to remain optically clear [8,9]. Judged 
on this basis from figure 3 and from the examination 
of polished samples which have been stored in the 
laboratory for approximately four years, both 
samples of As.S,; would be expected to maintain 
optically clear surfaces over a considerable period. 

Although the above indicates that arsenic sulfide 
glass has satisfactory serviceability, it rates very 
low in chemical durability as judged from figure 4, 
plotted from data in table 2. In fact, it is inferior 
to the very hygroscopic, pH responsive Corning 015 
glass. Among the most interesting features of the 
curve for the chemical durability of As,S, glass, is 
the appreciable swelling in the acid range below 
pH 6. Experimentally, the data between pH 2 and 
pH 6 were readily reproducible. Ai pH 8, however, 
the attacks were variable, showing wide disagree- 
ment in results on duplicate samples ranging from 
difficultly detectable to large attack, as shown by 
the arrow in figure 4. Results in this region were 
highly dependent on just when the swollen laver 
sloughed off or dissolved. This is characteristic of 
glasses at pH values near which the surface altera- 
tion is passing through the transition from swelling 
to attack (approximately pH 7). The As,S 
exhibited swelling of 1/2 fringe (comparable with 
the swelling in acid buffer solutions) when exposed 
to distilled water for 6 hr at 80° C At pH 11.8, 
the attack was so rapid that the exposure time had 
to be reduced to 1 min at room temperature. The 
value plotted was calculated on the assumption that 
the rate of attack doubles for each 10 deg C rise ia 
temperature [6]. This rapid attack in alkaline solu- 
tions is in line with the high solubility of As.S; in such 
solutions [10]. The swelling of the As,S, glass in the 
pu range of 2 to 6 is similar to that observed on the 
Corning 015 glass. 

Table 3 gives the results obtained by progressively 
increasing the concentration of electrolyte in’ the 
ambient solution by the use of H,SO,. Figure 5 
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glass 
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shows the over-all curve for chemical durability 
obtained by splicing these results onto the data 
furnished by the Britton universal buffer mixtures. 
Below pH —1, the swelling is greatly repressed, with 
little or no surface alteration detectable in the 
majority of tests between pH—2.45 and pH —8.2/ 
However, 2 specimens (plotted as dashed lines in 
fig. 5) out of more than 15, showed vigorous attack 
in this ‘‘super acid” region. As similar attacks were 
obtained upon the addition of small amounts of 
ferric salts as an oxidizing agent,* it was concluded 
that residual amounts of the ferric compounds pres- 
ent from the polishing rouge were probably respon- 


Although such observed value is pH 8.2 cannot be rationalized in the Arrhe- 
nius sense, i. e. they cannot possibly represent hydrogen-ion concentration in 
moles per liter [11], nevertheless these activities as indicated by the hydrogen 
electrode : saturated KC] cell are accepted and plotted for convenience 


* Additional tests using HCl and FeC'; mixtures gave attacks of 4! nterference 
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The chemical-durability for Pyrex 7740 is 
included i) figure 5 as a reference standard. 


curve 


b. Arsenic Sulfide Glass in Hydrofluoric Acid 


Silicate glasses are generally attacked vigorously 
by hydrofluoric acid with the formation of soluble 
products and volatile Sify. The chemical durability 
of arsenic sulfide glass ia HF solutions is of interest 
hecause silica. Preliminary 
experiments showed no displacement of fringes or 
other detectable alteration of polished flats upon 
exposure to tS-percent Hk for 6 hr at room tem- 
perature, At elevated temperatures, a perceptible 
swelling Was observed after | hr of exposure. Figure 
§ shows the results obtained upon exposing speci- 
mens to 48-percent HF at 80° C for various periods 
of time up to 6 hr. Swelling increased with time 
of exposure up to and including 4 hr, after which the 
swollen laver either became mechanically weak aad 
sloughed off, or oxidized, forming soluble products. 
This attack, after 6 hr, was repeated on five different 
specimens, With similar results in each case. 
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3.3. pH Response of Arsenic Sulfide Glass Electrodes 


For studying the pH response, electrodes were pre- 
pared by coating an aluminum rod (approximately 
2-mm diameter) with a film of the As.S; glass. This 
was accomplished by dipping the rod into the molten 
withdrawing and cooling sufficiently 
prevent ignition of the arsenic sulfide 
\luminum was used because it has ap- 


and 
quickly to 
in the air 


glass 
glass, 


$2GTO0 57 » 


proximately the same coefficient of expansion 
as the arsenic sulfide glass (23.810 and 
23.7 X10-°/°C, respectively), and was the only 
material readily available for making successful 


electrodes of the Thompson metal-filled type [12]. 
This procedure was necessary because the conven- 
tional Cremer bulb type of electrode could not be 
blown {13}. 

All emf measurements were made at room temper- 
ature with a Beckman pH meter model G, using a 
calibrated glass electrode as the reference half cell. 

At the time of preparation, these arsenic sulfide 
electrodes had resistance values greatly exceeding 
1,000 meg, and were incapable of initiating or sus- 
taining readable voltages on the electrometer, i. e., 
they behaved as an open circuit. The resistance 
upon exposure to aqueous solutions fell quickly 
to a few megohms, and the electrodes then behaved 
essentially as punctured electrodes. Their initial 
performance was in accord with the very low hygro- 
scopicity of the arsenic sulfide glass, whereas the 
later electrical response was compatible with the 
excessive swelling exhibited by this glass in acid 
solutions between pH —2 and pH —6,and general over- 
all poor chemical durability in the alkaline pH 
range. 

Figure 7 and table 2 give the voltage departures 
(errors) obtained on three arsenic sulfide glass elec- 
trodes immediately after preparation, compared with 
the corresponding departures of a pH _ responsive 
glass (Corning 015) and the calomel half cell. The 
latter has no electrical response to hydrogen ions. 
From this comparison, it is obvious that the arsenic 
sulfide glass electrodes exhibited little or no response 
to hydrogen ion activity of aqueous solutions. This 
lack of pH response is in full accord with the low 
hygroscopicity and poor chemical durability of 
arsenic sulfide glass [10]. 

In order to establish that the data in figure 7 
represent the performance of the arsenic sulfide, 
and not merely the characteristics of the aluminum 
electrical connector, the voltage characteristics of 
the bare aluminum rod were followed over the range 
pH 2 and pH 11.8. Figure 8 and tables 4 and 5 
show a comparison of the electrode performance 
with the reported chemical durability of Al metal 
[11] and the pH range of precipitation of Al(OH), 
(14]. The interesting feature of the voltage depar- 
ture curve is the direction of drifts of the voltage 
below pH 4 and above pH 11. These drifts are 
toward increasing departure, indicating the electrode 
is being attacked as shown by the durability curve. 
Between pH 4 and pH 11, the aluminum electrode 
tended to indicate the H-ion concentration, with a 
steady reading appearing near pH 8. The latter 
point should be near the pH of minimum solubility 
of Al metal in the buffers used. Inspection of these 
data suffices to establish that the voltage departure 
recorded for the glass electrodes in figure 7 is char- 
acteristic of arsenic sulfide glass and not of aluminum. 
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4. Physical Properties 
4.1. Thermal Expansion 


Expansion data were obtained by means of the 
Saunders modification of the interferometer proce- 
dure for measuring dilation [15]. Figure 9 shows a 
typical curve for the thermal expansion from room 
temperature to the deformation point of arsenic 
sulfide glass. The coefficient of expansion of this 
glass, 23.7 < 107°/° C between 50° and 175° C, is over 
7 times that for Pyrex 7740 (3.3 107°/° C), and 
approximately 2, times that for common window and 
container glass. 


4.2. Elastic Moduli 


The elastic moduli were determined by the pro- 
cedure used by Spinner [16]. Table 6 gives typical 
values for Young’s and shear moduli and Poisson’s 
ratio for annealed arsenic sulfide glass from 20° to 
190° €C Both Young’s and shear moduli decreased 
with increasing temperature, while Poisson’s ratio 
increased. A graph of Young’s and shear moduli 
versus temperature is given in figure 10. 

Arsenic sulfide glass has a Young’s modulus of 162 
kilobars at room temperature (20° C), or approxi- 
mately one-fourth that of ordinary glasses. This is 
very fortunate in respect to the thermal-shock re- 
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Young's and shear moduli versus temperature for 
NBS arsenic sulfide glass. 
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, Young's modulus; —@-, shear modulus. 


endurance, F’, is defined by the expression, 


FR P K 
ak Y se’ 


where / is the tensile strength, @ the coefficient of 
expansion,  Young’s modulus, A thermal conduc- 
tivity, s the density, and ¢ the specific heat [17]. The 
low Young’s modulus helps to counteract, to a cer- 
tain extent, the high coefficient of thermal expansion, 
because the coefficient of thermal endurance varies 
inversely with their product. In spite of the low 
Young’s modulus, the thermal-shock resistance is low 
in comparison with the regular glasses of commerce. 


Young's modulus, shear modulus, and Poisson’s 


TABLE 6 
versus temperature for arsenic sulfide glass 


ratio 


Annealed sample of glass made at N BS [4]. 


Young's Shear Poisson's 
Tempera modulus medulus ratio 
ture (ke (G 
4y Nilohars Nilohars 
1 162.0 62.0 0. 306 
) 160.8 61.9 209 
60 159. 5 61.1 305 
SU 158. 1 60.5 307 
100 156. 7 59.8 310 
120 155. 2 0.2 311 
140 153. 4 58. 7 307 
169 152. 3 48. 1 $l 
180 150. 9 A7.4 314 
100 149.9 57.0 315 


Although flat specimens of the As,S; glass could 
only withstand quenching thermal shocks of less 
than 20° C, the same specimens could be placed on 
en aluminum plate at 300° C without breaking. 
Observations with the polariscope revealed that in 
the latter situation, both the upper and lower sur- 
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faces of the specimens were placed in compression. 
The large expansion of the hot Jower surface caused 
upper 
the same time, as long 
as the main bulk of the slab remains below the defor- 
mation temperature, the lower surface was also held 


the specimen to curl upward, placing the 
surface in compression. At 


in compression. As soon as the deformation 


no breakage can occur. 


4.3. Modulus of Rupture 
The modulus of rupture was determined on 
transverse strength apparatus with midpoint loading 
Values for annealed and 
vylass specimens made at 
Standards are given in 


the National 
table 7 
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ments were made on specimens epproximately 0.4 
by 0.7 by 6.0 in. long, using a 4-in. span. The 
process used in “tempering” or “toughening” arsenic 
sulfide glass consisted of heating the specimens to 
250° C, quenching them in oil at 130° C, and leaving 
the specimens to cool in the oil. As can be seen 
from the table, the “toughening” treatment increased 
the arithmetic mean of the modulus of rupture from 
2,096 to 6,505 psi, or 3.1 times. Although the 
standard deviation was higher for the “toughened”’ 
samples (934 against 461 psi for the annealed 
samples), the coefficient of variation was less (14.36 
to 21.99% 


Lem- 
perature Is reached, the stresses are releas “l and still 


“toughened” arsenic sulfide 
Bureau of 


These mersure- 


4.4. Hardness 


, : - 
The Knoop hardness of annealed arsenic sulfide 


glass was found to be 120 to 125. However, wheg 
these same samples were “‘toughened”’ and agaip 
tested, the Knoop hardness was only lOO to 1]] 
The Knoop hardness of commercially annealed 
soda-lime glass is 520; that for lead glasses varies 
from 370 to 440. depending on the lead content 118] 


4.5. Optical Properties 
a. Transmittance 


The infrared transmittance curve for two thick. 
Hesses of arsenic sulfide glass made at the National 
Bureau of Standards is shown in figure 11. These 
measurements were made by the Radiometry See- 
tion of the Bureau Figure 12 shows the absorbance 
index and surface loss (one surface) calculated from 
the results plotted in figure 11 


19}. As can be 
seen, the maximum transmittance obtainable js 
about 70 percent because of the high reflection 
loss. This high loss is due to the high index of 


refraction, 2.4 1p 
iden 


the infrared. Figure 13 
of the difference in transmittance between a 
sample made at the National Bureau of Standards 
and a commercial sample made by the American 
Optical Company. These same samples were then 
“toughened” and their transmittances again meas- 
ured. 


Ives an 


\ comparison of the transmittance curves 


before and after “toughening” indicated that the 
“toughenipg” process did not affect the trans- 
mittance in the tafrered Also, though there are 


differences between 


samples 


the transmittances of the two 
NBS and commercial), the over-all trans- 
mittances of the two glasses are practically the same 

It will be noticed that there is appreciable dif- 
ference between the transmittance curves 1a figures 
11 and 13. The glass sample used for figure 1] 
resulted from a small melt (230 ¢ of batch); the 
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rransn fiance curves for commercial and 


VAS arsenic s fide glass 


Fictrt 
Annealed, 


CS Tat N BS glass, t=6.00 mm 


sample used for figure 13 resulted from a large melt 
5.000 @ of bateh) made in the 
described prey iously Also, in the case of the small 
melt, doubly distilled arsenic was used; in the large 
melt, singly distilled was used. In the 
small melt, 1 ml of water was added to the batch 
to displace any residual air; in the large melt, none 
was added. The small melt was not cast, but was 
annealed in the reaction tube after homogenizing. 
The differences lil the procedures enumerated above 


one-si ep process 


ursenic 


ean easily account for the differences in the two 
transmittance curves. 
b. Indices of Refraction 
To our knowledge, only one measurement of the 


indices of refraction of arsenic sulfide glass in the 
infrared has been made. That measurement was 
made at the National Bureau of Standards, by the 
Optical Instruments Section, on a prism furnished 
by Walter A. Fraser, at that time working in the 
Research Laboratories of the American Optical 
Co., Stanford, Conn. The dispersion curve 
this prism is shown in figure 14, plotted from the 
results given in table 8 
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lhis glass was produced by Walter A, Fraser while in the employ of the 
American Optical Co., Southbridge, Mass Measurements were made at N BS 
PARLE & Index of refraction of a prism of arsenic sulfide 
glass produced by Walter A, Fraser while at the American 
Optica C'o., Southbridae. Mass. 
Determinations made at NBS 
W elength Index of 
! ron refraction 
( TOO 2. 65822 
H43S4 2 5SM413 
1. OL30S ? 47230 
1. 1280 ” 45897 
$5703 2.44244 
1. 36728 2. 44100 
1. 39506 > 44050 
1. 52952 2.43474 
1. H032 2. 42065 
1. 7092 » 42913 
1. SL307 $20.56 
1. 97009 » 42341 
2. $2542 » 41848 
$. 3033 » 41134 
3. 41588 ) 410608 
$44 2. 39900 
6, 2328 > SU718 
s. HR 2. 38734 
4.724 2. 38029 
11. 035 » 37055 


5. Conclusions 


The very low hvgroscopicity exhibited by arsenic 
sulfide glass is a characteristic favorable to its use 
as an optical element. It is better than any other 
material of which we know, within its useful trans- 
mission range out to approximately 13 yu, for the 
optical elements of laboratory infrared instruments, 
for that reason. Infrared spectrometer cells made 
from it would be satisfactory for neutral or slightly 
acid organic liquids and water solutions. However, 
because of its low Knoop hardness, precautions 
should be taken to protect it when used in optical 
systems or as windows. It is also very easily at- 
tacked by alkaline and oxidizing materials. 

The verv high coefficient of expansion and low 
modulus of rupture tend to limit the thermal shock 
that can withstood. Hence, caution must be 
exercised in molding procedures, mounting design, 
and usage. The low value of Young’s 
modulus, however, does, to a certain extent, counter- 
act the effect of the high coefficient of expansion. 
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Thermal Conductivity of Beryllium Oxide From 
40° to 750° C 


David A. Ditmars and Defoe C. Ginnings 


The 


from 40° to 750° C, 


thermal conductivity of beryllium oxide has been measured by an absolute method 
The apparatus employed steady-state longitudinal heat flow along a 


rod of high-fired beryllium oxide, surrounded by a tube with matching temperature gradient 


to minimize radial heat loss. 


The estimated accuracy of the measurements is about 3 percent. 


However, the values of thermal conductivity of the ideal beryllium oxide crystal are probably 
considerably higher than the values given because of the lower density (87 percent theoret- 


ical) of the sample used, 


1. Introduction 


Not only is beryllium oxide useful as a moderator 
in the utilization of atomic energy, but it has an 
unusually high thermal conductivity, much higher 
than other nonmetals and even higher than most 
metals over a limited temperature range. At room 
temperature, its thermal conductivity is about that 
of aluminum, whereas its electrical conductivity is 
extremely low. It was the purpose of this investiga- 
tion to measure the thermal conductivity of beryl- 
lium oxide in the high-temperature range. 


2. Sample 


The beryllium oxide was originally fabricated by 
the Norton Co. by hot-pressing. A rough sample was 
taken from this material and machined, at the Oak 
Ridge National Laboratory, to a cylindrical rod 
about 0.5 in. in diameter and 6 in. long. From in- 
formation obtained on this material, together with 
that obtained from a spectrographic analysis, it 
seems likely that its impurities (other than carbon) 
were less than 0.2 percent. It is possible that carbon 
was present in larger amounts in the sample, although 
it was white, with only occasional dark inclusions. 
The sample was fired in the NBS Mineral Products 
Division at about 1,700° C and machined to the form 
of a true cylinder having a diameter of 0.4524 in. at 
room temperature and an average density of 2.62 
gem® (87 percent of single crystal). The method of 
original fabrication by hot-pressing may have caused 
a variation in density in the sample of several percent. 


3. Method and Apparatus 


The method and apparatus have been described 
briefly in technical reports {1, 2].'. The method used 
was absolute in that the results were obtained without 
comparison with another material. A longitudinal 
heat flow was used to establish a temperature gradi- 
ent in the sample. From the measured values of 
heat flow, temperature gradient, and the cross section 
of the sample, the thermal conductivity of the sample 
was calculated. Longitudinal rather than radial heat 
flow was used in order to obtain a reasonable tem- 


1 Figures in brackets indicate the literature references at the end of this paper. 
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perature gradient in a conveniently shaped sample. 
The temperature gradient on the sample was deter- 
mined by measuring the temperatures along the 
sample. The longitudinal heat-flow method, as 
applied to relatively long samples, has the inherent 
disadvantage that radial heat losses to the surround- 
ings may reduce the accuracy of the results. How- 
ever, in this case, the conductivity of the BeO was so 
high that radial heat losses did not seriously limit the 
accuracy of the results. 

A seale diagram of the essential parts of the appa- 
ratus is shown in figure 1. Measured electrical heat, 
introduced in the ‘‘sample heater” at the top of the 
sample (BeQO), flowed down the sample and _ its 
“adapter” to a heat sink. The sample heater con- 
sisted of six small helices of No. 38 Nichrome wire 
located in holes in the top of the sample. The adapter 
was used here to position the sample, as well as to 
fill in the needed length, because the apparatus was 
built to accommodate samples longer than 6 in. 
Anhydrous boric oxide was used to give good thermal 
contact between the sample and adapter, and_ be- 
tween the adapter and the sink. This compound has 
a very low vapor pressure and has excellent wetting 
properties, maintaining good thermal contact at 
temperatures far below its melting point. The sink 
was cooled with either water or air, depending on the 
temperature range, and was equipped with a heater 
and thermocouple so that it could be automatically 
kept at a constant temperature. 

The temperatures along the sample were measured 
with three thermocouples (No. 36 AWG platinum- 
platinum-rhodium) having reference junctions at 
0° C and principal junctions on the sample at the 
three levels shown in figure 1. In addition to these 
three absolute thermocouples, a differential thermo- 
couple was also used to ascertain directly the tem- 
perature difference between the upper and lower 
levels on the sample. All of the thermocouples on 
the samples were made with junctions peened into 
small holes (about 0.6 mm in diameter and depth) in 
the evlindrical surface of the sample. In order that 
the temperature gradient measured on the sample 
would correspond to the electric heat put into the 
top of the sample, precautions were taken to minimize 
radial heat loss along the sample. For this purpose, 
the sample was surrounded by a “guard tube” 





























OIC CC 
| nie of 
. ve all ‘ 
| wea r4 ad vé 
| er 
° ans 
| . e< a - 

. ih a ffere 
, an the 
. = we 4 
° 
: ple hea 
° | ilies tay gt Notas 
> AR a3 . 
+ MEE BBUy ? 
> | YEAR): 
cilia ¥oK RE ard tub 
: PT] 
. o | 
: | |, ie | 
c » wic * | 
. } | 
; r | 
a | ttorea ne 
| ° rn 
° te ° 
ICA 7: » “ xe. . 
Vté r ne »* 
1: 7 —A r 
. . R . | 
| . mk enter | 
: : | 
| es . . 
t 
° 
AZ| . nk 
tania . 
- . 
. ‘ . 
’ 
{ 2 3 
~ ae cm - : Zz “ 6 nm 
Figure |. Thermal-conductivity apparatu 
(Inconel, 0.03-in. wall), with temperatures corre- 


sponding to those in the sample. Temperatures along 
the guard tube were measured with three platinum- 
platinum-rhodium thermocouples, with principal 
junctions attached to the guard tube at levels corre- 
sponding to those of the sample thermocouples 
Three additional thermocouples, at the same levels 
on the guard tube but at different azimuths, were also 
provided. The top of the guard tube was heated by 
the “guard heater’, and its temperature was con- 
trolled automatically, using the ‘control couple’’ 
close to the heater. Another heater with a control 
thermocouple not shown in fig. 1) was provided ai 
the bottom of the guard tube, but it was found un- 
necessary to use it when the furnace temperatures 
were suitably controlled. The top portion (“‘temper- 
ing ring’) of the guard tube was made of thick nickel 
to which the electrical leads to the sample were 
thermally connected, so that it served as a temper- 
ing region for bringing the leads to the temperature 
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of the top of the sample, and thus to reduce heat 
conductior alone them. A “thermal shield”’ nickel 
and heater, placed above the sample and tempering 
ring, also served for this purpose and to prevent heat 
transfer upward through supports and insulation, 

Although the apparatus shown in figure 1 has heep 
used only for measuring the thermal conductivity of 
solids, it was designed so that it could be used also for 
high-conductivity liquids. For this application, the 
liquid would be contained in a_ thin-walled tube 
equipped with a suitable heater and a liquid expan- 
sion chamber extending up into the shield. To obtain 
the thermal conductivity of the liquid, it would be 
necessary to account for the thermal conduction ip 
the contamer tube by making another experiment 
with the liquid replaced by a powder of very low and 
known thermal conductivity. 

The whole assembly, supported by rods extending 
down to the bottom of the guard, was filled with fine 
aluminum oxide powder for insulation and enclosed 
in a 0.010-in. Inconel tube, which served to hold the 
powder. This tube, fitting closelv in the “container” 
0.035-in. Inconel tube) centered alone the axis of the 
furnace, facilitated assembly. Thermocouple and 
heater leads were brought out from this container 
through an insulating seal at the top (not shown in 
fig. 1), so that the conductivity apparatus could 
either be evacuated or filled with argon. 
experiments with the apparatus had shown the elee- 
trical insulation to fail when it was evacuated at high 
temperatures, the apparatus was filled with argon 
after evacuation at moderately high temperatures to 
outgas the aluminum oxide. The thermal-conduetiv- 
itv apparatus was maintained at the chosen tem- 
perature by the surrounding furnace, which was 
equipped with numerous taps on the heater winding 
to give the desired vertical temperature gradients 
Two automatic regulators, actuated by thermo- 
couples in the furnace, were to maintain 
constant temperatures. 


Because 


used 


4. Experimental Procedure 


Two basic types of experiment were performed for 
each measurement of thermal conductivity In the 
tvpe the 
furnace temperature was controlled the desired 
value, a known constant electric power was put into 
the sample heater, the temperatures of the guard and 
shield were adjusted to match those on the sample as 
closely as possible, and the sink temperature was 
adjusted to a constant value. After a steady state 
was obtained, several consecutive sets of observations 
of thermocouple emf and sample power were. re- 
In the second type of experiment (called 
power was put into the 
sample heater, but in other ways, the experimental 
procedure was similar to the first. The purpose of 
this experiment was to correct for errors that did not 
depend on the power transmitted through the sample. 
Errors of this tvpe are those due to differences, in 
thermocouples and those resulting from unknown 
heat leaks, which presumably were the same in both 
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types of experiments. Several other experiments 

| were made to detect other errors and to determine 
their importance. Some of these experiments are 
described later. 

Several hours were usually required to bring the 
various parts of the upparatus to the desired tem- 
peratures and to make sure that these temperatures 
were not changing significantly. The final data were 
usually obtained in a period of about 30 min, sub- 
sequent to an interval of about an hour during which 
the temperatures were observed 
Automatic thermoregculators were used to control the 


to be constant. 


temperatures of the furnace, guard, and sink, 


5. Calculation of Results and Uncertainties 


The conductivity values were calculated from the 
observed quantities by means of the equation 


QAN 
Ag(1+at,) At 
where 
i-=thermal conductivity (watts em™' deg 
(") at temperature f 
f,--average temperature of sample between 
, thermocouples 
(J=heat-flow rate (watts 
AX, =thermocouple spacing at 0° C (em). 
A, = cross-sectional area at 0° C (em*). 
a—coefficient of linear thermal expansion 
deg > 
Aft=temperature difference (deg C) between 


thermocouples 


This equation is valid for steady-state longitudinal 
heat flow over a small temperature interval. The 
determination of these factors in the conductivity 
equation, together with a consideration of their un- 
certainties, will now be discussed individually. The 
uncertainties referred to in this report are the authors’ 
estimates (based on their judgment) on the basis that 
the observed quantity would have about an equal 
chance of being within that limit as being outside it, 
and that the sign of the uncertainty is just as likely 
to be positive as negative. 


5.1. Heat-Flow Rate (Q) 
a. Electric Power 


Heat was generated by direct current in the sample 
heater, and power was measured in a conventional 
manner, using a potentiometer in conjunction with a 
high-resistance volt box to measure the potential 
drop across the heater, and a standard resistor in a 
current lead to measure the current. The errors in 
these electrical measurements were almost negligible. 
In measuring the potential drop across the heater, 
the potential terminals were located to evaluate 
properly the heat that went to the sample. Because 
the sample heater was made with very high resistance 
(relative to that of the heater leads), the uncertainty 
in the location of the potential terminals resulted in 
only about 0.1-percent uncertainty in the measured 
thermal conductiv ity. 


b. Heat Flow From Sample to the Tempering Ring 


In addition to the electric-power input to the top 
of the sample, it is necessary to consider the heat 
transfer between the sample and its surroundings. 
Because of the excess temperature of the sample 
heater, heat flowed along the heater leads to the 
tempering ring; the uncertainty in this heat flow 
(taken to be 50 percent of the correction) averaged 
about 0.24 percent of the total heat flow in the sam- 
ple. In addition to this heat flow, there was the heat 
flow between the top part of the sample and the 
tempering ring. The evaluation of this was difficult 
because of the configuration and the temperature 
distribution on the top part of the sample where the 
sample heater was located. Using differential ther- 
mocouples, observations were made of the tempera- 
ture difference between the isothermal tempering 
ring and a point on the top part of the sample. The 
location of this point was determined by calculation 
so that the net heat flow from the top of the sample 
to the tempering ring would be proportional to the 
emf of the differential thermocouple. The height of 
the top of the sample relative to the guard was made 
so that the bottom of the sample heater was at the 
seme level as the bottom of the tempering ring. In 
the actual experiments, the tempering-ring tempera- 
ture depended on the power in the guard heater and 
was indicated by the differential thermocouple to be 
2 or 3 deg higher than the temperature of the sample, 
making necessary a correction for the resulting heat 
flow. It was found conveaient to evaluate this cor- 
rection experimentally by making two conductivity 
experiments in which only the tempering-ring tem- 
perature was changed and the furnace temperature 
adjusted to maintain the match between the guard 
and sample thermocouples. Using the resulting 
change in sample gradient, the data were corrected 
to correspond to no difference in temperature as 
indicated by the differential thermocouple. How- 
ever, there still remained an uncertainty in heat flow 
due to some uncertainty in the proper location of the 
thermocouple junction on the sample head. It was 
estimated that the uncertainty in the heat ‘flow be- 
tween the sample and tempering ring, excluding heat 
flow along leads mentioned above, resulted in an 
uncertainty in measured conductivity averaging 
about 0.5 percent. 


c. Heat Flow Down the Insulating Powder 


Even when the guard temperatures matched 
sample temperatures, some of the heat input to 
the sample necessarily went to maintain some of 
the longitudinal heat flow in the insulating powder 
between the sample and guard. The conductivity 
of the aluminum oxide powder with argon gas was 
determined approximately by a few experiments as 
it was used in the apparatus. The temperature 
distribution in the insulation had been previously 
estimated with a resistance analog computer, setting 
up boundary conditions corresponding to the config- 
uration and assumed temperatures of the sample 
and guard. It was estimated from these results 
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that when the thermocouple on the top part of the 
sample indicated the same temperature as the guard 
ring, the sample heater contributed only 16 percent 
of the total longitudinal heat flow in the powder, 
the remainder of the heat being furnished by the 
guard. Under this condition, about 0.2 percent of 
the heat of the sample heater flowed down through 


the insulating powder. The uncertainty of this 
correction was estimated to give less than 0.1- 


percent uncertainty in the measured thermal con- 


ductivity. 
d. Heat Flow Between Sample and Guard 


Because it was found impractical to match the 
sample and guard temperatures exactly during all 
conductivity experiments, experiments were made 
that permitted calculation of corrections for imper- 
fect matching. For each guard thermocouple, two 
thermal-conductivity experiments were made, vary- 
ing only the difference between that guard thermo- 
couple and the corresponding thermocouple on the 
sample. From the resulting change in temperature 
gradient on the sample in these two experiments, it 
was possible to estimate a correction for small 
differences in matching guard and sample thermo- 
couples It was calculated that there was an uncer- 


tainty of about 50 percent in correcting for heat 
flow between the sample and the guard. This 
uncertainty resulted in an average uncertainty in 


measured conductivity of about 0.3 percent. 


e. Heat Flow Into Heat Capacity of Sample 

If the temperature of the sample were changing 
with time, some of the heat input would to 
produce this change, and the temperature gradient 
on the sample would not correspond to the heat 
input at the top of the sample. In all experiments, 
the rate of temperature change was less than 0.8 deg 
C/hr, corresponding to an effect of 0.5 percent in 
the calculated value of conductivity. The average 
uncertainty in the correction for this was negligible. 


- 
vO 


5.2. Thermocouple Spacing (1 X,) 


The distance between the principal junctions of 
the upper and middle sample thermocouples was 
4.97 em at 0° C, whereas the corresponding distance 
for the middle and lower thermocouples was 5.05 
em; this gives 10.02 em for the distance between 
the extreme absolute thermocouples. The distance 
between junctions on the differential thermocouple 
was 10.01 em. These distances, taken as the lengths 
between centers of the thermocouple holes, were 
measured to better than 0.01 em with a traveling 
microscope, but because the thermocouples were 
peene “l into holes 0.06 em in diameter, the possibility 
of nonuniform thermal contact makes a tolerance of 
0.03 em appear more realistic. This tolerance 
responds to 0.5-percent uncertainty in the 10-em 
spacing between the sets of thermocouples used in 
the conductivity calculations. The effect of thermal 
expansion on both thermocouple spacing and cross- 
sectional area is lumped into the correction (1+ af, 
which is described later. 


cor- 
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Cross-Sectional Area (Ay) 


The sample was ground to have a uniform diameter 
of 0.4524+0.0003 in., corresponding to a cross. 
sectional area of 1.038 em*®. The uncertainty jn 
this area was estimated to be less than 0.1 percent. 


5.4. Thermal-Expansion Correction (1+ at,) 


Thermal-expansion changes the thermocouple 


spacing by the factor (1+ af,), and the cross-sectional 
area by the factor (1+ af,)*, resulting in the (1-4 at,) 
term given in the conductivity equation. The 


coefficient of linear thermal expansion (a) has been 
determined by White and Schremp [3]. At the high- 


est temperatures (747° C) of the conductivity ex- 
periments, the correction for expansion amounted 
to about 0.6 percent, with negligible uncertainty 


in the measured conductivity. 


5.5. Temperature Difference (A?) 


The accurate measurement of the temperature 


difference on the sample was difficult, requiring a 
number of tests to eliminate or evaluate certain 
errors. As described in section 3, two different 


thermocouple svstems were used to measure the 
temperature difference over the 10-cm length on the 
sample. The two independent thermocouple sys- 
tems served to check on each other, usually agreeing 
on the measured temperature difference to better 
than 1 percent. 

All temperatures were measured with platinum 
platinum-rhodium thermocouples of No. AWG 
wires. A sample thermocouple was calibrated at 
several points between 0° and 1,000° C in the Pvrom- 
etry Laboratory the National Bureau of Stand- 
ards. No significant difference was he- 
tween the thermoelectric power of this sample and 
that given in the standard tables [4]. Even though 
the thermocouples were all made from wire off the 
same spools, it was possible that they had _ slightly 
different thermoelectric powers. Although these 
differences might not be serious for measurement of 
absolute temperatures, they could be significant for 
measurement of small temperature differences at 
high temperatures. These differences were auto- 
matically accounted for by the calibration experi- 
ments mentioned previously, in which no power was 
put into the sample heater. These calibration ex- 
periments gave differences in thermoc ouple readings, 
which increased regularly up to 470° C, and 
amounted to as much as 7 uv at this temperature. 
At temperatures approaching 750° C, the differences 
became both larger and more irregular, so that after 
each conductivity experiment, a calibration experi- 
ment was made to evaluate the differences. 

In addition to the differences described above, the 
evaluation of Aft was uncertain because of the occa- 
sional effect of high humidity on the potentiometer 
used to measure the thermocouple emf. This effect 
was observed as a reading on the potentiometer, 
even When the potential across its terminals was 
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gro. When the value of At was observed, using 
the differential thermocouple, the value was subject 
to the full potentiometer uncertainty (about 2 yv) 
hecause only one reading was involved. When the 
value of At was determined by using the two absolute 
thermocouples, most of this potentiometer uncer- 
tainty was reduced (to about 0.5 wv) because the 
value of At was obtained from a difference of two 
readings. It is for this reason that the uncertainties 
in the value of Af, using the differential thermo- 
couple, were larger than the uncertainties when 
using the absolute thermocouples, averaging about 
1.5 percent as compared to 0.4 percent. That the 
two thermocouple systems usually agreed to better 
than 1 percent is evidence that the error due to the 
humidity effect on the potentiometer was _ not 


excessive. 
6. Results 


The results of the individual thermal-conductivity 
measurements on bervllium oxide are given in table 
1. Values of observed conductivity (4) are given as 
determined by using each of the thermocouple sys- 
tems (absolute and differential) at the average tem- 
perature (f,) of that portion of the sample measured 
by the thermocouples. At the lower temperatures, 
where the thermal conductivity of the beryllium 
oxide changes rapidly with temperature, small cor- 
rections to conductivity were made for the curvature 
of the conductivity-temperature function. In this 
table, the quantities given are corrected for all known 


errors. In the previous discussion, each uncertainty 
has been estimated by the authors on the basis that 
the observed quantity would have an equal chance 
of being within that limit as being outside that 
limit. ‘These uncertainties have been combined 
(square root of the sum of the squares) and arbitrar- 
ily increased by over a factor of 2 to give more 
realistic values of estimated error listed in table 1. 

Table 1 indicates that the results, using the 
absolute thermocouples, seem to be reliable to 
about 2 percent. The experiments under 60° C are 
not as accurate as the other experiments. This is 
due to the smaller temperature drop in the sample, 
first because the limitations of the heat sink made 
it necessary to use lower power, and second, because 
the thermal conductivity of the beryllium oxide was 
so high in this low-temperature range. The larger 
error estimated by using the differential thermo- 
couple in this low-temperature region is due to the 
humidity trouble mentioned previously. No results 
are given for measurements with absolute thermo- 
couples above 500° C because of failure of their 
electrical insulation. 

A smooth function of thermal conductivity was 
obtained graphically from the observed values in 
table 1, giving greatest weight to those values hav- 
ing the smallest estimated errors. Table 2 gives 
smoothed values of the conductivity at even tem- 
peratures as obtained from the graph. Figure 2 
gives the deviations of the results (obtained with the 
two different thermocouple systems) from the smooth 


TABLE 1 Ex pe rimental results 
Absolute thermocouples Differential thermocouples 
) 4 Pow 
At Observed Estimated At Observed | Estimated 
error error * 
( / ( m-° ¢ ( wicem-~ ay A 

Aug. 2 is. 2 0. 8960 3.9 2.19 4 4.11 2.10 19.2 
Aug. 2¢ HH). 2 1. 4452 6.84 2.04 3.2 03 1. OS 10.7 
Aug. 2 2.6 1. SS04 9. 41 1.04 2. ¢ +61 1. 90 7.9 
Aug. 2 9.7 2. 3620 12. 15 1.S87¥ 2.1 12. 33 1. 849 6.2 
Aug. 24 85.38 $. 3179 19. 16 1. 674 lt 19. 28 1. 661 4.0 
2) s SH, 2 3. SOOT 21. 88 1.679 1.7 22.09 1. 660 3.6 
Sept. 20 aH. u 3. 6073 20. 40 1. 668 > F 21.07 1. 652 3.4 
Sept. 2 s ; 3. 6084 20). 86 1.672 1s 21. 08 1. 652 as 
Sept. 3 87.3 3. 6138 21.01 1. 662 1.7 21.17 1. 647 3.7 
Au ; t 2 5003 15. 20 1. 647 2.0 15. 10 1. 655 14 
Se} l 23.8 2.4474 16. 26 1. 454 Le 16. 33 1. 446 5.4 
Sept. 1 3.0 2. 4804 18. 27 1. 316 1.7 18. 36 1. 308 4.1 
sept. 2 3.4 2.479 18. 27 1. 311 1.8 18. 37 1. 302 4.2 
Oct. 11 202.0 3. Y603 33. 80 1. 131 2.0 33. 65 1. 135 2.6 
Aug. 23 202. 1 3.874 $3. 29 1.124 1.6 33. 10 1. 128 2. 4 
Au Is 241 $. 7312 36. 20 0. G95 & 35. 87 1. 002 23 
Au 1v 241.4 3. 7396 36. 51 YSU 7a 36.19 0. 996 = 2 
Aug. 20 251.4 0. 8739 &. 05 W42 2.2 SAS 952 6.5 
Se} 2 287.3 2. 3441 5. W3 S72 1.8 25. 54 SS4 2.8 
Se} li 287.7 2. 3378 25. 87 S72 1.9 25. 71 S76 2.9 
Oct 379. 5 $. S169 18. 12 704 1.5 47.77 709 2.0 
Oct $80.4 3. 5378 8. 64 701 1.Y 17. 60 715 2.3 
n 12 $39.2 3. 3054 51.83 614 1.2 1. 82 614 1.7 
0 Is 17.2 2. 9239 3.63 524 2.7 
On 19 7 7 2 5664 1. 54 478 3.5 
Oct. A 646.7 > 7048 59. 04 433 2.4 
N 2 747.9 ? 6370 65, 84 384 41.3 

Es " ! iuthors’ estimate, considering only the various uncertainties mentioned in the text, 

rt u \ ist 20 and 23 represent averages of 2 experiments on ¢ wh day 
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function. At 50° C and below, the observed values 
deviate markedly from the smooth function, giving 
values at 40° C that are 4 percent different. These 
deviations are probably due to the very small tem- 
perature difference of 4° C in the sample, so that 
errors in the measurement of this temperature differ- 
ence have greater influence on the result. Figure 2 
also shows the results at 251° C to be about 3 per- 
cent lower than the other results. It seems probable 
that this departure is also due to the lower tem- 
perature difference on the sample resulting from the 
lower power; the power here was only one-fourth the 
power in the other experiments in this temperature 
range. If it is assumed that there existed an un- 
known constant absolute error in either heat flow or 
temperature difference, the deviations of about 3 
percent in the low-power experiments would indi- 
cate that the experiments with the higher power 
might be in error by about 0.8 percent 
quently, the authors believe that the over-all ac- 
curacy of the results is more likely to be about 3 
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percent instead of the 2 percent indicated by the 
estimated errors listed in table 1 
TABLE 2. Thermal conductivity of BeO (densil 62 q 
Sn , 
ft ‘ f f ( 
l iu) ‘ 
4 1 7 im ~ 
iM s| uN 4] 
| l us 4 {US 
Ji 1. 131 un Wi2 
™y 0. “72 “ 14 
wn a4 iw) iw 
“ 74 is 4 
3 lites al cel pial ‘ail biti ttle 
4 
} 
2 
* 
: | 
z 
- «2 coo 
~ Oo 0 
< as ’ 
- } 
Mad 
ao 
/ 
! j 
> C + 
2 & 
J |! 
| ; 1 
4 
4 é 8 


Frat RE 2 be ations of ber yll im onide thermal-cond we 


“s N BS smoothed dat 
N BS (absolute ce ! 


98 


No attempt has been made to correct the thermal. 
conductivity values for the NBS sample to corre. 
spond to zero porosity. Because the density of this 
sample was only 2.62 g/em* compared to about 3.9 
for the ideal crystal, the thermal conductivity of the 
crystal should be significantly higher than the values 
given in table 2. From the investigation of Frane} 
and Kingery [5], it would appear that the condue- 
tivity of the ideal crystal would be about 15 percent 
higher. However, the measurements of Powell [6] 
on beryllium oxide specimens (densities 1.85 to 2.89 
¢/em’) would indicate a much larger correction. The 
authors feel that the correction for porosity is un- 
certain and that there are other factors beside poros- 
itv that also should be accounted for. One of these 
bonding of the individual 
particles by the firing process. 


factors is the degree of 


7. Comparison With Other Results 


Figure 3 gives a comparison of the results of the 
NBS measurements on BeO with the 
measurements at other laboratories on other samples 
At the lower temperatures, the agreement with 
Scholes [7] is probably as good as the plivsical states 
of the two samples permit. Scholes used a sample 
having a density of about 2.97 g¢ ¢m* as compared 
with the NBS sample having a density of 2.62. At 
higher temperatures, the Francl and 
Kingery [5] are consistently higher than the NBS 
results They used a sample having a density of 
2.86 g¢/em*, but this would probably 
only a small part of the difference. The results of 
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A Self-Balancing Direct-Current Bridge for Accurate 
Bolometric Power Measurements 


Glenn F. Engen 


Until recently, the most accurate microwave power measurements of the bolometric 


type have required the use of a manual d-e bridge. 


A self-balancing d-c bridge has been 


developed that preserves the inherent accuracy of the manual bridge, extends the dynamic 


range of operation, and greatly simplifies the operating procedure. 


A general description 


of the equipment and operating techniques is given, followed by a comprehensive survey of 
the sources of error accompanying the method and the accuracy achieved. 


1. Introduction 


A large amount of the research effort expended in 
the field of low level microwave power measurement 


in the past few years has been directed toward a | 


determination of bolometer-mount. efficiencies, or 
evaluation of the validity of the r-f-d-e substitution 
principle.'?"* Until recently the only accessory 
instrument of sufficient accuracy to utilize the results 
of these investigations has been a manually operated 
de Wheatstone bolometer bridge. The use of this 
device is, however, both tedious and time con- 
suming; the dynamic range over which the desired 
accuracy is achieved is rather small, and the danger 
of barretter burnout due to a misstep on the part of 
the operator is ever present. 

A self-balancing d-c bolometer bridge has been 
developed at the NBS Boulder Laboratories, which 
to a large measure relieves these problems, while 
preserving the accuracy and extending the dynamic 
range of operation. ‘To date, the technique has been 
employed primarily with barretter-type elements, 
although on the basis of preliminary results, the 
technique — should equally useful with 
thermistors. 


proy c 


2. Manual Bridge 


One of the simplest forms of the manual bridge 
previously emploved at the National Bureau of 


Standards is shown in figure 1. The r-f power ° is 
given by 
/ ° . 
P., 17 — 2). (1) 
; 
where 7, and iy are the total currents required to 
halance the bridge without and with r-f power 


present, respectively. A precision potentiometer- 
standard resistor combination is used to measure 


) and 7. If it is assumed that this measurement 

R. W. Beatty and Frank Reggia, An improved method of measuring efficien 
cies of ultra-high-frequency an t bolometer mounts, J. Research 
NBS 54, 321 (1955) RP2594 

A.C. M wpherson ind D. M Ker { microwave microcalorimeter, Rey 
Sei. Instr. 26, 27 (1955 

H. J. Carlin and Max Sucher, Accuracy of bolometric power measurements, 
Proc. Inst. Radio Engrs. 4@, 1042 (Sept. 1952 

‘D. M. Kerns, Determination of efficiency of microwave bolometer mounts 
from impedance data, J. Research N BS 42, 579 (1949) RP1995 

More correctly, ‘‘the retracted d-c powe! 


contains a random error of +0.005 percent, it can 
be shown that for a typical bolometer operating at 
200 ohms and requiring 8.5 ma of bias (7; =17 ma) 
the error in measuring P,,; may be as large as 0.02, 
0.3, and 3 percent at the 10, 1, and 0.1 mw levels, re- 
spectively, if all other sources of error are neglected. 
The small difference between 7, and 7 at low levels 
accounts for the deterioration in accuracy. <A 
variety of alternative techniques have been proposed 
and employed with some degree of success, but it is 
not within the scope of this paper to examine the 
manual technique in detail or to establish the 
maximum accuracy obtainable. In practice, this is 
often determined by the stability requirements 
imposed on the generator and measuring system by 
the time-consuming nature of the technique, rather 
than by more fundamental limitations. 
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FiGURE | Manual direct-current bolometer bridge. 


3. Self-Balancing Bridge 


The self-balancing d-c bolometer bridge employs a 
d-c amplifier in a feedback loop to balance the 
bolometer bridge. The amplifier is connected in 
such a manner that the bridge unbalance signal is 
amplified to provide the bridge current (fig. 2). The 
use of this self-balancing circuit eliminates many of 
the time-consuming features of the technique, but 
the power-measurement problem is essentially the 
same. This problem, at low power levels, is basically 
one of trying to accurately measure a small change 
in a comparatively large d-c power. 
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where Ai ) 


The difficult factor to measure accurately is A 
A constant-current generator of high stabilitv: has 
been developed for use in conjunction with a differ- 
ential current-measuring technique, which enables a 
direct and accurate determination of Aj. The circuit 
is illustrated in figure 3 and operates as follows 
First, with no microwave power present in_ the 
bolometer, 7, is determined by ammeter A,, and the 
output of the constant-current generator is adjusted 
with the aid of ammeter A, to equal 7, (ammeter 1 
reads zero Then, with microwave power applied, 
A, will read Ai directly, and the power may be 
determined from eq (5 In practice, a potenti- 


ometer and standard resistor are emploved in place | 
of A, and A, to achieve high accuracy. 
Alternatively, a voltmeter, 
will read (79/4 ii is , and the power Is given by the } 
product of ! and Abo. : ‘hese operations Mla also he 
combined in a single instrument of the dvnamom. 
eter type to achieve a direct reading of the power 
The degree of success obtained with this technique 


, connected as shown 


is, of course, contingent upon the degree of stability 
achieved in the constant-current generator. In this 
case the change in generator output for a complete 
cutoff of bridge current is less than | part in 10% 
resulting in i neghgible error from this source 


4. Performance 


It is shown in a following section that the over-all 
accuracy In power measurement achieved, using a 
potentiometer to measure the current, is about 0.1 
percent, whereas substitution of a suitable dyna- 
mometer to achieve direct reading is expected to vield 
The instrument 
thus provides an increase in accuracy of 10 to 50 
times that provided by the commercially avaiablel 
audiobridges. Other advantages over the tvl pica 
audiobridge include: (1) Comparatively simple tech- 
niques are avatlable for establishing, verifving, and 
monitoring the specified accuracy. (For example, 
the bridge balance may be monitored by a suitable 
calvanometer. 2) The performance is, to a large 
degree, independent of the nominal deterioriation of 
the vacuum tubes and other electronic components, 
3) The elimination of the audio bias power avoids 
its interaction with the pulse-repetition frequeney in 
the measurement of pulse power 

At low power levels the r-f power is approxi- 
mately proportional to Ay (eq 3), and the low-level 
performance may be displaved by recording this cur- 
rent with the r-f power alternately on and off. The 
system response to a l-uw signal is given in figure 4. 

Barring failure of the electronic components, the 
self-balancing operation eliminates the danger of 
bolometer burnout due to excessive current, 
whereas protection against an r-f overload may also 


an accuracy of about 0.3 percent 


bias 
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hiGgurReE 5 Prototype 


be provided through use of a suitable relay to turn 
off the r-f source when the power input to the 
bolometer exceeds a prescribed amount. 

The d-c amplifier emploved in this application is a 
commercially available instrument of exceptional 
gain, zero stability, and phase response, which was 
modified by the addition of a cathode follower to 
provide the required output current. The current 
amplification of the amplifier, as thus modified, is 
approximately — 100,000 Compensation for the 
thermal lag of the bolometer is provided in the feed- 
back loop by means of a condenser, as indicated in 
figure 2, but the achievement of stable operation 
with this large amplification is primarily a tribute to 
the amplifier phase response. The constant-current 
generator emplovs a second amplifier of the same 
tvpe, which, with modification, maintains a constant- 
output current by reference against a mercury cell. 
A prototvpe model of this equipment is shown in 
figure 5 


5. Analysis of Errvurs 


The practical problems involved in making an ac- 
curate determination of the power, using eq (1) to 
have been noted. But use of these formulas 
implies that the bridge is exactiv balanced, the bridge 
arms are equal to ro, ete. Li is shown in the appendix 
that a more complete expression for the power is 


(.4 


mode lol self-hbalanc ing bridge, 


front and rear views. 


where 


current amplification of 


a { don ino) (igs - 
amplifier 
“—fractional deviation of the lower bridge arm 

from the value ro 
r,=resistance of bolometer leads 
e=thermal (or other) emf in amplifier input 
total bridge current 
(, amplifier input current 
amplifier input resistance 
y=“ohms per watt” bolometer coefficient. 


For a typical bolometer operating in the 0 to 10- 
mw region, the following values may be assigned: 


200 ohms 
y=—4.5 ohms/mw 
i», == 17 ma 

ino >S ma 


lo 


Substituting these values in the above expression 


a) 


thi 
ts. ° ) 

9 ’ 
a - 


where the last three terms in the bracket reflect the 
inability of the d-c amplifier to balance the bridge 


vields the result 


' lo,: - = ee l 
Puy M1 mk +Z—- +-5.6¢e-4 LS(] + @)( -— 


— 8094 


exactly. 
The characteristics of the d-c amplifier employed 
in this application are such that the following values 
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obtain: 


g=1/4 (amplifier input resistance 50 
ohms) 
«>10 (current amplification) 


(amplifier zero stability). 
e|<10-° J 


Substituting these values in the last three 
gives a value for their sum of about 0.0002. 

The bolometer lead resistance can be measured 
with little effort and the correction applied, so the 
error from this source is negligible, whereas the value 
of |Z! can ordinarily be held to 0.0001 or less. 

Finally, it is possible to measure current by means 
of a precision potentiometer-standard resistor com- 
bination to an accuracy of a few parts in 10*; thus 
the first factor in eq (4) may be determined by the 
technique discussed earlier to an accuracy of about 
0.06 percent. 

In summary, the total limits of error are as follows: 


terms 


Deviation of bridge arms from 
value ro 0. O1 percent 
Failure of bridge to be exactly 


balanced 
Measurement of (4; + th 
Measurement of i»; — tp 
Bolometer lead resistance 
Amplifier and instability 
(experimentally observed 
Total 


. 02 percent. 
. 03 percent 
03 percent 
. OO percent 
noise 
0.1 ww 


. O9 percent +0.1 uw. 


6. Appendix 


From figure 6 the power may be expressed in 


the form 


r 


Pres 2) LAN YZ ,2,,€,77) (5) 


(t51—2s 


where Y, Y, and Z are the dev lations of each of the 


bridge arms from the value 7, 7, and e are the 
detector current and thermal (or other) emf, re- 
spectively, and r, the resistance of the bolometer 


leads. 

Assuming a knowledge of the bolometer resistance 
law, an exact expression for the second term could 
be derived, but the labor required suggests the fol- 
lowing approximate treatment. By hypothesis, the 
variables Y r, and f(X . r,) are small, and 
S(X ... r)) may, toa good approximation, be written 
asjthe sum of the first-order corrections for each of 
these quantities, where the correction for each of 
the variables is determined with the other variables 
equal to zero. 

Application of this technique to the circuit in 
figure 6 vields the following first-order expressions 
for the d-c power dissipated in the bolometer, and 
the bolometer resistance 
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The failure of the expression for power to have a 

first-order dependence on \ or Y is explained as | 
follows: A prescribed, small value for X or Y will 

cause a first-order variation in bolometer resistance 

r, as indicated in eq (7), but also produces an unequal } 
division of current in the two sides of the bridge. 

These two sources of error cancel to a first approxi- 

mation. The only general requirement on the upper 

bridge arms is that they be equal, and if they havea 

nominal value of ro, it can be shown that they may 

differ from each other by as much as 2 percent } 
and produce an error in the power measurement 
of only 0.01 percent. 


FIGURE 6 


rrors 


0 'nlo 


ro=rp I a y | ; a 4e 
r 


At a constant ambient temperature and assuming j 


equivalence of d-c and r-f heating,® the power, P, 
dissipated in the bolometer will be a function of its 


resistance, 7: 


P=xf (r). 


Because the bolometer is nearly a square-law 
detector, one may in the neighborhood of an ar- 
bitrarily selected operating resistance ry, write 


oP 
or 


/ lo), 


P= P,+ 


where P. is the power required to bring the bolometer 
to the resistance 7) at the given ambient temperature, 
P is the sum of the d-e and r-f power, P=Py.+ Pry, 
and OP /dr is the reciprocial of the “ohms-per-watt” 


issumed to be negligibl this discussion 


footnotes 2 and 3 


Che r-f-—d-e substitution error is 
For a treatment of this proble m, se¢ 
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bolometer coefficient, i. e., (OP/dr)=(1/¥). Then is subject to change with the application of r-f power. 
Substituting eq (6) and (7) in eq (8) gives the desired 


| result 
P. +P pe Pote () Io). 
) ! . . 9) X 
lo > +9 “9 , ry at < 
| Pye? (i i) | 1+ 2-24 — (1+ 
P, may be determined by measurements in the ; Po POT foe Ybor boa 
absence of r-t power: : 16(14 ) # 24 3 
+ ; : . ( —s rh ; + Lge ’ 
| Ytn2\ 51 7 Une) \@ Un) (lpi tT too) \ @ 
Pom Par (ni-n 
i 
where 
; ° ° : . 
where the subscript 1 has been used to indicate the ate 
value of quantities measured in the absence of r-f bei — Ugo 
power, Whereas the subscript 2 will indicate values 
with r-f present. 
The formula for r-f power is then in , 
The author extends his thanks to R. W. Beatty 
: : and John E. Lauer for their valuable criticisms and 
} P..=P)—Paet— (r2—19) =P aaa —Paee——(r,— rs). (8) | Suggestions in the preparation of this paper, and to 
rors Y 7 | Rav R. Rumfelt for his assistance in the design and 
construction of this equipment. 
In the given application, only 7, of the variables 
6 X...vr, has a functional dependence upon 7, or Bou.pEerR, Coro., November 28, 1956. 
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Some Observations on Hydrated Monocalcium Aluminate 


and Monostrontium Aluminate 


Elmer T. Carlson 


Monocalcium aluminate hydrate, 


was prepared by precipitation from calcium aluminate solutions at 1° C, 
aluminous cement. It 


of pastes of monocalcium aluminate or 


probably having the composition CaQ-AlO,-10H,O, 


hydration 
as hexagonal 


and by 
was obtained 


prisms, very weakly birefringent, with mean index 1.471. Three molecules of water were 
expelled by drying over CaCh, without apparent change in crystal structure; however, the 
structure was destroved by heating to 175 Monostrontium aluminate hydrate was 
obtained by analogous methods in the form of minute needles or prisms with mean index 
1.478. The maximum degree of hydration was not successfully determined, as the com- 
pound underwent decomposition to 3S5rO-Al,0,-6H,O and gibbsite. X-ray diffraction 
powder patterns indicate a close similarity in structure between the monocalcium and 


monostrontium aluminate hydrates. 


1. Introduction 


The existence of a monocalcium aluminate hvdrate 
of some uncer- 


has been, for many vears, a matter 
taintv among investigators in the field of cement 
hydration. Over 20 years ago, Assarsson [1] ' re- 


ported the a ation of a compound having the 
composition CaQ-Al,O.-10H,O It was obtained 
the form of a “gel” by the action of water on an- 
hydrous calcium aluminates and aluminous cements. 
Although Assarsson gave an X-ray diffraction powder 
pattern for the material, this in itself was not con- 
elusive evidence of the existence of a monocaleium 
aluminate hydrate. A number other hydrated 
calcium aluminates give somewhat similar patterns, 
and the gel might conceivably have been a mixture 
of two or more phases. 


of 


More recently, Longuet [2] reported the prepara- 
tion of a compound for which he indicated the 
formula CaO-AlLO,-7H,O as most probabie. He 
also gave an X-ray diffraction pattern, differing in 


several respects from that given by Assarsson. 

It has been pointed out in a previous paper [3] 
that the two patterns probably can be reconciled; 
that patterns agreeing in general with that of 
Longuet have been obtained at the Bureau on 


lhvdrated pastes of aluminous cements and of mono- 


also, 


calcium aluminate. Because of the probable im- 
portance of this compound in the setting and 
hardening of aluminous cements, some additional 


work was undertaken in an attempt to obtain the 
hvdrate in better erystallized form. 

While this study was in progress, 
important investigations bearing on the subject 
were reported in the literature. In the revised 
edition of “The Chemistry of Cement and Concrete” 
[4], Lea gives an X-ray diffraction pattern obtained 
by H. G. Midgley for CaQO-Al,O,-10H,0. The 
pattern agrees fairly well with the more abbreviated 
patterns given by both Assarsson and Longuet. 
Farran [5] developed a novel technique for the study 
of crystal formation at the interface between cement 


two additional 


Figures in brackets indicate the literature references at the end of this paper 


a gh aggregate, and by this means obtained crystals 
of a monocalcium aluminate hydrate large enough 


iy a determination of optical properties and shape. 
The crystals grew from a paste of aluminous cement, 
and the hydration product contained small amounts 
of one or more other phases. Farran also published 
an X-ray pattern for his preparation. 

The existence of a monocalcium aluminate hvdrate 
thus seems established. Because of the close rela- 
tionship between calcium and strontium, it was con- 
sidered of interest to determine whether a mono- 
strontium aluminate hydrate also can be prepared. 


2. Apparatus and Procedure 


From the work of Assarsson [1] it appeared likely 
that formation of the monocalcium aluminate hy- 
drate would be favored by low temperature; conse- 
quently, in the present study, most of the reaction 
mixtures were held near the freezing point. A 
household refrigerator, equipped with a special 
control unit and a blower, was used as a cold cabinet. 
The temperature ranged between 0° and 2° ¢ 

Supersaturated solutions were prepared by acltie 
synthetic anhydrous monocalcium aluminate with 
water, filtering, and diluting the filtrate with dis- 
tilled water or lime water as required. The solu- 
tions were kept in polyethylene bottles. A similar 
procedure was used in the study of the strontium 
analog. 

The solutions were allowed to stand in the cold 
cabinet until precipitation was pelle complete, 
which in one case required several months. The 
precipitate was then separated by filtration with the 
aid of an aspirator and a suction train to mini- 
mize contamination by atmospheric carbon dioxide. 
Washing was omitted, because the solutions were 
extremely dilute, and the resulting contamination 
was negligible. 

In addition to the experiments based on precipita- 
tion from supersaturated solutions, a few were con- 
ducted with pastes made from the anhydrous 
aluminates with small amounts of water. 
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3. Results and Discussion 
3.1. Monocalcium Aluminate Hydrate 


to three preparations of mono- 
calcium aluminate hydrate are given in table 1. 
The precipitate in No. 1 was so finely divided that 
it formed a flocculent suspension that failed to settle 
out completely. Under the light microscope it ap- 
peared amorphous, but its crystalline character was 
revealed by the electron microscope, as well by 
X-ray diffraction. Figure 1 is an electron micro- 
graph of this preparation. Other areas in the field 
showed some poorly defined matted fibers, possibly 
of amorphous hydrated alumina, and a very few 
hexagonal plates, much larger in size, were also 
noted. Neither of these phases is visible in figure 1. 
The X-ray diffraction pattern was similar to those 
reported by other investigators [1, 2, 4, 5] for mono- 
calcium aluminate hydrate. No additional diffrac- 
tion lines were observed. 

The precipitate from solution 2 formed aggregates 
of fine crystals, which settled rather rapidly after 
the mixture was shaken. When first examined 
microscopically, 11 days after the reaction mixture 


Data pertaining 


was prepared, the precipitate was well formed. 
There was no visible change in it at 74 days, when 
it was filtered off. The crystals apparently were 
isotropic, with a refractive index about 1.47. The 


index was observed to vary with the degree of drying. 
Under the electron microscope the preparation ap- 
pe ared to consist of aggregates of prismatic crystals, 
much larger than those of preparation 1. Some of 
the prisms showed terminal pyramids. <A few hexag- 
onal plates were also observed. The X-ray diffrac- 
tion pattern was similar to that of the first prepara- 
tion, but the intensities were stronger. A very 
weak line corresponding to the strongest reflection 
of 2CaQ- Al,O,-8H,O was also observed. 

Preparation 3 was initially much more dilute than 
the others. Gibbsite was added to the solution, in 
the amount of 1 g for 1 liter, to hasten the precipita- 
tion, because it was supposed at the time that gibb- 
site would be the stable phase in this area. Precipi- 
tation was negligible during the first month, but after 
6 months a new solid phase was observed, mixed 
with the gibbsite. It consisted of short hexagonal 
prisms with terminal pyramids (see fig. 2), occurring 


TABLE | » ALO 


-1OH () fron 


singly or clusters. Also, in many cases, it ap. 
peared as a chainlike aggregate of crystals, which 
appeared to have grown around a long, slendep 
needle. When vie wed in contact with the solution 
in which they were grown, the prisms showed ex. 
tremely low birefringence. The index of refraction 
of the prisms agreed with that of preparation 2. 


and the X-ray diffraction patterns showed that they 
were the same phase. The nature of the needles jg 
unknown; they had the same index of refraction ag 
the prisms, but showed slightly higher bire fringence 
A few of the needles ma \ be observed in figure 2. The 
prismatic crystals are similar in appearance to thoge 
described by Farran, except that Farran’s crystals 


are, 1 comparison, much shorter in the direction of | 


the prism axis. 
In addition to gibbsite and monocaleium alumi- 
nate hydrate, another phase was present in minor 


amount in preparation 3. It consisted of large 
hexagonal plates, and was at first assumed to be 
2CaQ-Al,O;-8H,0. The X-ray diffraction pattern, 


however, showed no trace of the dicalcium compound, 
It did contain three weak lines that could be attri- 
buted to tetracalcium aluminate hydrate, along with 
strong patterns of monocalcium aluminate hydrate 
and gibbsite. The appearance of tetracalcium 
aluminate in this region of the system, however, 
would be contrary to expectations based on earlier 
work at 21° C [6]. 

The molar oxide ratio, CaO/AlLOs, in preparation 1, 
determined by chemical analysis, was 1.04. In 
preparation 2 it was 1.01. Greater weight is attached 
to the latter value, this preparation was 
microscopically almost homogeneous. The only 
observable impurity was dicalcium aluminate hy- 
drate. present in amount estimated at than 
5 percent. Preparation 3 was not analyzed, as there 
were three solid phases present; also, the amount of 
material available was very small. Additional 


because 


less 


gibbsite may have been precipitated (as was orig- 
inally expected), in which case it would be indis- 
tinguishable from that initially added. Hence, the 


oxide ratio 0.57, calculated from the change in con- 
centration of the solution, not even approxi- 
mately represent the composition of any one solid 
phase. 
Accepting the 
the question of 


does 


there remains 
water of hydration. 


oxide ratio as 1 to 


the amount of 


; 


mately 1°C 


Data pe rlaining to preparation of Cal calerum aluminate solutions al approx 
Init concentration Fir concen i Mola 
Prepa olution Duratior ut l rat ! 
ration of storage Solid phases present at precipita API I I 
CaO ALO 
Ca) Ala) ( ‘) Ab 
lite Da 
\‘ 10-AbO, OHO | 
] 5 1. 278 44 plate probat 0. 108 0. 148 104 a 
oe mis hydrated Al»O 
. . . , f{CaO-AlyOs 1OHLO;: trace | f 
- S69 1. 1) 4 | 2CaO0-ALO. RHO { AA 14 I ) Fir : 
™ Gi sits CaO ( iL hey 
; 327 464 210 /Gibbsi i \ by LOHLO iv | my) ' ‘ f \ A ate ta 
| onal plate nh minor amount - ae . 
*» Gibbsite added (1 @ in 1 liter Calculated from change [th lut 
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FicurRE l. Electron micrograph of monocalcitum aluminate 
h ydrate pre paration 1). 
Magnification 32,500 diameters 


Assarsson |1] concluded from his analyses that the 
hvdrate contained 10 molecules of water. Longuet 
[2] reported approximately 7 molecules of water in 
his preparations after drying at 45-percent relative 
humidity. 

In the present study, preparation 2, after standing 
in a desiccator over CaCl, for 3 months, had an 
ignition loss equivalent to 7.13 moles of H,O per 
mole of AlO,. When placed over a saturated solu- 
tion of CaSO, at room temperature (relative humid- 
itv assumed to be 98%) it rapidly took up water 
until the ratio of H.O to Al.O, reached 9.9 at the 
end of 2 days, remaining essentially constant there- 
after. A similar preparation confined over saturated 
NH,CI (relative humidity about 79%) attained a 
ratio H,O/Al,O,=9.8. The results of these experi- 
ments strongly suggest that the fully hydrated 
compound has the composition CaQ-Al,O,-10H,O. 
(This formula, though not definitely established, is 
used in this report to designate the compound.) 

Three of the molecules of water apparently are 
very loosely bound, as they can be removed by 
drving over CaCl. 

Monocalcium aluminate hydrate was also prepared 
by mixing anhydrous monocalcium aluminate with 
a little water and storing the paste in the cold 
cabinet for 10 days. The X-ray diffraction pattern 
indicated the presence of CaQ-Al,O,-10H,O, with 
smaller amounts of 2CaQ-ALO,-SH,O and unhydrated 
CaQ-AlLO,. Under similar conditions, a white alu 
minous cement, consisting essentially of ALO, and 
CaQ, likewise, was partially hydrated to mono- 
calcium aluminate hydrate 

It has shown that monocalcium aluminate 
hydrate may be prepared at approximately 1° C 


been 








100m 


Photomicrograph of monocalcium aluminate hy- 


drate (preparation 3). 


FIGURE 2. 


Magnification 300 diameters 


There is evidence that it may also form at slightly 


higher temperatures. Farran’s preparations were 
made at 15° C, whereas Midgley (as reported by 
Lea [4]) obtained CaO-Al,O,-10H,O at 1° and 20° C, 
but not at 25° C. In the present study, a few 
experiments with pastes led to the formation of 
monocalcium aluminate hydrate at 25° C. Thus, 
when monocalcium aluminate paste was hydrated 
at 25° C for 10 days, the X-ray diffraction pattern 
indicated the presence of monocalcium aluminate 
hydrate, together with a somewhat greater amount 
of dicalcitum aluminate hydrate. As reported in a 
previous paper [3], similar experiments with the 
pastes of three aluminous cements at 24°C also 
gave evidence of the formation of monocalcium 
aluminate hydrate. These experiments have since 
been extended to include a white aluminous cement, 
known to consist largely of monocalcium aluminate 
and corundum. On curing for 7 days at 25° C, the 
cement was partly hydrated to CaQ-Al,O,-10H,O 
and a smaller amount of 2CaQ-Al,O,-8H,O. 

From the foregoing discussion, it is evident that 
monocalcium aluminate hydrate can be formed at 


temperatures as high as 25° C. In view of this fact, 


it is rather surprising that Wells, Clarke, and 
MeMurdie [6], in their study of the system 
CaQ— AL,O,—H,0 at 21° C, failed to observe any 


such compound. Perhaps the explanation may be 
found in the difference in experimental conditions. 
The experiments of Wells et al. were performed not 
with pastes but with relatively large volumes of 
solutions. 

There is evidence that in some earlier work of 
Wells [7], monocalecium aluminate hydrate was 
obtained at least once by precipitation from solution 
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at room temperature. The actual temperature, 
unfortunately, was not recorded. The solution was 
an aqueous extract of an aluminous cement. The 
crystals obtained were well formed and similar in 
appearance to those shown in figure 2. Another 
precipitated phase (presumably hydrated alumina) 


was also present; hence the composition of the 
crystals could not be determined. The following 
optical properties were recorded: 1.477; e=1.492; 
uniaxial positive. Apparently the phase was not 


observed again, although the study [3] included a 
large number of similar experiments. It is reason- 
able to assume, however, that the unknown crystals 
found by Wells were CaQO-Al,O,-10H,O. 

The completely hydrated crystals rr d in the 
present study had a refractive index of 1.471, with 
birefringence too faint to measure. After she crystals 
had been dried to a water content equivalent to 
CaO.-Al,O,-7H,0, the birefringence was stronger, and 
the minimum and maximum indices were 1.477 and 
1.480. The birefringence is still much lower than 
that reported by Wells for the crystals described in 
the preceding paragraph. No explanation for the 
discrepancy can be offered, because the degree of 
drying of Wells’ crystals is unknown. Farran [5] re- 
ported the indices, 1.489 +O. 002, « 1.507 
t0.002, but here again the amount of preliminary 
drying was not stated. 

Although the optical properties of the monocalcium 
aluminate hydrate were measurably affected by the 
degree of drying, no significant change was observed 
in the X-ray diffraction pattern. The pattern given 
in table 2 is for the material rehvdrated to the 
decahydrate stage. The pattern for the material 
dried to the heptahydrate stage showed lines agreeing 
closely with those in the table, with no apparent 
trend in either direction. The pattern in table 2 is in 
excellent agreement with that obtained by Midgley, 
as repor ted by Lea [4]. The difference between corre- 
sponding d-spacings does not exceed 0.02 A, except 
for the first line, for which the measurement error is 
necessarily high. Also, a few lines are reported 
doublets in one pattern and as singlets in the other; 
and a very few weak lines are present in one pattern 
and not in the other. Good general correlation is also 
observed between these patterns and those of Assars- 


son [1], Longuet [2], and Farran [5], if allowance is 
made for the lines due to other phases present. 


Nevertheless, a few discrepancies are worthy of men- 
tion. The pattern given by Assarsson fails to show a 
spacing near 14 A, although the patterns reported by 
the other investigators show a strong line in. this 
region. Limitations of the earlier X-ray equipment 
may account for this difference. The largest spacin 

reported by Assarsson was 7.68 A, whereas the le 
closest to this in the other patterns ranges in position 
from 7.16 to 7.3 A. No explanation for this discrep- 
ancy can be advanced at present. In other respects, 
the pattern of Assarsson agrees well with that given 
in table 2. Farran’s pattern is likewise in agreement, 
except that a few of the weaker lines are absent. 
-Longuet’s pattern shows the same series of lines, but 
for some reason they are all displaced toward higher 
d-spacings in comparison with the other patterns 


TABLE 2 V-ray powde ? diffr action patterns of 
CaO-Al,Oy-lOHLO and SrO-ALO,-nH.O* 
CaO0-AbO.1lOHLO SrO-ADO HO Ca0-AhkO»-lWHLO SrO-AbOy-n Hor 
d I l | 1 I d I 

i 1 i 
14.2 “) 14 1) 2. 18 7 2. 22 29 
lf 100 7.27 wy 2.17 7 
7 s 5.4 4 2.11 7 2.14 » 
4.73 10 1.8 Ds 2. 07 5 2. 11 23 
4. 51 i +. 61 8 2 03 l 
1.16 ‘ 4.23 6 1. 969 2 2. 02 23 
+ OH 2 $ (4 5 1. 04 6 1.97 21 
4.70 7 3. 74 15 1. 927 l 
$. 56 13 4. 65 MH) 1.879 2 1. 91 10 
$. 26 7 1.3 $3 1, Sat 1, 8S l 
3.10 6 3.14 15 1. 786 $ 1.84 s 
3. OS - 1. 757 2 
2. 91 7 2. OS 45 1.713 2 1.74 10 
2. Sh s > yl 31 1. 60S ] 
2. 70 7 20 1, 650 $ 1.68 7 
2. 11 2 2 1. 637 l 
2.48 6 2. 53 1. 630 l 
2 46 ‘ 2.48 s 1. 604 1 1. (4 s 
2. 37 s 2.43 10 1. 555 l 
2.34 6 2. 38 14 1. 524 l 
2. 30 l $3 1.472 
2. 27 11 2. 30 2 1, 388 1.42 t 
* Obtained on a Geiger-counter diffractometer, using copy Ka radiation 
8.4 in sample used for determination of X-ray pattert i may rang 
7tolW 
interplanar spacit [=relative intensity 


3.2. Monostrontium Aluminate Hydrate 


In an earlier paper [8] on the strontium aluminates, 
it was reported that only the cubic tristrontium 
aluminate hydrate was formed by the various meth- 
ods emploved in that study. However, none of the 
experiments had been conducted below room tem- 
perature. It was, therefore, considered of interest 
to carry out some additional reactions in the cold 
cabinet. 

In one experiment, 
aluminate was prepared by 
strontium carbonate and alumina at 1,460° C, fol- 
lowed by fusion in an oxygen blast. A small amount 
of the material was ground and placed in a vial with 
just enough water to moisten it. It was then stored 
in the cold cabinet for 10 days. An X-ray pattern 
of the hardened mass showed a series of lines very 
the pattern of monocalcium aluminate 
Lines of unhyvdrated SrO-Al,O; were also 


anhydrous monostrontium 
heating a mixture of 


similar to 
hvdrate. 
present. 

In another experiment, the new hydrate was pre- 
pared by a precipitation method analogous to that 
used for the calcium aluminate hydrate. ' For this 
purpose, cir g wig monostrontium aluminate was 
first prepared by heating a mixture of the hydrated 
oxides at 1.400° C. (¢ ‘ombination was not complete, 
but further heating was avoided because increased 
crvstal size would tend to lower the rate of solution. 
The product was ground, shaken with water, and 
filtered. and the filtrate was stored in a glass flask in 
the cold eabinet. The data relevant to two such 
experiments are given in table 3. 

Because of the extreme fineness of the precipitate, 


the optical properties could not be accurately 
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TABLE 3. Data pertaining to preparation of SrO-AbO,-10H,O 
from strontium aluminate solutions al approximately 1° C 


Initial concen Dura- Molar 
tration of tion Solid phase ratio 
Prepa solution ol present at end SsrO Appearance of 
ration stor of storage AleO;in precipitate 
ipt prepa 
sro ALO ration 
g/liter g/liter Days 
| 443 2 62 14 SrO-AbO+ WHO 100 Needles 
> 6. 92 5 6 2 SrO-.AbDO»1l0HLO 0. OS Slender prisms 
- about + 


long 


The crystals were acicular and slightly 
with a mean index of refraction about 


determined. 
birefringent, 
1.478. 

The water content was determined by methods 
similar to those used for the calcium compound 
After exposure to high humidity (over saturated 
Ca(OH). solution) at 1° C, the H,O/AI.O, ratio was 
99. whereas in material dried over CaCl, the ratio 
was 6.9. Thus, the strontium and the calcium 
compounds appear to be analogous in their response 
to humidity changes. However, it was found neces- 
sary to keep the strontium aluminate hydrate in the 
cold cabinet while exposing it to the higher humidity, 
because it decomposed at room temperature. The 
products of decomposition, identified by microscopic 
and X-ray methods, were 3Sr0-Al,0,-6H.O and 
gibbsite (Al,O,-3H.O). 

The X-ray pattern given for the strontium alumi- 
nate hvdrate in table 2 was obtained on a portion of 
preparation 1, which had been dried to a water con- 
tent of about 8.4 moles per mole of SrO-Al,QOs. 
Because of the acicular habit of the crystals, the 
relative intensities indicated may be somewhat 
affected by orientation. The similarity between the 
patterns of the strontium and calcium aluminate 
hydrates is readily apparent. 

The question may be asked whether barium might 
also form a monoaluminate hydrate analogous to 
those of calcium and strontium. In an earlier paper 
9} a barium aluminate hvdrate, BaO-Al,O,-7H.O, 
was described. The amount of water of hydration 
is the same as that of the monocalcium and mono- 
strontium aluminate hydrates dried CaCl,. 
However, the barium compound differs from those of 
calcium and strontium in ervstal habit, indices of 
refraction, and birefringence, and gives a distinctly 
different X-ray pattern. It seems clear, therefore, 
that BaO-Al,O,-7H,O is not an analog of the mono- 
calcium and monostrontium aluminate hydrates. 


over 


TOO 


4. Summary 


Monocalcium aluminate hydrate, probably having 
the composition CaQ-Al,O;-10H,O, was prepared by 
precipitation from calcium aluminate solutions at 
approximately 1° C, and by hydration of pastes of 
monocalcium aluminate and of aluminous cement 
at 1° and 25° C. Under conditions favorable for 
slow crystal growth, this hydrate crystallized as 
hexagonal prisms with terminal pyramids. They 
were very weakly birefringent, with a mean refractive 
index of 1.471. Drying over calcium chloride re- 
duced the water of hydration to 7H,O, and the indices 
and birefringence increased, but the X-ray pattern 
was apparently unchanged. 

An analogous monostrontium aluminate hydrate 
was prepared by similar methods, but it was obtained 
only as minute needles or prisms. Like the calcium 
analog, it probably crystallizes with 10 molecules of 


H.O, 3 of which may be removed by drying over 
calcium chloride. The mean index of refraction is 
about 1.478. When placed in a humid atmosphere 


at room temperature it decomposed, with the forma- 
tion of 3Sr0-Al,0;-6H,O and Al,O,-3H,O (gibbsite). 

X-ray powder diffraction patterns are given for 
both monocalcium aluminate hydrate and mono- 
strontium aluminate hydrate. The similarity be- 
tween the two is readily apparent. 


The author acknowledges with thanks the assist- 
ance of several members of the Constitution and 
Microstructure Section of the Bureau, who prepared 
the X-ray diffraction patterns and the electron and 
photomicrographs. 
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Turbidity and Viscosity Measurements on Some Cationic 
| Detergents in Water and in Sodium Chloride Solutions 


Lawrence M. Kushner, Willard D. Hubbard, and Rebecca A. Parker 


Light-scattering and viscosity measurements were made on solutions of three cationic 


detergents in distilled water and in various solutions of sodium chloride. 
hydrochloride, dodecyltrimethvlammonium chloride, and tetradecyl- 
The first of these was investigated at 30° C 
The micellar weight and the intrinsic viscosity of the micelles of each detergent 


were dodecvlamine 
trimethvlammonium chloride. 
; at Ze CG 


The detergents 


and the others 


were calculated, and the dependence of these quantities on the concentration of sodium 


chloride is discussed 


1. Introduction 


The anionic detergent sodium dodecyl! sulfate has 
been the subject of a number of physical chemical 
investigations. Studies have been made on dilute 
aqueous solutions of the detergent with respect to 
light scattering [1,2],' viscosity [3], and electro- 
phoresis [4]. As a result, the micellar behavior of 
the material has been well characterized. 

A particularly interesting feature of these investi- 
gations is the marked dependence of the physical 
property being measured (turbidity, viscosity, etc.) 
on the quantity of simple electrolyte in the solutions. 
In the case of light-scattering measurements, it is 
found that the addition of sodium chloride up to a 
concentration of about 0.05 to 0.06 M causes a large 
increase in the turbidity of the solutions. Accord- 
ingly, the micellar weight calculated from the data 
by the method of Debye [5| or Mysels [6] increases 
rapidly with increasing salt concentration in this 
range. The relative viscosity of the solutions de- 
creases markedly in the presence of similar concentra- 
tions of simple electrolyte. This is reflected in the 
intrinsic viscosity which decreases rapidly in_pre- 
cisely the same range. 

This paper describes the results of turbidity and 
viscosity Measurements made on solutions of three 
cationic detergents in distilled water and in solutions 
of sodium chloride. Dodecylamine hydrochloride 
was chosen because its surface active ion is the posi- 
tively charged analog of the dodecyl sulfate ion. 
Because of the poor solubility of the primary amine 
hydrochloride in salt solutions, measurements were 
also made on a comparable quaternary salt, dodecyl- 
trimethylammonium chloride. Tetradecyltrimethyl- 
ammonium chloride was investigated to give an 
indication of the effect of increasing the chain length 
of the hydrophobic portion of the surface active ion. 


Se 


' Figures in brackets indicate the literature references at the end of this paper, 
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2. Experimental Procedure 
2.1. Materials 


Dodecylamine hydrochloride was prepared from 
dodecylamine in the following manner. Fifty grams 
of dodecylamine (Armour & Co.) was dissolved in 
anhydrous benzene. Gaseous hydrogen chloride 
was bubbled through the solution until a slight excess, 
as indicated by a methyl redspot test, had been 
added. The solution was cooled to approximately 
0° C, and the dodecylamine hydrochloride crystals 
were filtered off at that temperature. The material 
was then recrystallized twice in the same manner 
from benzene and finally dried overnight in a vacuum 
oven at 60°C. As a check on the purity of the salt, 
the specific conductance ef a 0.01354-N aqueous 
solution was determined at 30° C. A value of 
13.62 10-* ohm! was obtained. This to be 
compared with the value 13.64107-* ohm” given 
by Armour & Co. The chloride content was deter- 
mined to be 15.96 percent. The theoretical value 
is 15.99 percent. 

Subsequent measurements of the turbidity of 
aqueous solutions of the dodecylamine hydrochloride 
indicated the presence of traces of a water-insoluble 
impurity, suspected to be benzene. ‘Two recrystal- 
lizations of the salt from cold methyl alcohol elimi- 
nated the contaminant. 

Dodecyltrimethylammonium chloride and _ tetra- 
decyltrimethylammonium chloride were obtained as 
the quaternary salts from Armour & Co. Both were 
purified by precipitating them from solutions of 
ethanol by the addition of anhydrous ethyl ether. 
The resulting slurry was cooled as before, and the 
crystals were removed by filtration. They were 
then dried under vacuum at room temperature. 


Is 


2.2. Refractive Index 


Measurements of the refractive-index increments 
(An/Ac) of the three detergents in distilled water and 
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TABLE |. Light-scatte ringand viscosity data for solutions of dodecylamine h ydrochloride, dodecyltrimethylammonium chloride , and 
tetradecyltrimethylammonium chloride in distilled water and in solutions of sodium chloride | 
ep=centipoise; [y]=intrinsic viscosity; An/Ac=refractive index increment; C MC =critical micelle concentration; [Z]=intrinsic dissymmetry 
; 
DODECYLAMINE HYDROCHLORIDI wr 
Solvent Concentra Density Viscosity ” Ceneentra Turbi tity An; XM Micellar 
tion tion weight j 
yidl yiml cp a/mlX10 cm <10 nlig 
0 0. GUS565 0. 7975 0 4 76 
0. 100 GO556 | 0. O89 500 
200) QU548 | 222 4. 95 | 
300) Y953Y SOM) | 25 5. 33 
| Kn) Q94531 SOUS 343 6. 45 
50M) 99523 si48 | 363 7.81 | 
Hoo WO515 SISSY | un) 12. 23 
iM) YO4US 8327 439 15. 28 | 
Distilled H.O | 1. 000 YWU482 S40 » ond OS *). 42 oO 161 1 @ 
CMC =0.323 g/dl 1. 100 9473 S519 Hl 24.18 : . 
1. 200 W465 Su? Hao) 31. 42 
sw) $8. 70 
| | “OT 1). ON 
| W27 4). 89 
1. 027 43.03 
1. ISS 45.17 
1. 428 iS) 
Has 1. 31 i 
( 4+) \ 
0.0 4 
12] 5. 09 
1) 17 
21s 3s 
| 275 U5 | 
0.01—M NaCl 320 12.00 \ 16 >» mM 
CMC =0.272 g/d 10 21.13 —_ 
510 24. 04 | 
| “Hlo 44.84 
735 40. 41 | ' 
| S15 43. 32 ' 
| | 1. 46) 
1. 21 4. 36 
t 0. Wun 0. 74u90 “ $76 
0. 100 GUB3S 0. 123 5. 24 
aM) 10630) S037 2h 12. 42 
$M) aue21 SO7T5 300 17.37 
1M) uuBnLS SIIS in) 30.5 } 
0.02 VW NaCl wi) WUHOS S16] 0) Hi). 1 
CMC =0.227 g/d wn 44546 8205 0. O4S foo is. | Oo 10 4 100 
| GOSS) S2US TO) 55.2 j 
1 (HM YWOSHS S307 wi) H1.4 
wow) ih 
ooo 71.9 
1. 100) 76.6 
1. 25) 83.3 
a 0. WUT2QA 0. 8006 0 5. 24 
1H) W717 SUAS 0.103 5. 38 
a) v9709 SOS! IST 1. 57 | 
sim) “VyT02 SOY aH) 12. 44 
Hn) GU8O4 S114 300) 30.5 
004—M NaCl uM) YUBS6 S150 un) 45.2 
CMC =0.165 ed | How QOBTS slau mL a 
iM) GUHHS SINT Ogu TL 70.8 ( Ta Is 0) 
ooo uunay sl 7) S16 
si) 91.7 
GOO) w1.2 
1. (4M) 1.4 
| 1. 100 1.0 
1 20 Ww 
’ 2) 133.8 
HUST 0. S02] it 5. 24 
i Win WU7TUS SOU] | oO] 5. 24 
AM) YUTSU mwoyl an) AS | 
40M 780 sli si) 12s | 
Hx QuTTl 8135 HM) 62.8 
0.06 VJ Nat TD w7H2 SI1S6 | vn) si) 4 | 
CNC =0.145 g/d Hw WUT? S28) | Oro 105.7 
. a $4 = 0. 1 $4, (m0 
) sii VuT34 stun 4) 271.9 . 
uy va7 te Suu | sin) 129. 5 
con 145 
1 oO) 1S 
Lon 16s 
| 1. 100) 171 
1 





in the sodium chloride solutions were made with a 
differential refractometer similar to that described 
by P. P. Debye [7] 
at 23° C 


index increments determined are given in table | 


All measurements were made 


using blue light (436 mu The refractive- 


2.3. Turbidity 


Measurements of turbidity were made in semi- 
octagonal cells in an absolute light-scattering photom- 
eter |S]. Semioctagonal cells were used so that 
determinations of dissymmetry, 7 (the ratio of the 
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TABLE |. Light-scattering and viscosity data for solutions of dodecylamine hydrochloride, dodecyltrimethylammonium chloride, and 





























( 
ind tetradecyltrimethylammonium chloride in distilled water and in solutions of sodium chloride -Continued 
cp=centipoise; [n]=intrinsic viscosity; 4n/ic= refractive index increment; C MC =critical micelle concentration; [Z]=intrinsic dissymmetry) 
; 
Solvent Concentra- Density Viscosity n Concentra- Turbidity An/Ac Micellar 
tion tion weight 
ht | g/dl g/ml cp g/ml X 102 em 105 ml/g 
0 0. 99890 0. 8036 \ 0 5. 24 
0. 100 YYRS | 8057 0). 037 5.14 
200) YOSTS SOS 097 5. 28 
300 YUsti4 8124 155 17. 47 
0 YyYR55 S165 200 32.58 
0.08-M NaCl 500 YUs46 8207 300 67.8 . . 
; CMC =0.118 g/dl 600 GOR3S 8254 0). O38 340) 814 0. 160 51, 100 
S00 YES R355 461 122 
i 1. 000 YesOS 847Y 612 167 
794 222 
926 261 
} 1.000 283 
iM 
0 0. 99972 0. 8049 0 5. 42 
0. 100 G4 SOY 0. O99 5.47 
200 gO955 S004 1) 31.9 
300 W946, S148 315 27.6 
- HM Y943S S202 385 167 
; 0.10: VeNaCl FM) 99929 N20) , 0.044 { 440) 225 0. 160 85, 500 
CNM « 0.109 ——— poee > 
600 49921 8327 ESS 284 | 
aO0 4yudgo4 S484 755 3S1 
j 1. 000 YOSST SoS? SAL) 454 
¥70 506 
1. 000 524 
0 1. OO178 0. 8090 f 0 5. 38 
0. 100 1. OO169 S108 0. 050 ». 47 
200) 1. 00159 S167 O74 5. 52 
300 1. 00149 S261 107 11. 57 
ww 4100 1. 00139 S3N4 130 47.6 
0.15—-M NaCl soo 1. 00130 S528 205 162.58 9 
CNC =0.101 gal 600 1. 00120 S604 0, OU 330 353 0. 160 231, OO 
' S00) 1. OO101 v101 is) 55S 
1. 000 1. OOO81 VOSS 592 700 
780 919 
| 916 1076 
1. 000 1176 } 
0 1. 00367 0. 8129 0 5. 52 
0. 100 1. 00359 S167 0. O50 5. 42 
200 1. 00350 8350 107 61.9 
} 300 1. 00342 S610 171 346 
0.20--M NaCl / 410) 1. 00333 SUH 0.15 281 #43 0. 160 655, 000 
mM) CMC =0.098 g/dl 500 1. 00324 “401 40S 1257 
' 600 1. 00316 yuoy 67 1950 
S00 1. OO2U8 1. 1204 769 2220 
| 1. 000 1. OO2S1 1. 2065 | 901 2550 
1.000 27M 
{ 0 1. OOSS0 0. S156 
| 0. 100 1. 00570 S270 | 
, 200 1. OO55Y SOTO 
P aie a a 300 1. 00548 9230 0. 33 7\™1. 20 ~1,125,000 
: . 1K) 1. 00537 1. 0060 
| How 1. OO516 1. 2701 | 
SOO 1. 00494 1. 7631 
0 1. OO785 0. 8201 
we) 0. 100 1. 00772 S3s4 ; 
0.30-M Nat 200) 1. OO75S V102 0. 43 7\™1. 40 ~ 1,800,000 
C\C =0.085 g 300 1. 00745 1. 0446 : 
1) 1. 00731 1. 2009 
500 1.00717 1. 6907 
DODECYLTRIMETHYLAMMONIUM CHLORIDE (23° ¢ 
0 0). O74 0. 9323 0. 183 5.47 
| 0.404 W703 9491 225 5. 62 
| As GUTS QAI5 27 >». 06 
WOT GO6GS u546 478 ». 71 
TOS QU687 GPRS 545 5 71 
L) st Ik 1 H20 SOO) UU67S UiSl 0. 040 20 &. 81 0. 156 9, 900 
- i CMC =#0.50 1.013 QO655 a744 709 13. 09 
1. 202 Y9636 {S853 7¥3 16. 66 
1. 402 GQU615 uu74 O55 21.75 
i 1. Ov2 24.09 
} 1. 197 26, 89 








scattered intensity at 45° to that at 135°), could be | ing concentration were prepared by the addition of 

made as a check on the cleanliness of the solutions. | a known weight of detergent crystals to clean solvent 
.emi- The solvents, distilled water or solutions of sodium | or solution already in the light-scattering cell. After 
tom- chloride, were easily obtained with no dissymmetry the detergent crystals had dissolved, the solution 
that by filtering two or three times through ultrafine | was filtered a number of times until the dissymmetry 
Pthe ¢  fritted glass or porcelain filters. Solutions of increas- | was at an acceptably low level. Except for dodecyl- 


115 








TABLE 1. 
tetradecyltrimethylammonium chloride 


(ep=centipoise; [y]= intrinsic viscosity; An/Ac=refractive 


Solvent Concentra Density 
tion 
g/dl 
0 a 
0. 360 
nw) 
0.02-M NaCl OOO) 
CMC =0,.420 g/d SOU) 
1. O00 
1. 200 
1. 400 YOBYS 
0 0. 99914 
2m) YYUSS2 
300 YYRSO 
0.04-M NaC] po pons 
CMC =0.310 g/dl YR55 
sin) YYRSH 
1. 000 YOR 17 
1. 20 Q9799 
1. 400 G9TSO 
{ 0) 1.00149 
0.150 1.00134 
200) 1. 00130 
1) 1. OO113 
0.10-M NaCl 600 1. 00006 
CMC =0.190 g/dl \ SOK) 1. 00079 
1. 000 1. OOOB82 
1. 200 1. 00045 





TETRADECYI 








0 0. GU7M4 
0. 100 vu743 
| a) W732 
Distilled H,O } 300 9721 
CMC 0.120 e/d wn “9710 
HOw GUBSS 
si) YYhbbD 
1. O00 W444 
( 0. YUR2Z5 
0. 100 YURLS 
An) UYSOIH 
300 “y7e7 
0.02-M NaC un) GUTS7T 
CMC =0.070 e/d \ Hoo QuT7HY 
SOO) WoT 
1.000 W731 
0 0. vu914 
0. 100 YOGA 
20) GUSG4 
300) GUSS4 
0.04-M NaCl TTD YOR7T4 
CMC =0.040 g/dl OO YUR54 
SO) YUSS4 
1. 000 YYR 13 
0 1.00149 
0. 100 1. 00142 
200 1. 00133 
! 300 1.00124 
0.10-M NaCl no) 1. OO115 
CMC =0.030 g/dl 600 1. OOOO 
S00 1. 00080 
1. O00 1. OOOB2 
amine hydrochloride in the presence of high con- 


centrations of sodium chloride, a small dissymmetry, 
generally 1.05 or less, was observed. In the case of 
dodecylamine hydrochloride, the dissymmetries ob- 
served are caused by the presence of very large 
micellar structures and are of significance. This is 
discussed further in a later section. 


116 


index increment; CMC 





tIMETHYLAMM¢ 


V iscosity 


9349 
9452 
U4 
9491 
9571 
9653 
9744 
9R31 
¥923 


9349 
¥459 
O476 
U540 
v621 
WOOT 
O779 
YS5Y 


INIUM CHLORIDE 


S23 
W355 
YSY6 
0466 
W554 
9732 
WSs 
0121 


¥S2y 
367 
9406 
9451 
OAD 
YODS 
9712 


YS20 


W349 
9379 
9422 
464 
502 
9592 
YOST 
O7S86 


YS3Y9 
9438 
G44 
Y50Y9 
9546 
9616 
Y6YS 
9779 
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0. 040 


0. 035 


0.042 


0. 040 


0. 040 








the 
hydrochloride in salt solutions at 23° 
t 30° 


Concentra- 


tion 


measurements on it were made 


turbidity 


out. 


Turbidity 


a/mlX102 cm! X10 
f 0 4.581 
0. 390 5. 52 
431 6. 24 
178 9.14 
583 15.4 
737 22. 8 
935 30.9 
1. 139 37.8 
1. 249 41.4 
0 1.8] 
0.179 5. 43 
323 6. 52 
373 9.71 
478 17.1 
638 27.0 
765 34.3 
SOO) 41.2 
1. 084 OY 
0 ». O5 
0. 144 5. 24 
24 6.71 
51 11.4 
342 20.3 
445 29. 5 
5) 38. 2 
738 MS 
¥2h 2.9 
1. 048 71.1 
0 1S! 
0.044 im 
143 6.47 
104 9.43 
334 14.9 
431 17.8 
513 19.2 
616 2.8 
0 4.1 
0. 054 5.14 
105 Ww. 1 
1tit) 17.0 
a) 22.7 
2 | OS 
344 32. 
407 37.2 
{Us 2.7 
WS 45.2 
0 4.72 
0.043 5. 52 
ous 1? 
154 20.3 
229 2v 
308 38.3 
41 47.1 
493 55. 0 
S85 63.3 
680) 69.3 
0 1 uw 
0. 054 9. 76 
112 18.8 
178 28.7 
255 39.4 
329 5.6 
iy 61.3 
$4) 72.1 
io2 84.7 
692 93. 4 
low solubility 
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C, the turbidity 
C. All other 
The 


436-my line of a mercury-arc lamp was used through- 


measurements were fash 
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=critical micelle concentration; [Z]= intrinsic dissymmetry 


An; AX Micellar 
weight 
ml/g 
0. 156 10, 200 
0.156 11, 000 
0. 155 16, 400 
0. 157 10. 100 
0. 157 21, 800 
0. 157 22, 200 
0.157 25, SOO 
dodecylamine 
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Light-scattering and viscosity data for solutions of dodecylamine hydrochloride, dodecyltrimethylammonium chloride, and 
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2.4. Viscosity 


r The experimental details and necessary corrections 
to be applied to the measured efflux time have been 
described previously [3]. The only change in pro- 
cedure was the use of smaller pycnometers for the 
density determinations [9]. The measurements on 
the dodecylamine hydrochloride solutions were made 
at 30° C. All others were made at 23° C. 


3. Results and Discussion 


All of the pertinent experimental data are given 
in table 1. Typical plots of turbidity and viscosity 
as a function of detergent concentration at various 
concentrations of sodium chloride are shown in 
figures 1 and 2. The pronounced effect of the ad- 
dition of sodium chloride is evident in both figures. 

In order to evaluate the micellar weight from the 
turbidity data, and the intrinsic viscosity of the 
micelles from the viscosity data, it is necessary to 
evaluate the critical micelle concentration of the 
detergent in each of the solvents used. Debye [5] 
has shown that this can be obtained conveniently 
from the dependence of the turbidity on concentra- 
tion (ef. fig. 1). The detergent concentration at 
which the turbidity begins to increase rapidly 
clearly represents the concentration at which aggre- 
gation of the individual detergent molecules or ions 
into micelles is beginning to occur to a significant 
extent. These concentrations, as obtained from the 
turbidity data, have been chosen the critical 
micelle concentration, and are given in table 1. 
The dependence of the critical micelle concentration 
of each detergent on the amount of sodium chloride 
present is shown in figure 3. The data reported by 
Corrin and Harkins (dye titration method) [10] for 
dodecylamine hydrochloride at 26+2° C are also 


as 


shown. It is seen that the quaternary dodecyl- 
ammonium chloride has a much larger critical 


micelle concentration than does the corresponding 
primary amine salt. This is simply another mani- 
festation of the greater degree of ionization of the 
quaternary salt. The addition of two CH, groups 
to the long hydrocarbon chain of the surface active 
ion, as in tetradecyltrimethvylammonium chloride, 
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Dependence of the turbidity on the concentration of 
chloride in water (()) and in 
various sodium chloride solutions 
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CIGURE | 
dodecyltrimethylammonium 


| 
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increases the hydrophobic nature of the ion and 
results in a lower critical micelle concentration. 

The intrinsic viscosity of the micelles was obtained 
from plots of the experimental data according to 
eq (1). 

(rei 1) 


( =[n]+ACu. 


(1) 


m 
[n] is the intrinsic viscosity. Cy, is the concentration 


of micelles in grams per deciliter and is assumed to 
be equal to the total concentration of detergent 
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Figure 2. Dependence of the viscosity on the concentration of 
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letradecyltrimethylammonium chloride in water (QC) and in 
various sodium chloride solutions. 
@, 0.02 M; @, 0.04 M; ©, 0.10 M. 
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FiGuRE 3. Dependence of the critical micelle concentration of 


dodecyltrimethylammonium chloride (@©), dodecylamine hydro- 
chloride (CQ), and tetradecyltrimethylammonium chloride (@) 
on the concentration of sodium chloride in aqueous solutions. 


@, data of Corrin and Harkins [10] for dodecylamine hydrochloride at 26° C, 
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minus the critical micelle concentration. Nrey IS the 250 r-— 10.5¢ fel 
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viscosity of a solution of concentration C,, relative a SS: 
to the viscosity at the critical micelle concentration ’ | 9 . 
(i. e., at C,=0). The constant A is related to 200} lose \ a lo 
hydrodynamic interactions in the solution. The we (9) | oo 
relationship between [n] and a shape factor, A, for | 7%, ne i 
suspended particles of density d, is B al 1S pose 
0.30 mol 
- : for 
2 2 | | met 
LOOd 4 100} 1o.2¢ the 
¥ , } carr 
For spherical particles, K=2.5. Values of K as a | * we 
function of axial ratio for oblate and prolate sphe- ” 040 oT 
roids of revolution have been calculated by Mehl, ? —_ « “ 
Oncley, and Simha [11]. ore —e—e—* ct 
Calculation of the micellar weights from the light- O% v aes 4 wit 
scattering data was done with the simple Debve MOLARITY OF Ne : - 
equation , 
( 1 Figure 4. Dependence of the micellar weight (©)) and den 
Be ax +2BC,,. 3) intrinsic viscosily (@) of dodecylamine hydrochloride micelles 1) 
T M on the concentration of sodium chloride in aqueous solutions diff 
where ma’ 
H lee 4) a —_ —_ _y as 
» 4 0.1 
Once again C, is the concentration of micelles, | » - 
although this time in grams per milliliter. 7 is the | 2 | the 
micellar weight. £ is an interaction constant; its | 2 9), ——— 15 
significance in solutions of charged colloidal particles ” | da 
has been discussed [12]. 7 is the turbidity due to os a rn 
the presence of micelles. It is the observed tur- — 1 spl 
bidity of a solution minus the turbidity at the critical | of 
micelle concentration. The refractive index of the | .°* - ¥ —4 res 
solvent is No. (An/Ac) is the refractive-index increment Ss bri 
of the solutions. A is the vacuum wavelength of oser 7 me 
the light used, 436 mu. N is Avogadro’s number. | 
In the case of dodecylamine hydrochloride at —— “oon ——a 6 
sodium chloride concentrations above 0.20 A/, the MOLARITY OF Ne Ci or 
micelles are large enough to scatter asymmetrically. |, : on 
mM: a. ‘ . cha Rea Figure 5 Dependence of the micellar weight ) and . 
rhis means that the scattered intensity observed at a le ts sca A al ne: eR BS ol 
90° is not indicative of the turbidity of the solutions micelles on the concentration of sodium chloride in aqueow } ev 
and that the micellar weight calculated therefrom is solutions ms 
in error. Because no dissymmetry corrections were + om 
applied, turbidities for these systems are not given uy 
in table 1. An estimation of the intrinsic dissym- i eh A ne ee ee oe ee ae vi 
metry has been made, and these values are given ad 
The micellar weights listed are approximations “ m 
calculated from the data as if there were no dis- 2 20 m 
symmetry effects. s / ds 
Plots of the intrinsic viscosity and micellar weight f eX 
as a function of sodium chloride concentration for io4 | () 
each of the detergents investigated are given in t cl 
figures 4, 5, and 6. erie | th 
Considering both .V/ and [y]| over the total range of - rae | dl 
sodium chloride concentrations covered, each of the ~ 004} . oa 4 m 
three detergents investigated shows different be- 0.02! ” 
havior. However, considering only the dependence ea ae me 
of {[n] on sodium chloride concentration up to about 3 005 O10 bs 
0.05 \/, the three detergents behave similarly, i. e., MOLARITY OF NoCI 
the intrinsic viscosity decreases as sodium chloride ; 
. ' FIGURE 6 De pe ndence of the micellar weight | ) and intrinsic 


is added. his 1S the same behavior as that exhibited riscosily | e of tetradec yltrimeth ulammonium chloride mice lles 
by sodium dodecyl! sulfate, and reflects the diminution on the concentration of sodium chloride in aqueous solutions ‘ 
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ef electroviscous effects caused by the high degree 
of ionization of the micelles.” 

At a concentration of about 
chloride, the electroviscous effect 
virtually eliminated. An interesting result 
electroviscous effect for 
chloride is substantially 
for dodecylamine hydrochloride or tetradecyltri- 
methylammonium chloride. The inference is that 
the micelles of dodecyltrimethylammonium chloride 
carry a lower charge than the micelles of the other 
two detergents. However, no explanation for this 
can be given at this time. 

The minimum value of [yn], 0.033 to 0.040, obtained 
with all of the cationic detergents is to be compared 
with the value 0.035 obtained with sodium dodecyl 
sulfate. All of these values are consistent with the 
concept of spherical hydrocarbon micelles having a 
density of about 0.75 g/ml. According to eq 2, 
| for such micelles would be about 0.033. The 
differences among the intrinsic viscosities observed 
may reflect a difference in the density of the micelles 
in each case. 

At higher concentrations of sodium chloride (up to 
0.1 M) the intrinsic viscosities of the two quaternary 


0.05 Af sodium 
appears to be 
is that 
dodecyltrimethylam- 
than that 


the 
monium 


less 


salts remain constant in spite of the increasing 
micellar weight. This merely serves to emphasize 
that small micelles (those containing about 50 to 


regardless of their pre- 
cise shape, will necessity have axial ratios close 
to unity and will behave hydrodynamically as 
spheres. The growth of the micelles on the addition 
of neutral electrolyte has been interpreted [5} 
resulting from a reduction in the work necessary to 
bring together the ionic heads of the detergent 
molecules when a micelle is formed. 

With dodecylamine hydrochloride the addition of 
sodium chloride beyond 0.05 7 results in rapid 
growth of the micelles. This is made evident not 
only by the micellar weight but also by the intrinsic 
viscosity, which increases greatly. Once again it is 
evident that the light-scattering method is much 
more sensitive for measuring the growth of the 
micelles than is the viscosity method, particularly 
up to micellar weights of about 250,000. The high 
values of [n| that are obtained as .V/ increases are 
additional evidence for the high asymmetry of large 
micelles. The quantitative correlation between 
micelle dimensions calculated from the viscosity 
data and from the light-scattering data is poor. For 
example, at a sodium chloride concentration of 
0.3 M, |n]}=0.45, and the intrinsic dissvmmetry is 
close to 1.40. Taking a rodlike micelle as being 
the most likely [13], one calculates from the intrinsic 
dissymmetry [14] a length of about 1177 A. Approxi- 
mating the rod to an elongated ellipsoid of revolution 
whose minor axis is equal to the extended length of a 


150 detergent molecules), 


For a brief 
by H.R 


effect, see Colloid science, edited 


Co., New York, N. Y., 1952 


discussion of the electroviscou 
Kruyt, vol. I, p. 348 (Elsevier Publishing 


seVaETOO y 


-! 
= 


dodecylamine hydrochloride molecule, 20 A, one 
calculates [11] a length of 720 A from [ny]. It must be 
mentioned, however, that for this system both 


Z and (n-i:—1) 
of Cy, as C,->0. 
the intrinsic 


(\, are rapidly changing functions 

Hence the probable errors in both 
dissymmetry, [Z], and in the intrinsic 
viscosity, [n], are quite large (+10 to 15% would not 
be unreasonable). In addition, with such large 
asvmmetric micelles, one must consider the possi- 
bility of non-Newtonian flow. This was not investi- 
gated. 

The rapidly increasing size of the dodecylamine 
hvdrochloride micelles is probably the first step in 
the salting-out process. On standing at 23° C, 
visible crystals were observed in 0.5-percent solutions 
of the detergent in a 0.04-M sodium chloride solution. 
At 30° C it was possible to prepare 1-percent solutions 
of the detergent in sodium chloride solutions as 
concentrated as 0.30 MM, but at higher salt concentra- 
tions it was insoluble. The quaternary salts, being 
inherently more soluble, show no such effects in this 
range of sodium chloride concentrations. 


In the case of tetradecyltrimethylammonium 
chloride in distilled water, because it has a low 
critical micelle concentration, it was thought that 


it might be possible to observe dissvmmetries of less 
than unity, as discussed by Doty and Steiner [14]. 
Such effects should exist, being the result of external 
interference arising from the nonrandom fluctuations 
in concentration of the highly charged micelles. 
However, dissymmetries of less than unity were not 
observed. Apparently, even in this case, in which 
the shielding electrolyte is unaggregated detergent 
at its critical micelle concentration (0.0043 .), 
shielding is sufficiently good to permit random fluctu- 
ations, at least up to Cy™1 g/dl. 
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Relative Stress-Optical Coefficients of Some National Bureau of Standards 
Optical Glasses’ 


Roy M. Waxler and Albert Napolitano 


The relative stress-optical coefficients of 27 optical glasses made at the National Bureau 


of Standards have been determined 


The glasses were loaded in compression, and the 


amount of birefringence was measured with a polarimeter, using a quarter-wave plate com- 


pensator. 
2.95 and 
coefficient 
l 


1.82 brewsters 
with 
16 bre wasters, 


1. Introduction 


The double refraction produced in isotropic solids 
by the application of stress was discovered by 
Brewster [1]? who, incidentally, used glass for the 
greater part of his experiments. From Brewster’s 
early work, the theory of photoelasticity has been 
developed [2] and is the basis of a recognized method 
of stress analysis for isotropic transparent solids. 

For a given material, the relationship between the 
optical-path difference and the applied stress is 

p= Cid. (1) 
Cis called the relative stress-optical coefficient and 
is characteristic of the type of glass under investiga- 
tion. The above equation holds true for both com- 
pression and tension [3], and it is presumed to be 
valid up to the breaking stress of glass [4]. When 7, 
the retardation, is expressed in angstroms; ¢, the 
stress, in bars; and d, the thickness, in millimeters, 
the resulting coefficient, C, is measured in brewsters. 

The effect of the chemical composition of glass on 
the photoelastic effect was first investigated by 
Pockels [5], who came to the conclusion that an in- 
crease of lead oxide, of boric oxide, 
always lowers the relative stress-optical coefficient. 
In compression most glasses act as uniaxial negative 
ervstals, but| Pockels succeeded in observing a re- 
this effect for very heavy lead silicate 
The effect of the composition was more 
studied by Filon [6] and Adams and 


|, who reached the same general con- 


or decrease 


‘ 
< 


versal of 
glasses. 
accurately 
Williamson [7 
clusions. 

In more recent work, Vedam [8] measured the 
elastic and photoelastic properties of 18 optical 
glasses, and Schaefer and Nassenstein [9] made a 
very extensive study of these properties for 154 
optical glasses produced by Schott & Co. Unfor- 
tunately, the compositions of the glasses measured 
by Vedam are given only to a rough approximation, 
and these data are not mentioned in the paper of 
Schaefer and Nassenstein. 


‘his work was sponsored by the Air Research and Development Command 
U.S. Air Fores under Contract No. AF 18 (600)-139 
2 Figures in brackets indicate the literature references at the end of this paper 


weight pereent of lead oxide was found with values ranging from 
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The stress-optical coefficients determined for the crown glasses ranged between 
For the flint glasses, a regular variation of the stress-optical 


3.18 to 


As the relative stress-optical coefficient relates the 
stress exerted upon a glass to the optical-path differ- 
ence produced, a knowledge of this coefficient per- 
mits one to carry out a stress analysis of a glass speci- 
men on the basis of the optical evidence of the exist- 
ing strain. Furthermore, Lillie and Ritland [10] have 
shown that a knowledge of this coefficient for optical 
glasses is needed in order to develop fine annealing 
schedules that will keep the amount of birefringence 
within prescribed limits. In the course of developing 
annealing schedules for large optical elements, it was 
felt desirable to evaluate these coefficients for several 
glasses made at the Bureau. Because these are simi- 
lar to commercial types of optical glasses, the results 
are offered for their potential interest to glass tech- 


nologists. 
2. Optical System 


In most investigations of the relative stress-optical 
coefficient, the stressed specimen is placed in a 45-deg 
position between crossed nicols. A Babinet com- 
pensator, also placed between the crossed nicols and 
at the same angle as the specimen, is used to measure 
the optical-path difference. However, the use of a 
quarter-wave plate as a compensator affords greater 
accuracy, and was used in this investigation [11]. 
With this arrangement, the polarizer is set at a 45-deg 
position with respect to the principal directions of the 
stressed specimen, and the elliptically polarized light 
emerging from the glass is restored to a linear polari- 
zation by a properly oriented quarter-wave plate. 
The quarter-wave plate and stressed glass in com- 
bination then have the same effect as an optically 
active material that rotates the plane of polarization 
of the incident light. The angular difference in the 
extinction positions of the analyzer may be noted if 
the analyzing nicol is mounted in a graduated circular 
scale. The amount of rotation is equal to one-half 
the phase angle, and may be used to measure the 
optical-path difference in the birefringent glass. 

The optical arrangement and loading device used 
in this investigation are shown in figure 1. The light 
source consisted of a laboratory sodium are lamp. 
The light was collimated and filtered through a glass 
cell 15 mm in length that contained a 6-percent solu- 
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tion of potassium dichromate. The wavelength was 
then taken as the mean of the two sodium /) lines, 

5893 A. 

The polarizer, quarter-wave plate, and analyzer 
were each mounted in a brass tube 3 in. long and 1% 
in. in diameter. Nicol prisms were used for polarizer 
and analyzer, and the quarter-wave plate was a 
Babinet compensator adjusted to give the proper 
retardation. 

The circles for measuring the amount of rotation 
were 8, 1in. in diameter and graduated to 0.5 deg, with 
a vernier permitting readings to be made to 0.1 deg. 
A small telescope was mounted in the tube that held 
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the analyzing nicol and was focused upon a -in. 
diameter hole, which formed a stop for the light beam 
as it emerged from the polarizer. 


3. Preparation of Samples 


The optical glasses tested are shown in table 
together with their chemical compositions as com- 
puted from the raw materials in the batch. Column 

of table identifies the glass according to optical 
type, refractive index, and Abbe valve [12]. 

These glasses were made in clay pots holding from 
1,000 to 2,000 Ib of glass, depending on the type of 
A 5-lb piece of high quality glass,® which was 
representative of the type of glass under investiga- 
tion, was selected from each melt. A test specimen 
was cut from each piece of glass and annealed. Dur- 
ing heat treatment, each specimen was enclosed in a 
metal box in order to eliminate temperature gra- 
dients. The specimen was supported on small bits 
of refractory material so that it was not in contact 
with any part of the metal. The glasses were heated 
to a temperature about 20 deg C below the deforma- 
tion point [13], held at this temperature for several 
hours, and then cooled at a rate of 1 
through the transformation region. Once through 
the transformation region, the glass was rapidly 
cooled to room temperature. 

Each specimen of glass was annealed and _ then 
carefully ground to the shape of a rectanguls ar block, 


glass. 
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tention is called to the fact that the length of the 
specimen was made three times its width to avoid 
end effects [14]. 

After annealing, grinding, 
residual strain was measured. An optical path dif- 
ference of no more than 2 or 3 myu/cm was observed. 
This nominal residual strain was presumed to remain 
constant during the tests. <A zero correction was 
unnecessary as the relative stress-optical coefficient 
is determined from differences in applied load. 


and polishing, the 


4. Method of Loading 


The load was applied by means of a vertical shaft 
which supported a horizontal voke, and suspended 
loading platform, as shown in figure 2. The shaft 
and the two vertical guides were machine polished 
to insure a smooth and close fit, thus allowing only 
vertical motion to the shaft. <A *-in.-diameter steel 
ball was seated in the center of a metal disk at both 
ends of the glass specimen. Seats were also made 
on the end of the shaft and the bottom steel support 
to accommodate the ball. When this assembly was 
inserted under the shaft, the axis of the shaft and 
glass coincided. By taping the glass near both ends, 
the specimen was centrally located within each disk. 
A light film of oil was placed at the glass and metal 


interface. This method has been described in a 
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previous paper [15], and it has proved quite satis- 
factory in producing a uniform stress.‘ 

By means of leveling screws and a bubble level, the 
steel cylinder was adjusted into a vertical position, 
making it possible to rest the glass on a horizontal 
surface. The glass, being then in vertical position, 
could be placed in alinement with the direction of 
force created by the load. This line of force es- 
tablished the directions of the planes of polarization 
of light transmitted by the glass. 

To apply the stress, 50-lb weights were supplied 
by the Mass and Seale Section of the Bureau. Upon 
calibrating these weights, it was found that they were 
accurate to within 0.007 lb, which was felt to be more 
than satisfactory because it was desired to determine 
the values of the stress-optical coefficients only to 
three significant figures. 


5. Procedure 


Each part of the optical train was adjusted so as 
to be perpendicular to the light path. This was 
accomplished by noting the reflection of light from 
each surface and making certain that the reflected 
beam of light coincided with the incident beam. 

It was necessary also to have some base point on 
the graduated circle of the analyzer in order to 
orient the polarizer and quarter-wave plate in 45- 
deg positions relative to the optic axes of the stressed 
glass. For this reason each specimen was loaded to 
300 Ib, and one of the optic axes was found by 
alternately rotating the polarizer and analyzer until 
a position of minimum intensity was reached. Ex- 
tinction positions noted on the analyzer and polarizer 
could be reproduced to within +0.1 deg of are. 
Both axes were found in this manner, and a point 
midway between the two, as noted on the analyzer, 
was used as a base point for subsequent settings of 
the polarizer and quarter-wave plate. These ele- 
ments then could be set to within +0.1 deg. This 
method was felt to be satisfactory, considering the 
magnitude of errors that might possibly be intro- 
duced by missettings [11]. , 

The birefringence produced in stressed glass may 
be used to ascertain the uniformity of stress distri- 
bution, as well as its magnitude. Before making 
anv measurements, several of the specimens were 
loaded and examined in a broad beam of white 
light between sheets of Polaroid in a crossed position. 
With the aid of a sensitive red-tint plate it could be 
seen that the whole central portion of the specimen 
had a uniform color. 

In making measurements, the glass was weighted 
with a tare consisting of the steel shaft, voke, and 
loading pletform, and the extinction position of the 
analyzer was noted. Then a 50-lb weight was added 
to the loading platform, and another extinction 
position noted. This process was repeated with 
50-lb weight increments until the maximum load 
was 300 lb. Waiting about 5 min before readings, 

‘With this method of loading, the vertical thrust passes exactly centrally 
through the block, and it is a well-known principle in mechanics that at a suf 


ficient distance from the point of application of a load, the stress system depends 
only on the statical resultant of the load. 








the measurements were made in a darkened room 
so that the eve might be more sensitive to the in- 
tensity of the light beam viewed through the analyzer. 
The room was kept at a temperature of approx- 
imately 20° C. 

It was noted that a specimen would exhibit the 
same extinction position before and after the applica- 
tion of the 300-lb load. This was checked with 
several glasses, and no difference could be detected 
outside the error of the instrument. 

In investigating a particular specimen, 10 measure- 
ments were made for each extinction position, making 
a total of 70 measurements. 


6. Results and Discussion 


The results of this investigation are given in 
table 1, which lists the values of the relative stress- 
optical coefficient and their standard errors. The 
values were determined from the measurements in 
the following manner. 

If Aw is the load increment in pounds, d the thick- 
ness of the specimen in inches, and 6 the width of 
the specimen in inches, we obtain 


1.752C (Aw)d 
. bd 


~ 


where r is the retardation in angstroms, and ( is 
the stress-optical coefficient in  brewsters. The 
number, 1.752, is a conversion factor for the units 
emploved. 

As discussed in section 2, the retardation can also 
be expressed in terms of the angle through which the 
piane of polarization is rotated. The equation reads 


2(Ay)A 
360.0 


where 2Ay is twice the rotation of the analyzer in 
degrees of arc and X is the wavelength of the light 
in angstroms. 

The angular rotation of the plane-polarized light 
is a linear function of the applied load. The incre- 
ment in angular rotation can be represented by the 
relation 


Ay= B(Aw) (4 


from eq (2) and (3), where B is the slope of the 
straight line. It is appropriate to consider incre- 
ments rather than accumulated total rotations, 
because, for each glass, the measurements correspond- 
ing to the various loads were all made on the same 
specimen, and consequently a cumulative effect may 
be expected in the errors of the accumulated rotation 
values. Errors in the optical measurements, on the 
other hand, are effectively reduced through the com- 
pensating effect of averaging 10 duplicate measure- 
ments at each load. A detailed statistical analysis 
shows that this reduced error in the optical measure- 
ment is negligible in comparison with the experi- 
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mental errors due to the method of applying the 
load. Under these conditions, the correct least- 
squares estimate of slope B is obtained by averaging 
the five values of Ay, corresponding to the five con- 
secutive 50-lb increments. From these five values 
it is also possible to calculate the precision of the 
evaluated slope. Each slope estimate was converted 
into a C value, expressed in brewsters, and _ the 
standard error of © was derived from that of the 
corresponding slope. Table 1 lists the C values and 
their standard errors. A measure of uncertainty 
(95% confidence interval) can be obtained for these 
values by considering the interval C + (2.8 times 
standard error). 

As shown in table 1 most of the glasses studied in 
this investigation have values of stress-optical 
coefficient that are fairly close together, ranging 
from 2 to 3 brewsters. Two glasses in table 1, 
Nos. 16 and 17, vary only slightly in composition, 
and may be seen to have almost the same value of 
stress-optical coefficient, 1.85 and 1.82 brewsters, 
respectively. 

Because optical glasses are in general made up of so 
many different chemical components, it is difficult 
to correlate changes in physical properties with a 
particular constituent. As the series of flint glasses 
consists principally of silica and lead oxide, it is 
interesting to note that, when the values of the 
stress-optical coefficient are plotted against per- 
centage of lead oxide, a regular variation is indicated, 
This confirms the observations of other investigators, 
notably Pockels® [5] and Filon [6], whose results are 
shown in figure 3, together with the values obtained 
in this investigation. Furthermore, flint glass 
920/211 is seen to depart from the general rule that 
glass in compression acts as a uniaxial negative 
crystal. 

Pockels determined the absolute changes in the refractive index of a stressed 
glass for light polarized, parallel, and perpendicular to the direction of force, and 
expressed his results by the coefficients p/r and q/r. With knowledge of these 
values, the relative stress-optica! coefficient, C, may be calculated for purposes 
of comparison. This was done by Filon, and his results are shown in figure 3, 


Adams and Williamson also calculated C for Pockels’ glasses and found that the 
results showed good agreement with their own determinations. 


| ONBS GLAsses ___| 
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Pockels [5] and Filon [6] also concluded that | 


lowering the amount of B,QO, in a series of glasses 
decreased the relative stress-optical coefficient. It 
was impossible to verify this from the results of the 


present investigation because the glasses studied 
here were multicomponent and contained com- 


paratively little boric oxide, the amount being, at 
most, only 12.5 percent. 


The authors are deeply appreciative of the statis- 
tical analysis of the data by John Mandel and 
Mary N. Steel. They are also grateful to Joseph 
\. Nivert, Jr., for preparation of the glass specimens 
for this work. 
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Physicochemical Studies of the Destructive Alkali- 
Aggregate Reaction in Concrete 


Robert G. Pike and Donald Hubbard 


In a further study of the 


alkali-aggregate 
and nonreactive experimental aggregates were examined by several procedures. 


reaction in high-alkali cements, 11 reactive 
On com- 


paring the results of each method with the standard expansion-bar test, it appears that each 
property studied has some bearing on the expansion characteristics of each aggregate. For 


example: (1) The chemical durability 


poor for the 


aggregates that showed excessive 
shown to be chemically reactive by this procedure will cause expansion, 


as determined by the interferometer procedure was 


but not all aggregates that are 
(2) all aggregates 


expansion, 


are sufficiently hygroscopic to attract the water necessary for the reaction if other properties 


are also favorable, (3 
determinations; however, free 


surface electrical-resistance measurements confirm the hygroscopicity 
ions in glasses may have considerably more effect on the sur- 


face resistance of certain glasses than do their hygroscopic properties, and (4) the uneven 


distribution of migratable ions between the 


mechanism for the osmotic 
expansion in concrete members 


1. Introduction 


In a previous publication [1]! concerning the 
mechanism of the destructive expansion of concrete 
structures caused by opaline and glassy silicate types 
of materials when used with high-alkali cements 
IZ, Ds 4] there were included some data on the chemi- 
eal durability of a few natural and synthetic aggre- 
Those results, obtained with an interfero- 
metric procedure under controlled exposures of time, 
temperature, and pH, illustrated surface alterations 
of the aggregate materials during the incipient stages 
of attack. In light of the varied and unexpected 
dimensional changes found for various reactive and 
nonreactive aggregates, it became apparent that 
much remained to be learned about the chemical and 
physical nature of these and kindred materials before 
any completely acceptable picture of the destructive 
alkali-aggregate reaction could be formulated. 

The present investigation was undertaken to gain 
further insight into the physicochemical principles 
operating in the alkali-aggregate reaction by extend- 
ing the previous work on chemical durability to a 
broader range of aggregate materials. Further, in 
order to obtain a more comprehensive picture of the 


grates. 


nature of these materials, their relative hygro- 
scopicities, heterogeneous equilibria, electrical sur- 
face resistance, and reactivity in expansion bars 


were studied. 


2. Experimental Techniques and Results 


2.1. Techniques 


For the sake of clarity, the tee hniques for each of 
the procedures are outlined briefly or covered by 
references in the sections where the results for each 
of the respective methods are presented and dis- 
cussed, 
it the end 
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2.2. Expansion Bars 


The expansion of standard mortar test bars 
[5, 6]* brought about by 11 different aggregates 
when used in a high-alkali cement is shown in figure 
1. The three forms of silica, viz, quartz (crystalline 


The expansion bars were prepared and measured according to the procedure 
of Mielenz and Witts [5]. A water to cement to standard Ottawa sand ratio of 
0.4 to 1 to 2 was used, and the bars were then stored over water at 100° F, as 
described in a previous paper [1]. Six percent of Ottawa sand was replaced by the 
aggregate under test, using the size fraction passing a No. 50 and retained on a 
No. 100 U. 8. Standard Sieve. This size gave maximum expansion with Pyrex, 
ind 6 percent seemed a fair compromise between the expansion caused by the 
addition of 8 percent of Pyrex, as reported previously [1], and that caused by 
lower percentages of opal, as reported by other investigators [4, 8]. An analysis 
of the cement used is given in table 1. However, it should be kept in mind that 
the curves in figure 1 would be altered if the test conditions, such as percentage 
4) or size [8] of aggregate, temperature of storage [9], and percentage of alkali 
[10], were changed. 
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SiO,), fused SiO, (glass), and opal (hydrous SiQ,), 
form an interesting series, ranging in activity from 
quartz with no abnormal expansion of the bars in 8 
months, to opal with approximately 2 percent of 
expansion in 4 months, followed by mechanical 
failure of the specimens. Fused SiO, occupies an 
intermediate position, with an expansion far in 
excess of permissible tolerance [5]. Other aggregates 
besides quartz causing no appreciable expansion are 
F 40, BaC 572, and fiber glass type E. Pyrex 
7740, the glass adopted as a standard reactive 
aggregate [6], is included for comparison. 

The induction period, which appears in the curves 
for the reactive aggregates, can be partially attributed 
to the shrinkage that takes place during the early 
stages of the setting of cements. However, the 
great difference in the induction periods for different 
materials implies that the induction features are 
primarily characteristic of the individual aggregates. 
These composite characteristics are undoubtedly a 
function of the chemical durability and other 
chemical and physical properties of the materials 
and the storage environment. 


2.3. Chemical Durability * 


The varied chemical durabilities over an extended 
pH range, characteristic of the 11 aggregates evalu- 


ated in figure 1, are shown in figures 2, 3, 4, and 5 
The chemical-durability curves for the three 


forms of SiO, (quartz, fused SiO,., and opal) are 
shown in figure 2. These three forms of silica have 
distinctive individual durability characteristics, with 
quartz showing no detectable alteration between 
pH 2 and pH 11.8. Over the same range, fused 
silica shows surface alteration (attack) only at 
alkalinities above pH 10, whereas opal exhibits 
dimensional changes over the entire pH range, with 
liberal swelling in the acid buffers and severe attack 
at pH 11.8. The relative order for destructive 


swelling of mortar bars by these three forms of 


relative chemical 
and opal the 


their 
the least 


as for 
quartz 


silica is the same 
durabilities, with 
most reactive. 

In figure 3 are grouped the silicate aggregates 
that exhibit pronounced swelling over the acid range 
and well into the alkaline region, similar to opal. 
These three greatly dissimilar silicates (table | 
having the common characteristic of swelling, also 
cause conspicuous destructive expansion of mortar 
bars, as shown in figure 1. Continuing, in figure 4 
are presented the glassy aggregates that showed 
unusual attack near neutrality by the buffers between 
pH 4 and pH 10, in contrast to the aggregates shown 
in figure 3, with extensive swelling over the same 
pH range. A glance at figure 1 reveals a pronounced 
difference for these groups in the expansion-bar test. 
For example, of those with the attack near pH 7, 
one (BaC 572) shows no destructive expansion, and 

The chemical durability of the various in 
interferometer procedure that has been particularly useful in the study of the 
chemical durability of optical glasses, voltage anomalies of the glass electrode, 
und in practical applications in various other fields [11, 12 The procedure 


consists in observing the surface alteration of the specimen after exposure to 
Britton-Robinson universal buffer mixtures [13] for 6 hr at 80° C 


iggregates was determined by 
f 


128 


FIGURE 


FIGURE 





‘T = T Tv T T 7T 5 | 
} 
4 
| 
| 
j 
oa | 
3} 
| ) 
r 
| 
| } 
| N 
i iN 
“ | Ca 
oO 2r Ba 
2 | pt 
= | | / 
we 
oe Zn 
W 4 
2 B; 
S as 
= r st 
< ss! 
a Si 
Ww 
e 
- 4! 
=z} | Fi 
‘ M 
w £ - 
> | + USED ] ae 
<= «| SILICA 
uo! I 
@ s} 4 i 
25] 
wal _*f QUARTZ + 
ee I 
r) } 
i 1 H 
=| . 
4 7 
, | Cc 
r 1 Cc 
| Publi 
} ~ “ 4 | : 

a | jecter 
| | sed 
a a ee ee eee ee reel 

DH preci} 
2 Chemical durability (surface alteration ol three 
forms of silica. 
Exposure 6 hr at 80°C. 
' 
7 —~ _ 
f Q 
| i 
5 j : 
| 
3} 
} 
} 
} 
| 
2 
7) | 
ny 
So Tf 
= | RNIN 
« | $ 
= | 
z | 
- 
da 
a 
= } 
a / 
<4 
ao 
1S) “4 ; 
a 
+ ‘| 
~ | 
Y o oe — —— ray 
zt RNIN s 
=} 
y 
rT j 
| PAL 
b oO 
SODA FLINT 
r & 5 —<tr 4 
ze ae =a J Fre 
9) 2 a 6 a 4 ‘ 
pr } 
3 Chemical durability surface alleration of three D 
aqgreqates that exhibit swelling in the acid ranqe 


Exp 


u 


re6Ohrat soc 


TaBLe 1. Chemical compositions | the other two (BSC 517 and F 620) much reduced 
| or delayed action. 
; An additional group of experimental aggregates is 
given in figure 5, in which the chemical durability of 


O15 
la 
Kk 





Oxides ei6izi+ idle leal.|. = Pyrex 7740 glass, which serves as a standard reactive 

} Fig i$ /8/\4\ & \4"|9/2/ 27) =° aggregate for testing high-alkali cements [6], is com- 
sail Hhendll adil Tradl Ihodll Hadi soc head! Realll ies pared with F 40 and fiber glass (commercial type E). 

e igi eae | The latter two aggregates show attack in the acid 

KO. me 271s BB | Od io» | region at pH 2, but are dissimilar in their reaction 

Cad "lon 4 ane =. in the alkaline range, with F 40 being unaltered, 

PbO 45.1 20 whereas the fiber glass is severely attacked. These 
an0 7.2 05 | two aggregates, along with quartz and BaC 572, 

BiOs thier sy | showed no destructive expansion of the test mortar 

Sh:0 sal aos | aac alaical eal ocala |" “soleaconleargal ace bars. As F 40 and quartz are chemically resistant 

- : wie Be Se | to attack by the high-alkaline buffers, their accept- 

AO: 22 145 "lt os 23 | | able behavior in the mortar-bar tests is to be ex- 

MgO 5.0 19 | pected. The behavior of BaC 572 and the fiber- 

TO 20 ~" | glass aggregates, with their poor chemical durability 
| itil = 40 532 11 | In the alkaline range, seems to be anomalous. How- 
Insoluble > | ever, these two exceptions furnish the most inter- 
HF residue « 151 430 —«| esting and probably the most useful data of the 

| entire series. For example, cross sections of the 

Composition from references [1,1] » Compositions of batch material. | expansion bars reveal that BaC 572 and fiber-glass 
and “lee tek TE 4 Commie han ee ee aggregates had not been attacked by the high- 
ct Ce een has. t The HF residue was fused and dissolved and'k.o, | @lKaline conditions [14] existing in the hydrated 


eeipitated with NH,OH and ignited. ‘The FezOs was determined in the ignited | high-alkali cement. Thus, the inhibiting effect of 

recipitate and the remainder of the RO; assumed to be AhO ‘ ~ ° 
—_ ae ” saturated Ca (OH), [15] from the cement, and certain 
individual constituents of the aggregates, such as 
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2.4. Hygroscopicity * 


As the alkali-aggregate reaction cannot take place 
without the presence of water, and also because 
water is one of the chief constituents in the products 
resulting from the reactions, the relative hygro- 
scopicities of the aggregate materials were deter- 
mined. Figure 6 illustrates the results obtained for 
the 11 aggregates of figure 1, after 1 and 2 hr of 
exposure. Some of the groupings applicable to the 
chemical-durability data persist for the hygro- 
scopicity in a superficial fashion. For example, the 
three forms of SiO, fall in the order, opal, fused 
Si0,, and quartz, although the differentiation 
between fused SiO, and quartz is not very pro- 
nounced. The three chemically dissimilar aggre- 
gates, opal, Corning 015, and soda flint, characterized 
by swelling in the acid range, are also the three most 
hygroscopic materials of figure 6. Likewise, 6 of 
the 7 aggregates that caused unacceptable destruc- 
tive expansion of mortar bars are found to be the 
most hygroscopic of the 11 materials tested, fused 
silica being the only destructive aggregate that does 
not conform. BaC fiber glass, quartz, and 
F 40, the 4 nondestructive aggregates, are among 
the 5 least hygroscopic materials. The fact that 
fused SiO, constitutes such a glaring exception once 


572, 





again furnishes the most interesting data of the 
series and needs further study. 
‘ The procedure for determining hygroscopicity was that reported in a previous 
publication [1], which consists in determining the change in weight of sample 
water sorbed [16], ground to pass a No. 140 U.S. Standard Sieve) when exposed 
to the high humidity maintained by a saturated aqueous solution of CaSO, 2H,0 
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2.5. Heterogeneous Equilibria 


The uneven distribution of migratable ions, ag 
shown in figure 7, using excess Ag(NHs,).* as the 
indicator, reveals that the aggregate materials studied 
possess a negative nonmigratable ionic charge after 
leaching. Early inquiry into the destructive expan. 
sion of concrete structures caused by the presence of 
reactive aggregates, when used with high-alkalj 
cement, revealed that osmotic forces were operatiye 
(3, 17]. In addition to differentiating between aggre. 
gate materials, these data also demonstrate that the 
frameworks of these materials (primarily silicates 
behave as negatively charged nonmigratable ions, 
This condition is favorable for initiating and aceel- 
erating reaction within the silicate framework. The 
increased sodium ion concentration within the silj- 
cate framework over that in the ambient phase 
would reduce the solubility of Ca(OH), so that it 
precipitates at the aggregate surface [14], inhibiting 
reaction there while the inner portion of the particle 
was being completely reacted and disintegrated as 
alkali silicates [1, fig. 14]. The destructive capacity 
of any aggregate would be partially determined by 
this heterogeneous equilibrative factor and by the 
chemical and physical nature of the reaction products 

The three forms of SiO, fall in their usual order of 
reactive and destructive capacity, with opal occupy- 
ing the least, and quartz the most favorable position 
of the series. 

In the experiments on heterogeneous equilibria, the uneven distributior 
the migratable ions Ag(N H,).* and Br- between the aggregate material and t! 
ambient phase was determined by potentiometric titration of the excess Ag 
ons remaining with the aggregate or the excess Br- ions appearing in the ambient 


solution after the leached aggregate sample was treated with full-strength aqueou 
N H,OH saturated with AgBr. The procedure consisted in treating 10 g of aggre 


gate material (—5O0-, +140-mesh particle size) with %) ml of the Ag(NH,))B 
solution for 5 minutes, decanting and acidifying, followed by potentiometr 
titration, using Ag-AgBr and glass electrodes as the indicator and reference 
electrode respectively Standardized solutions of KBr or AgNO whicheve 


the situation required) were used in the titration 
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2.6. Surface Resistance ‘ 


Figure 8B shows the surface hearse resistance 
obte ained for the same five as gveregates whose hetero- 
geneous e quilibri ium characte risties are given in figure 
- The surface resistances for the three forms of 
silica constitute a consistent series comparable with 
the heterogeneous equilibrium data, but the positions 
of the opal and Corning 015 aggregates are reversed. 
They follow, instead, the order of their hygroscopicity 
characteristics in figure 6. Thus, it appears that elec- 
trical surface-resistance Measurements mi: L\ give an 
indication of hygroscopicity. However, there are 
marked differences in the characteristics of the curves 
for these two types of data. The electrical-resistance 
data come to apparent acceptable equilibrium values 
for periods exceeding 2 hr, whereas the hygroscopicity 
curves (figs. SA and 6) In some cases, hygro- 
scopicity values have continued to increase with the 
vapor until the weighing bottles 


do not. 


acquisition of water 
have been filled to overflowing, as dictated by the 
lowered vapor pressure of aqueous solutions. On the 
other band, electrical-resistance measurements give 


——_——— 


The surface-resistance data were obtained by measuring with an eight-decade 
rithmic ohmmeter [18] the resistance of a cell (fig. 10) containing a measured 
lume of powdered sample 30, +50 U.S. Standard Sieve) when the cell was 
placed in the grid circuit and exposed to a controlled relative humidity. This 
cial meter Was necessary because the usual d-c megohm bridge caused excessive 
larization, and the readings by the available a-c meters were vitiated by the 
city of the cell, 
x 
so 2 
4 
7 
q 
Fi RE A Electrical existance and water sorption of several 


aggregates 


\, Typical hygroscopicit y curves illustrating continued sorption of water with 
! 4 exposure; B, typical surface electrical resistance curves for five aggregates 
‘posed to 9S percent relative humidity at flow rate of 0.5 liter per minute, illus- 


n equilibrium state after 2 hr 


constant values characteristic of the aggregate mate- 
rial.’ 

Although the surface ohmic resistance and hygro- 
scopicity “of the aggregates might be expected to 
parallel each other, inspection of ‘the results in figures 
6 and 9 reveal reversals of position for several of the 
aggregates, and some of these even occupy different 
positions at different humidities. These latter re- 
versals are undoubtedly determined by the nature 
of the hydrated and hydrolyzed reaction products 
following the sorption of moisture. 


? This constant resistance can be attained much more rapidly at faster flow 


rates through the resistance cell. 
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FIGURE 9 Surface electrical resistance of 11 reactive and non- 


reaclsve aggreqales al various relative humidities 
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Figure 10. Diagram of the cell used in measuring the resistance 
of aggregate samples while being exposed to controlled relative 
humidities. 
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3. Conclusions 


Although the properties studied and the procedures 
used lead to a broader understanding of the aggre- 
gate materials and of the part played by various 
physicochemical principles operative in the destruc- 
tive expansion of concrete structures, it is evident 
that no single property was determined that can be 
used as a serviceability criterion for all materials. 
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675 Insulated Loop Antenna Immersed in a 





-alkal . . 
a Conducting Medium 
» iy | 
oS James R. Wait 
Sci. 2, ] 
roper. \ solution is given for the fields of a circular loop in a conducting medium. The loop 
earch is assumed to have a uniform current, and it is enclosed by a spherical insulating cavity. 

The impedance of the loop is also considered. It is shown that the power radiated from the 
meter loop varies approximately as the reciprocal of the radius of the cavity for a specified loop 
rabil. current. Furthermore, if the cavity is electrically small, relative to the external medium, 
NBS | the radiation field is not significantly affected by the presence of the cavity. 
. van | 
ning, | 1. Introduction 
‘ining 
1943 
udies There has been some interest shown in recent vears on the subject of radiating systems 
yp immersed in dissipative media [1 to 9].!. Such a condition exists when transmitting antennas 
ptio are located in sea water, buried beneath the ground, or surrounded by an ionized medium, such 
» 916 as the ionosphere. Most investigators have considered only the propagation aspects of the 
K.T general problem. That is, attention was devoted only to the radiated field for a specified 
rete, antenna current. An associated, and a usually much more difficult problem, is the calculation 
ohm. | of the impedance of the antenna. 

Tai {1}, in an interesting analysis, showed that for a Hertzian electric dipole in a dissipative 
medium, infinite power is required to produce a finite field strength at some distant point in the 
medium. This difficulty is overcome when the dipole is placed in an insulated cavity. It has 
also been shown that a magnetic-dipole-type antenna was more efficient than an electric- 
dipole-tvpe antenna [2, 5]. 

It is the purpose of the present paper to reinvestigate the problem treated before. Spe- 
cifically, the impedance and fields are calculated for a circular loop with a spherical insulating 
cavity, all immersed in a homogeneous dissipative medium. In the previous analysis [5], the 
loop was represented as a magnetic dipole at the center of the cavity. It was implied that the 
dimensions of the loop were small compared to those of the cavity. In this sequel, these limi- 
tations are removed by treating the problem from another viewpoint. 

In terms of spherical coordinates (7, 6, @) the insulating cavity is defined by r<a@ and has a 
dielectric constant, «. The exterior homogeneous dissipative medium, r>a, has a dielectric 
constant, e, and a conductivity, ¢. The permeability, u, of both regions is taken to be that of 
free space. The loop, for the moment, is considered to be a single turn carrying a uniform 
current, 7. The coordinates of the loop are @=8 and r=), so that the axis of the (circular) loop 

passes through the center of the spherical cavity (see fig. 1). 

The dimensions of the cavity are taken to be very small compared to the wavelength, 
which suggests that a quasi-static approach to the problem should be adopted. As a conse- 
quence, the fields are a solution of Laplace’s equation inside the cavity and a solution of the 
Wave equation outside. Although it is not necessary to invoke the quasi-static assumption for 
the field calculation, it leads to the desired result in a direct fashion without becoming en- 
cumbered with a myriad of spherical Bessel functions. Furthermore, in any instance where 
the cavity was not small compared to the wavelength, the assumption of uniform current around 

: the loop would be a poor one for a localized generator. 
q Figures in brackets indicate the literature references at the end of this paper, 
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2. Formal Solution 


In view of the polar symmetry of the stated problem, the usual vector potential, .1, has 


only an azimuthal component, A. The nonvanishing field components are then given by 


kk ipwAA, 
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The (quasi-static) form of the vector potential of the primary field of a circular loop_in an 
insulating space can be expressed in terms of spherical harmonics (see appendix 


/ sin BS, bh b\’ 

S ( )P cos 6)P' (cos B 2 
9 hand ; 

2 pert 101 L)\7 

for r >b, where 7’) is the associated Legendre function To account for the spherical wall of 
the cavity, a secondary field, A’, must be added. Because it is also a solution of Laplace's 
equation and is finite ;=0, it must be of the form 


T sin BY . 8&8 7 
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} ) P' cos 6)P! (eos B) 3 
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2 n= iii 


for r<a, where S, ts an unknown coefficient 
The field in the external region is a solution of the wave equation, and furthermore the 

solution must give rise to outgoing waves at infinitv. Therefore, 

T sin Bey Tb kily 
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for r >a, where k,/2) is a spherical Hankel function defined by 


n (n+m)! 
~ ~ 
(5) 


k,(z) =e = 
a 


—m'(n — m)!(22) 


y= ((uow— eno*)’®, and 7), is an unknown coefficient. The boundary condition at r=a requires 


the continuity of #, and 77s. The solution of the resulting algebraic equations leads to 


n+a, b\2*+! ‘ 
S, ) (6) 
(n+1)—a, f\a 
_ 4 (2n+1) b\” on 
Ti ) ; (7) 
(n+l)—a, \\a 
d 
Q, f log k, | 2 | 
dz _ 


where 





ya 
When S, and 7), are inserted into eq (3) and (4), the complete solution is obtained. 


3. Discussion of Result 


It is now of interest to investigate the influence of the cavity on the external fields. With 


this in mind, 1. for ra. is written in the form 


ITsn Be, 6 k(yr)n! 
rA S = (2yb)"G, « P! (cos 0)P! (cos 6), (8) 
2 <in(n+1) (2n)! 7 
where 

9 Dy! 

' 2n+1 (2n)! 
G, . (9) 

(n+1)—a,,_] n!(2ya)"k, (ya) 

It can be readily verified that lim @,=1. When this limiting form is inserted into eq (8), 


~a—) 
the resulting expression would be the field of a small loop in an infinite dissipative medium 


without any cavity. Specific values of @, for finite ya are 


. oe 
Gy, = = ‘ (10) 
d+oste" 
and 
. L5e 
Gr — — ars ’ (11) 
15+152+62°+2° 
where z=ya. The value G, can be regarded as the ratio of the effective dipole moment for a 
loop in an insulated cavity to that of the loop without any cavity. It agrees with a result 


given previously [5] by a somewhat different method. @, is the corresponding ratio for the 
next higher multipole term. For most cases of practical interest, ya is somewhat less than 
unity, and the multipole terms are not important in the external region when yr is large com- 


pared to one. Furthermore, in this case 


and the external fields for a specified loop current are not dependent on the presence of the 


cavity. 
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Attention is now turned to the effect of the cavity on the impedance of the circular loop. 
The impedance of the loop in free space is denoted Z), which is assumed known. The im- 
pedance of the loop inside the cavity is then written 


Z=2,+4Z 


where AZ is the incremental change due to the finite size of the cavity. AZ is the secondary 
electric field “3 integrated along the circumference for a unit loop current. Therefore, 


»y) } 
, , 270 
AZ mwA’ ] (12) 
. . ‘ ff ei >I >I 
byw Sin B(rb) >> ( ) P'(cos B)P\(cos 86), (13) 
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where 7, is given by eq (7). The above expression is simplified somewhat in the important 
case where |ya! is small compared to unity, because then 


) 
T,.=— = (Ya) (14) 
(2n+-1)(2n 1) 
neglecting terms containing (ya)*, (ya)', ete. Furthermore, if @=8=7/2, it follows that 
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where S=rb*?. When displacement currents are negligible (i. e., w<e) and (b/a)?*<1, 


(uw)? S? 


AZ~ (16 


Oma 
For a specified current, J, on the loop, the power supplied to the surrounding conducting medium 
is then given by 


ya (ues)? 12S? 


, l2ra - 
which agrees with a previous result [5] derived from energy considerations. 

Although the preceding analysis refers to a single-turn loop, the results are directly appli- 
cable to a loop of N turns by replacing Sby SN. The appropriate expression must also be used 
for Zp, of course. 

It would appear that power radiated from the loop for a specified current varies approxi- 
mately as the reciprocal of the diameter of the cavity. This conclusion is not modified to any 
extent by the finite size of the loop in the cavity. 


4. Appendix 


It is well known that the vector potential can be obtained from a volume integration over 
the contained current. Furthermore, because the loop has a filamental uniform current with 
a component only in the ¢ direction, the vector potential has polar symmetry and has only a 
¢@ component A; it is given by the line integral 


A ah —_" b sin Bdd,, (18 
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where ¢ is the azimuthal coordinate for a point on the loop relative to the reference plane, and 
where 
R?—=r*+-b?—2rb cos Q 


and cos Q=cos @ cos B+ sin @sin 8 cos ¢. Now the factor 1/2 is a solution of Laplace’s equa- 
tion and can be written in terms of spherical harmonics, as is well known [10]. 
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pP = 2 Em, ; P™(cos 0=)P"(cos B) - cos md, (19) 
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i n= n / r” hr Ae r<b, 


where @=1, €,=—2 (m0), and where the P? are the associated Legendre polynomials. On 
inserting the above expression for 1/R, and utilizing orthogonality, eq (18) reduces to 
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which is the desired result. 
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Noncrystal Ionic Model for Silica Glass 
Leroy W. Tilton 


A model was constructed of pentagonal rings of symmetrical tetrahedra so connected 
that all five silicons lie in one plane and all angles Si—O—Si approximate 180°. From an 
initial tetrahedron the rings extend in six directions to include its six edges, and in the 
resulting three-dimensional network each ring forms a common interface between dodeca- 
hedral cages which necessarily have fivefold symmetry. Such symmetry is incompatible 
with formation of a crystal lattice and thus the structure can be extended in three dimensions 
only to a very limited degree and with unavoidable stress and distortional strain. A stressed 
network cluster of such dodecahedral cages is here called a viTRON. 

If the scale in this model for silica glass is taken as 1.6 A, from silicon to oxygen, the 
model is found to be in accord with other radial distances computed from diffraction data 
and the density is correct within 10 percent. The apertures of the pentagonal interfaces 
are of suitable size to explain data on the diffusion of noble gases through silica. Shared 
penta faces and intermittent oxygen bridges connect neighboring vitrons and constitute 
an interstitial tissue of relatively weakened structure that accounts for the observed low 
tensile strengths of silica glass and provides channels for viscous flow. 

The proposed cages would be less pliable at low temperatures and so the interstitial 
tissues between vitrons should expand in volume on cooling. This “negative’’ expansion 
would counter the normal expansion and account for the known extremely small (net) 
expansivity of silica glass. The accompanying changes in strength and volume of the glass 
at its weakest places, the interstitial tissues, can account for a number of known “anomalies” 
in thermal behavior of silica glass with respect to its volume, compressibility, elastic moduli. 
and viscous flow 

The vitron concept reconciles the crystallite and the network theories of glass by pro- 
posing nuclei that cannot grow extensively and a distribution of localized stresses that may 
constitute an acceptable degree of randomness and provide effective modulations in con- 
tinuitv. This concept of definite but limited microregularity in structure suggests that 
other nonerystal symmetries should be studied for possible interest in other fields such as 
liquids and high polymers. 


has been stimulating and fruitful in the field of in- 
organic glasses. A very limited degree of short- 
range order of arrangement of nearest neighbors and 
next nearest neighbors among the ions in glass is 
now definitely determinable by X-ray and electron 
diffraction, but the unbroken continuity and the 
very considerable randomness that is often assumed 
at longer distances seem incompatible with some 
experimental data. 

W. A. Weyl [5] has considered that the network 
theory “in its present form” is not satisfactory to 
account for some existing data, and he wrote a 
paper ‘‘to show the need for a revision in our present 
approach to the constitution of glasses.’’ In this 
paper he says: “One of the most important proper- 


1. Introduction 


Evidence from X-ray data shows that glasses and 
even liquids have some degree of ordered arrange- 
ment, although it definitely differs from that in 
crystals. Also, in studies of liquids and plasties it 
is recognized that limited ordered arrangements, 
often called “crystallinity,” can exist without forma- 
tion of macroscopic crystals. Under these circum- 
stances two important theories of the glassy state 
have developed. The erystallite theory was stated 
by Randall and coworkers [1]! in 1930. The 
presence of submicroscopic crystallites or crystalline 
nuclei in very large numbers in undercooled glasses 
was inferred, soon after the X-ray reports published 


| 


by Clark and coworkers [2] in 1929, in order to 
account for the breadth of the observed diffraction 
rings. The acceptance of this idea proved difficult 
because many investigators found no evidence of a 
normal growth in the size of such crystallites under 
conditions that do promote growth whenever nuclei 
of crystals are known to be present. 

The network theory resulted largely from the work 
of Zachariasen [3] and his coworkers in the field of 
chemical crystallography, and some associated ideas 
were more definitely formulated by Warren [4] and 
his coworkers in analyzing their X-ray data to estab- 
lish interatomic distances and in interpreting the 
breadth of the diffraction rings from glass as con- 
trasted with the well-defined rings from crystals. 
This continuous random network theory of glasses 


br ndicate the literature references at the end of this paper 


ties of glass, its mechanical strength, seems to resist 
all attempts to unveil its secrets’; again, ““‘We have 
no plausible explanation vet for the thermal expan- 
sion of glasses as a function of their composition”; 
also, “It is customary to talk about ‘flow units’ in 
describing viscosity about association-disasso- 
ciation equilibria in any discussion of the effects 
which thermal history has upon the properties of 
Nobody, however, has yet given an 
accurate description of a ‘flow unit’ nor pre- 
sented a reasonable concept of the unit which 
disassociates, he 

Confining attention to the simple case of silica 
glass, one finds numerous and well-established data 
that characterize this material as unusual or even 
unique in respect to some of its properties, such as 
expansibility, compressibility, rigidity, and varia- 
tions therein as effected by temperature and pres- 


rls _ 
glass. 
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sure. These relationships have been discussed in 


detail by Babcock, Barber, and Fajans [6], who say, loca 
“ _ . Use of the random network theory . . . has in z 
general been unsuccessful in interpreting the proper- vit 
ties of silicate glasses... .”’ Spinner [7], who meas- regu 
ured elastic moduli at elevated temperatures, says ; The 
“It seems clear that . . . the usual picture of the dist 
structure of silica, as a random 3-dimensional net- qua 
work of SiO, tetrahedra, needs to be amended or heat 
modified to explain the various phenomena asso- silic 
ciated with the material.”’ con 
Evidence of somewhat long-range order in silica | pro| 
glass was found by Melkonian [8], who wrote, fere 
“Fused quartz . . . was expected not to show micro- for 
crystalline interference effects because of its supposed | = 
lack of long-range crystalline structure. Thus, it a 
was expected to show a transmission nearly inde- } gis 
pendent of neutron energy and hence to be suitable 
for use as a secondary standard . . .. Fairly large 
interference effects are evident which make this ; Ss 
particular sample useless as a secondary standard.”’ pers 
In accord with this, Milligan, Levy, and Peterson | tetr 
[9], in their neutron-diffraction studies on vitreous con 


silica, found 5 or 6 well-defined maxima and 4 or 5 by 


View along fivefold axis of a single-element vitror 





others that were less well defined. Hoffman and | Ficure ! j asst 
. . _ struct f 20 tetrahedra 1 resenting 20 si 12 af © : 

Statton [10] used the low-angle scattering of X-ravs constructed of 20 tetrahedra representing 20 icon atoms at ain 
= » ‘ - . . . their centers and 30 orygen atoms at their connected corners 

in glass and found evidence of spacing in fused bon 


ro . _ ¢ . . l'wenty outer corners are available for attachments Only one-quarter . 
silica of 30 A or larger. In brief, there is ey idence each Si and one-third of each (connecting) O belong to this cage as an element o if al 


unit of a larger vitron In other words, 5 (SiOe) constitute a single element 


that the short-range order should be extended at | Unit of @ larger vitron. in other wor } can 
the expense of longer-range randomness. thes 
In order to study possible extensions of the fort 
accepted short-range order in simple silica glass, a Pon 
three-dimensional model was constructed starting tetr 
with silicon-oxygen tetrahedra connected at their the 
corners by the sharing of their oxygens. This ( 
cee 


model, described in this paper, gives the correct 
density for silica glass within 10 percent, provided 
the scale is determined in accord with the established 
average distance of 1.6 A from silicon to oxygen as 
found [11] for silica and the silicates. It will be 
shown that the interatomic distances are in very 
good agreement with all radial distances computed 
from X-ray, electron, and neutron-diffraction data, 
and the apertures are of suitable size to explain 
Norton’s [12] data on the diffusion of noble gases 
through silica. 

The proposed model is a network of pentagonal 
rings of tetrahedra with 12-sided regular cavities, 
as in figure 1 (approximate diameter, locus of silicon 
centers, 8 A). It cannot grow extensively without 
peripheral tensional distortions that ultimately pre- 
vent further maintenance of the silicon-oxygen 
bonds. A cluster of such dodecahedral cavities or 
cages, such as shown in figure 2, is here called a 
VITRON, & new concept that may have application in 
other fields such as rheology and high-polymer 
studies where noncrystal order may be evidenced. | 





Aside from the fundamental conceptual difference 
between a vitron and a crystallite, the proposed | Freure 2 
model does in many respects accord well with the 


2 Off-aris view of completed four-element  vitron 


type ,. tetrahedral 


aiid . : 1 .. - a . Coiled wire springs indicate extra O— charges and vacant tetrahedral corners 
cry stallite hypothesis. It IS also essentially the net- indicate Si+charges available for attachments to other vitrons by oxygen bridges 

, «Le . cid ; : sin a - At center foreground is a reentrant three-face that can unite with another vitror 
work hypothesis with a modulation of its continuity, | presenting s conver threeiace. See appendix 
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a reduced randomness, and an added feature of 
localized stress distribution. 

Partial unions and attachments between adjacent 
yitrons form interstitial channel tissues that are less 
regular and weaker than the intravitron structure. 
These ti together with the Jocalized 
distributions within vitrons, are important in 
qualitative explanations of the peculiar effects of 
heat and pressure on the volume characteristics of 
silica glass; also in understanding viscous flow and 
eonductance. Preliminary studies indicate that the 
proposed model will permit semiquantitative in- 
ferences concerning possible cavity-saturation limits 
for nonsilica that correlate with property-composi- 
tion curves of solubility, chemical attack, volatiliza- 
tion, and specific volume for binary alkali-silicate 


tissues, stress 


glasses. 


2. Pentagonal Dodecahedral Model 


Silica glass is here considered to have a high-tem- 
perature arrangement of connected silicon-oxygen 
tetrahedra that has been cooled so quickly that the 
conditions at high temperature have been “frozen”’ 
by rapid increases in rigidity. It is therefore 
assumed that the connected tetrahedra will have 
minimum density, and this can obtain only if the 
bonds form angles Si-O-Si of 180° and the rings, 
if any, are circular and planar so that the tetrahedra 
can have maximum distention. (See fig. 3.) Under 
these circumstances, probability greatly favors the 
formation of pentagonal rings, mentioned by 
Poncelet [13], because the interior angle of a regular 
tetrahedron is 109° 28’ and differs so little from 108°, 
the interior angle of a regular pentagon. 

Growth of this structure of minimum density pro- 
ceeds uniquely in 6 directions from an initial tetra- 


as 





hedron to form a network containing dodecahedral 
holes or cavities, each separated from its 12 some- 
what-distorted neighboring cavities by 12 planar 


interfaces of very slightly stressed penta rings. An 
important property of this model is its many axes, 
15 of twofold, 10 of threefold, and 6 of fivefold 
symmetry. Each of the 15 twofold axes makes an 
angle of 36° with each of its 4 neighbor axes, and 
similarly the 10 axes are equidistributed at angles of 
11° 48’ and the 6 axes at 63° 26’. These features 
mean close relation to crystal symmetry and a high 
probability of approximate alinement of different 
portions of the structure even if separated by less 
well alined interstitial regions or tissues. 

The most important property of this proposed 
model is its fivefold symmetry, which precludes 
formation of crystals. Just as penta tiles fail where 
hexagons succeed in covering a floor, so do regular 
penta dodecahedra fail to fill space completely and 
extensively. The interior dihedral angle of a 
dodecahedron is, however, 116.6°, and thus not far 
from 120°, which would be exactly right for filling 
space without stress and distortions. If one regular 
dodecahedron is attached or joined on each face to 
12 like structures, none in the outer shell can touch 
any of its 5 neighbors. By radial compression these 
neighbors may meet and the requisite tangential 
tension in the peripheral bonds may not be excessive 
for the first shell. Further increase in cluster size 
by the addition of 32 cages in the second shell would 
require much higher tensions in the outermost bonds, 
and some defective attachments or broken bonds 
might result in order for bonds to hold elsewhere on 
the periphery. Thus the great strength of the 
silicon-oxvgen bond is the factor that should cause, 
end vet automatically limit, the growth in size of 





FIGURE 3. 


View of 


regular pentagonal dodecahedral cage formed of planar 


pentagonal rings of tetrahedra as shown at lower center. 


The high-temperature low-density planar ring proposed for silica glass contrasts with the nonplanar 


hexagonal elemental rings in the denser silicas tridymite (left 
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and cristobalite (right 





vitrons. The automatically limited ability to grow, 
and a distribution of stresses localized and balanced 
within each vitron, are features that may distinguish 
a vitreous substance from a crystal. 

Silicon ions are situated at each of the 20 vertices 
of a pentagonal dodecahedron, and oxygen ions at 
the midpoint of each of the 30 edges. An isolated 
unit of this sort would consist of 20 silicon atoms and 
30 of oxygen with 20 plus charges, and thus be repre- 
sented as (SinQg)"*. <A cluster of 7 such units is 
(Sigs0,45)* but it is decidedly nonspherical. With 
13 such cages there is again an approximation to 
sphericity with a diameter of 20 A (see fig. 7). This 
13-cage cluster can be represented as (Sij30Q29)*, and 
the plus charges are distributed over the whole 
surface in 12 pentagonal groups. Upon addition of 
the next shell, consisting of 32 dodecahedral cages, 
to form a 45-cage cluster the diameter increases to 
30 A and the formula is (Sigg6Og)*. The plus 
charges indicate the degree of attachability of these 
vitrons to their neighbors through oxygen bonds. 
In addition, or alternately, they may actually unite 
or grow together with neighbors to limited ex‘ents, as 
further discussed in the appendix. 

Figures 4, 5, and 6 (also 9 and 10) are helpful for 
visualizing the arrangement of atoms in the proposed 
model. With dimensions based on the average 
X-ray datum of 1.60 A for the Si-O distance (be- 
tween centers), an element or cage can be inscribed 
in a sphere of 9 A diameter, and a sphere of 7 A 
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Figure 4 
element of a 


Schematic projection of pentagonal dodecahedra 
plane 


vitron on perpendicular to two of its 


pe nlaqgonal faces. 


Fieure 5 Schematic dodecahedra 
of its siz ares of fivefold 


ere rie ni riewed atong Orne 


symmetry through pairs of} opposite 


and paralle l faces 
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diameter can be inscribed inside a cage. Only a 
central cage of a vitron can be truly regular and 
symmetrical. The surrounding cages must be de. 
formed in order to contact their neighbors. Hence 
it is interesting to recall that Warren [14] estimated 
8S A as the maximum dimension of the region of 
periodicity in fused quartz, and that Valenkov and 
Poray-Koshitz [15] estimated 12 A as the exten 
of perfect order, with surrounding regions of inter. 
mediate degree of order. Also, measurements on 
thermal conduction in glasses have made it possible 
as shown by Kittel [50], to estimate 8 A as the mean 
free path distance over which plane elastic heat 
waves (phonons) may travel in silica glass at room 
temperatures before distortional scattering. 


Dodecahedral element viewed 
ares of five fold symmetry 


FIGURE 6 aiong one of its siz 


“ome outer atoms are 


omitted 





I Ist RE ¢ View along one of 15 identica twofold axes ot 
symmetry of a 13-element vitron fy pe ]. spherica 
In this model a tetrahedror represented by two wires ben 


49.5° and soldered at their symmetrical crossing to represent 
sented by friction-tight couplings Iwo other axes of tw 
normally to the line of ht 
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. hes = 3s 
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; These densities are only 10 percent lower than 2.20, | = 7 2 § 
the density of silica glass from room temperatures to | = : 2 ; 
1.400° C or higher. . a2 SYERS AS 3 
Possibly such low values of the densitv of fused | =) 2) EZ 3 ee < 
silica may be approximated in its liquid state at very | = SE | 3 F 
: high temperatures at or near vitron formation. How- | 3 "ia z 
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when silica is cooled, because of the noncrystal nature & 3 £ see 28 gE° 2 
of its structure with fivefold symmetry and the dis- | S| = 22 7. . 
tortions and fissuring inherent in the process of its p == F 2 
macro-growth. The tensions induced in the outer | & - Pe 8 
8 siz bonds, as the vitrons grow in volume, probably lead | = . e - 3 it 
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Aside from densification through such a decrease | 2) - 42 > 3 _ — oe > 
in the unit volumes of the outer and distorted cages | = ts 1 a 
of vitrons, there would be a variation in the effective = si < 2 
density of vitrons according to their size. The in- | *% - " aie ‘ 2 |e 
ternal portion of each elemental volume or cage con- | = Ee SASEZS2n =Baee Se gan 
, tains 5(SiO,) in 251 A®, and the volume of the atoms | =| = | Z& | 2 —— ee 2 
is only about 95 A®. Outside an isolated small vitron | = mS | 3 4 
there may be almost as many atoms as inside. The | = 7 t E 
outside atoms must project into interstitial or chan- | > - a me m * 2 | § 
nel tissues and add proportionately more to the mass | = 23 SZRass- S¢aas =e ze% | = 
of the vitron than they do to its effective volume. A | ‘=o |) 22) 3 iis ee & 
useful and attachable elemental unit volume to con- | & =s | 8 3 
sider is (SijQ) with 10+ and 10— associated | & ul fi. + 
charges. Such a unit with balanced charges is at- | & " ai cad sd wel 2 i 
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cluster vitron (SiggQy)*" The unit (SijoQOo) "| 4 || = | EF) 3 oo te a. 3 
has 5(SiQ,) inside and the same number outside its | 2 os | = 5 
251 A® limits. Considering this unit as a single- | & ae A = 
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a minimum thickness of 1.6 A (the Si-O distance). | Sis eB ig & FE 
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as computed on the continuum basis. In other 
words, the interstitial tissues can be denser than the 
internal portions of vitrons, even if the bonds therein 
be fewer, less regular, or unconnected, because these 
tissues cannot share in the extreme degree of rare- 
faction in space that characterizes silica in dodec- 
ahedral structure. 

The results of a study of vitrons of various sizes 
are given in table 1. It is considered that at the high 
temperatures of formation the vitrons grow simul- 
taneously in situ throughout a melt and that smaller 
vitrons will grow in every sizable region that might 
otherwise remain undeveloped between larger vit- 
rons. All these vitrons of various sizes and perhaps 
of somewhat irregular shapes can be actually, if im- 
perfectly, united with neighbors at some common in- 
terfaces, and by multiple oxygen bridges across inter- 
vening channels can be more or less securely attached 
to each other. 

It appears from table 1 that the maximum effective 
density exceeds the constant internal density of 2.0 
by various amounts depending on the size of the 
vitron and the degree in which its shape approximates 
the sphericity assumed in the estimates. It will be 
evident that nonspherical clusters must have greate! 
peripheral areas than assumed, and hence their esti- 
mated effective densities must be too high. However, 
there are several clusters for which the estimated 
maximum density is less than 10 percent above 2.20, 
the known density of silica glass. Further details 
concerning vitrons and their aggregation are given 
in the appendix. 


4. Diffusion of Gases in Silica Glass 


An important feature of a model for the network 
of glasses with high silica content is the minimum or 
limiting aperture through which atoms must pass 
when gases diffuse through glasses and when ions 
migrate under the influence of an electrical field. 
Experiments concerning the diffusion of gases 
through glass may not, in general, be considered 
solely on a mechanical basis, as a process of the 
squeezing of molecules through molecular holes or 
cracks. The more reactive hydrogen diffuses slower 
than the inert helium when oxygen is present, but 
their rates are about the same in a solid without 
oxygen anions. 

However, in the case of diffusion of noble gases 
it is probably valid to consider dimensions as much 
more important than chemical nature, and in seek- 
ing to understand the very marked diminution in 
the permeability of silica glass to argon, in compari- 
son with that for neon, it seems tenable and sufficient 
to regard size only. As measured by Norton [12], 
the diffusion rate for neon with atom diameter 2.4 A, 
was found to be more than 100,000 times as high as 
the rate for argon, atom diameter 3.2 A. Figure 8 
shows how well the penta faces of the elements of 
vitrons in silica are adapted to pass neon and smaller 
atoms such as deuterium, hydrogen, and helium, but 
‘ obstruct argon, oxygen, and nitrogen, which have 
larger diameters. 

Norton found, also, that the rates for helium at 
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i 
Frat RE &. Limiting aperture for diffusion of gases through 
silica glass. 
Che clear iperture of the pentagon il faces of the vitron elements is comparable } 
n diameter with the neon atom which is the largest that rapidly diffuses in 
quantity 


100° C are decreased by a factor of 1,000,000 as the 
percentage of modifiers is increased from glass to 
glass. It is here suggested that this difference is 
caused mainly by the physical presence of nonsilica | 
packed inside the elemental cages of the vitrons. 
It does not seem necessary to look to the interstices 
between vitrons as the only channels for permeability 


in silica glass. | 


5. Radial Distribution of Atoms 


All atoms in this model are found to be distributed 
in accord with published data on the radial distribu- 
tion of atoms with respect to any one ion in vitreous 
silica. That is, for very short distances to their 
nearest neighbors both Si and O atoms are so situated | 
that the arrangement is not distinguishable from that 
formerly assigned to high cristobalite. Warren [16] | 
was a leader in applying the radial-distribution 
method to the analysis of X-ray diffraction patterns 
obtained from glasses. The concept of a random 
network required such treatment instead of the | 
simpler computation of Bragg-d values as is done for 
crystals. Hartlief [17] also reported such distances ; 
for fused silica from X-ray data. More’ recently, 
neutron diffraction data for fused silica were treated 
in this manner by Milligan, Levy, and Peterson [9}; 
and Weber [18] et al. have reported similar unpub- 
lished data on two samples of silica glass, one of 
“high purity”? made by the Corning Glass Company 
and another of “ordinary” quality. These expert 
mental results are listed in table 2 for comparison 
with the atomic separations as computed from the 
model made by the writer. (The seale value of 
(Si-O) = 1.60 is, of course, assumed in the computa- 
tions. ) 
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TABLE 2. 


Radial distribution of atoms in silica glass 
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Distances to neighbor atoms in angstroms 


From O to From Si to 
Proposed | \ruttigan 
model lowe & 
y < 
Peterson [9] 
No. Kind No Kind 1951 
2 s { 0 Loo 1. SS 
ie) 200 2. oH 
t Ss 3. 20 3. 02 
t = l 2) Loo 4.0 
() + 1u 
~ © 4.92 | 
12 0 uy 1 
2 s 5.18 
12 0 US | 
u 
12 ~ { () 6H. 14 | 
o 2) 6.79 hit 
2 0 7.14 } 
t = 9 © 7. 15 
12 oO | 7 4 
1 = 7.33 
2 s 4 0 oe 
2 ) s.04 
‘) 2 ~ &. OA 
24 S s.4 
~ () 8. O4 
Is ~ 8.6 
{ ~ & U7 if 
Read in 1956 from ‘‘uncorrected”’ curves in Prof. Wa 
Norway, just as tl paper goes to press, lists his “‘une 
1 5A 
Not i 
No attempt is made to compute the angle 
‘* ~ . “7 ° rr 
Si—O-—Si from the diffraction data. The departure, 


a. of this angle from 180° might be estimated by 


means of the equation 

COS (a 2 Si-S1 , Si-O). 
but the precision is much too low, and even for 
averages there must be residual error because the 
peak corresponding to the O-O distance of 2.6 A is 
usually imperfectly and that makes the 
reading for Si-Si systematically too low and the com- 
puted departures from 180° too large. 

The proposed model agrees with the experimental 
data of table 2 somewhat better than does a random 
network model. A random orientation around bond 
axes soon leads to indefinite predictions of the radial 
The first uncertain case is the oxygen to 
second-oxygen distance. For random orientation 
there are 18 oxygens that might be found anywhere 
within the range 4.2 to 5.2 A. The last fixed dis- 
tance would be silicon to 12 second-silicons at 5.2 A. 
Thereafter 18 third silicons might vary from 4.9 to 
6.5 A from an oxygen, and 36 third oxygens be 
found anywhere within the same limits of radii from 
any silicon Incidentally the uncertainties start 
with nearer neighbors and rapidly become much 
greater if values other than 180° be assumed for the 
angle Si--O—-Si Weber [19] was evidently sur- 
prised to find only 6.9 instead of 18 oxygen atoms at 
4.46 A from an oxygen, and he remarks, “. .. A 
random orientation of these outlying tetrahedra 
probably causes 1 oxygen of each 3 to be closer in 


resolved, 


distances. 


Neutrons 


X-rays 


Weber [19] 1954 


Warren [15] Hartlief [17] Nerman * 


et al. 1936 1938 
Corning |‘ Ordinary” 
1 wa 1 1. 62 1.7 1. 64 
2. 62 2.71 b (2. 65 b (2.7) 2. 67 
3. 32 3. 22 3.2 3. 1 3.14 
3.04 4.2 4.1 4.17 
4.46 

5 ON 5.2 5.3 5. 25 | 

' 

| 

| 

6.3 6.5 | 

| 

| 

| 

en's laboratory 4 letter received from N. Norman of Oslo, 


wrected”’ and still unpublished values as 1.62, 2.62, 3.05, 4.08, 


vielding 6 as the number of oxygens at the O-2d O 
distance.”’ In table 2 it will be seen that the pro- 
posed model predicts that 6 of these second oxygens 
should be found at 4.19 A and 12 of them at 4.99 A, 
very close to 3 other distances, so that only the high 
radial density near 5.1 to 5.3 A is noticed by most 
observers. 


6. Planar Distribution of Atoms and Values 
of Bragg-d 


Computations of values of Bragg-d are not usually 
made for glasses if it is assumed that their structure 
is random. When a considerable degree of order 1s 
assumed, it then becomes pertinent to consider such 
tests. In contrast with the simplicity of crystal 
spacings, those in this model are very numerous and 
multivalued. Within a single dodecahedron, which 
seems useful as a unit volume, one finds the atoms 
all distributable in 6 equivalent sets of 7 planes, each 
set perpendicular to one of the 6 axes of fivefold 
symmetry. Also, the same atoms are distributable 
in 10 sets of 11 planes, each set perpendicular to 
one of the 10 axes of threefold svinmetry ; and in 15 
sets of 13 planes, each set perpendicular to one of 
the 15 axes of twofold symmetry. Figures 9 and 10 
show how the atoms may be projected to locate 
these three different systems of planes. 

Among the very numerous spacings of figures 9 
and 10 there are many duplications, and one finds 
only 31 different possible values of plane separations 
for a single cage between limits of 0.8 and 4.5 A and 
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The aton ire further prejccted to detern the spacing of 1 planes perpe ular to tl In twotold axe 


considering all three sets of projections. This can 
be compared with about 35 peaks (or indications 
thereof) that can be used for computing d-values 
from the curves of scattering of neutrons as found 
by Weber [18] for a sam ple that he labels as Corning 
fused silica; or 18 peaks for a sample called ordinary 
fused silica. 

Because of the numerous in intervening 
planes, no X-ray values of Bragg-d should be ex- 
pected for plane separations exceeding about 3 A, 
as judged from analysis of the data in figures 9 and 
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atoms 


10. Below 3 A it ts estimated by a rough process of 
atom counting and weighting that about 10 d-values 
might be detected by X-ray and electron diffraction, 
and 6 or 7 are found from such data. 

For vitron clusters, wholes, there 
other possibilities which seem chiefly to be distorted 
values of the separations as shown in figures 9 and 10. 
An example is the d-value near 4.5 A that 1s not 
predicted from the model on the basis of a single 
element but is prominent in X-ray data. In clusters | 
there takes place an overlapping of the vitron ele- 
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Figure 10. 


mms are further projected to detert 


ments. An important case seems to be the over- 
lapping of fivefold and twofold axes, and this was 
first noticed in a model of the 13-element vitron as 
viewed transversely to any of its fivefold axes. The 
distortions are such that approximate coincidence is 
forced for the 4.19 A separation along the twofold 
axes and the 4.40 A separation along the fivefold 
axes. The hybrid “planes” are curved with separa- 
tion varying from 4.40 A along the principal five- 
fold axes of the whole vitron to 4.19 A at distances 
of something like 6 A from the principal axis (see 
figure 11). 

In brief, because of the limit of precision in the 
experimental data one can say only that the number 
and values of Bragg-d separations predictable from 
the proposed model agree with the diffraction data 
about as well as the agreement among different ex- 
perimenters or for different samples of silica glass. 
It is suggested, however, that the large number of 
these distorted pseudo planes has an important effect 
on the level of observed scatter and the general 
fuzziness of diffraction patterns. 
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Projection of dodecahedral element on plane perpendicular to one of its 15 axes of twofold symmetry. 


sets of pseudo lattice spacings of 7 planes perpendicular to the fivefold axes and also spacings of 11 planes perpendicular 


7. Potential Flaws and Tensile Strength 


The properties that may logically be inferred from 
the nature of vitrons, as they have been here defined 
and described, seem to offer interesting explanations 
of some of the unusual and not satisfactorily under- 
stood properties of silica. For example, in the testing 
of bulk specimens the actual tensile strength of fused 
silica (and other silica glasses) is only a fraction of the 
computed “cohesive” strength. The Griffith [20] 
crack theory is the usual “explanation,” namely, that 
glass is completely riddled with flaws, but there seems 
to be no acceptable basis for the existence of the cracks. 
Stanworth [21] says that it seems as if the entire 
volume of glass may be permeated with cracks. 
Cox [22] considers some of the facts as suggestive 
of crack propagation in an elastic continuum after 
initiation by a flaw, but on the other hand finds 
that associated factors lead to an atomic view in- 
volving period and energy of vibration and size of 
ions in order to explain what may be “‘macroscopic 
manifestations of the discrete nature of glass.” 











View perpendicular to one of the 6 fivefold axes 
of symmetry of a 13-element vitron. 


FIGURE 


Strings show locations of prominent curved or hybrid “planes” in which nu me: 
ous ions are located. The separation of adjacent hybrid planes of like curvature 
is 4.4 A along the fivefold axes and 4.2 A along approximately parallel (secondary 
twofold axes at distances of 6 A from such principal axes. A similar forced 
approximation to coincidence is found at right and left of center along the 15 axes 
of twofold symmetry, to which secondary fivefold axes are approximately parallel. 
These separations of 4.2 and 4.4 A are of interest because the most prominent 
diffuse ring commonly observed in both X-ray and electron diffraction patterns 
indicates a Bragg-d of 4.2 +1 A. 


According to the vitron concept, the intrinsic 
tensions of distortional origin cause weak peripheral 
bonds that provide a basis for incipient cracks, and 
there are fewer bonds per unit cross section in the 
interstitial connective tissues between vitrons than 
within them. This explains the low glass strength as 
compared to a medium having bonds of uniform 
strength uniformly distributed. 

Murgatroyd [23] discussed the delayed elastic 
effects in glass and attributed the effects to gradual 
yielding of weak bonds. He proposed that the con- 
stitution of silicate glasses is such that the binding 
force between atoms varies from the maximum 
possible value to a relatively small value, and ‘re- 
gions occur where numerous bonds of the weakest 
type exist as neighbors.”” He concluded that “the 
low strength of massive glass is due to the numerous 
weak bonds it contains.”’ 

Murgatroyd also mentions the smallness of the 
difference in density between liquid silica glass and 
cristobalite, and points out that silica resembles a 
normal liquid in this respect. He observes that one 
way in which the great increase in fluidity may be 
accounted for is “the presence of regions of weak 
bonds interspersed amongst normally bonded groups 
of atoms.” 

On the network 


basis of a continuous random 


theory there is no satisfactory evidence that weak 
bonds would be sufficiently numerous to account for 
the very great discrepancies between theoretical 
cohesive strengths and maximum realized tensile 
strengths of glasses, and especially so for fused silieg 
glass; but the vitron concept, with very numeroys 
bonds that are seriously tensed and interstitial cop. 
nective tissues where the bond density is low, seems 
definitely satisfactory in this respect. 

Much of the testing of glass for breaking tensile 
strength has been done on fibers. The values 
obtained are higher than for tests on samples with 
large cross sections, and especially so if the fibers are 





a 


- 


drawn in approved manner and great care taken to | 


avoid damage to their surfaces. If ambient condi- 
tions are favorable during drawing and if average 
results are plotted against diameter of cross section, 
or formation temperatures as was done by Otto [40], 
the curves can be extrapolated toward smaller 
diameter or higher temperature for somewhat higher 
estimates of the breaking strengths, which however, 
are still much below expectations from the computed 
cohesive strengths of 2.5 to 510° Ib/in’ for com- 
mercial glasses [39]. 

Considerable attention has been given to reasons 
for higher breaking strengths of fibers as compared 
with massive glass. Chilling should introduce com- 
pressions in the surface layers and increase breaking 
strengths to some extent but this effect must become 
less important as diameters are decreased because 
there must be interior compensating tensions.  Ori- 
entation of structure into threadlike alinements, 
caused by the drawing process, has been much dis- 


——- 


— — 


~——— 


~ 


~ 


cussed, but Otto and Preston [48] found little differ- 


ence in strengths whether by tensile or torsional tests; 
more recently Otto [40] and Brannan [49] have 
reported high values of strength on fibers that do not 


show anisotropy; and Goldstein and Davies [47] did | 


not find anisotropy in silicate glass fibers although it 
was present in some others. Additional pertinent 
evidence has been presented by Lynch and Tooley 
[46], who find that their observed increases in 
strength (up to 62%) by down-drawing under in- 
creased load are not eliminated, although materially 
reduced, after abrasion of the surfaces of the fibers. 

The idea of change in either statistical or spatial 
distribution of types of vitron is not analogous to 
elongational distortion of a given type of vitron. 
The stresses and distortions involved may, however, 
be smaller and thus union with neighbors may be 
improved. In any event the stress systems would 
necessarily be balanced on a vitron-volume level, 
and it would seem that vitrons at and near the sur- 
face could more easily meet the requirements for 
strong unions and attachments with neighbors 
because of less exacting three-dimensional require- 
ments. This would mean a maximum of strong 
intervitron connections in all directions near the 
surface, and a progressively smaller proportion of 
strong connections inside until the distance from the 
surface becomes large compared to the vitron size. 
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8. Effects of Localized Stresses on Extensive 
Properties and Elastic Moduli Above Room 
Temperatures 


The unusual extensive and elastic properties of 
silica glass have been so well discussed by Babcock, 
Barber, and Fajans [6] that some details and many 
references will be omitted here. In order to e xplain 
these peculiar relationships those authors suggest the 
coexistence of two or more different ionic structures 
in equilibrium in ratios that vary with temperature 
and pressure. Also, they consider that the low- 
temperature heat capacity, as well as infrared and 
Raman spectra, may ‘indicate the coexistence of 
relatively weak and strong forces.”’ 

An experimental curve showing the temperature 
variations in the modulus of rigidity of fused silica 
indicates strong thermal variations in stresses within 


the components of its structure. This modulus 
decreases from 2° K to a minimum near 65° K, 
rises to & Maximum near 880° C, and then has an 


extremely rapid decrease near 1,020° C [24]. It will 
be noticed that the minimum near 65° K coincides 
with the minimum for Young’s modulus of silica 
glass, the maximum at 880° C is somewhat higher 
than the similar maximum for Young’s modulus 
(not far from a broad weak minimum in expansion 
coefficient), and at 1,020 the lower limit of the 
annealing range has been reached where viscosity 
approximates 10" and localized differences 
in stress are not maintained indefinitely. 

It is here suggested that the localized stress sys- 
tems consist of (1) intratetrahedral attractive and 
repulsive forces that promote and maintain angular 
symmetry within the tetrahedra, and (2) the tensed 
Si—O—Si bonds that make and maintain contact 
between tetrahedra through the oxygen bridges 
contacts that cannot exist without internal distortion 
in the tetrahedra. 

As hot glass cools toward the critical lower limits 
of the annealing range, the strain response to stress 
has become decidedly time dependent. For silica 
glass in particular (free from the complication of 
modifiers) the repulsive and attractive forces within 
tetrahedra increase and maintain greater symmetry 
therein at the expense of greater tensions in the Si 
O—Si bonds between tetrahedra. This necessitates 
more regular and symmetrical dodecahedra which 
can obtain only if some connections are broken and 
interstitial tissues or fissures are extended or widened. 
This is confirmed by the marked decrease in activa- 
tion energy of flow (slope of login versus 1/7) for 
silica glass at te mperatures below its annealing range 
where viscosity increases very slowly above logion 13. 

After the initial rapid increase In cohesive forces 
during cooling near 1,020° C has produced increases 
in the moduli there are slower progressive changes 
or relaxations in the balanced internal stress systems 
of vitrons that produce marked but gradual decreases 


poises 


in the moduli but have very little effect on volume. 
In respect to volume the normal contractions of the 
units of structure are nearly offset by the expanded 
channels between vitrons over a long range in 
temperature. Although these contractions permit 
some readjustments and relaxations of stresses in the 
tetrahedra that in turn check the initially excessive 
fissuring effect, there is still progress in the fissuring 
(because of increased attraction and consequent 
increased symmetry within the tetrahedra) so that 
interstitial (intrachannel) cohesion can decrease 
slowly and the concurrent increase in intravitron 
cohesion be negated by the mere existence of the 
weaker channel interstices. 

From the above discussion it is evident that the 
net coefficients of thermal expansion of silica glass 
can be very small and nearly constant between 1,000° 
C and room temperatures. Furthermore it is clear 
that these net coefficients at increased pressure should 
be larger because the negative components caused 
by stress-induced fissuring would be correspondingly 
suppressed by the increase in pressure. 

The analysis given also suggests that the compres- 
sibility of silica glass should be abnormal and in- 
crease with pressure because the tensed bonds in and 
near the interstitial or channel tissues should aid the 
applied pressure and do so more effectively as and 
until the atoms approach their equilibrium positions. 
This has indeed been observed up to pressures near 
36,000 kg/em? where Bridgman [25] found that the 
abnormal increases in ¢ -ompressibility ceased abruptly 
and were then followed by decreases, in a normal 
manner, at higher pressures. 

According to the automatically-limited-growth 
feature of the vitron theory, and to the extent that 
Bridgman’s pressures are equivalent to hydrostatic 
pressures, his 36,000 kg/cm? applied radially to each 
vitron would seem to substitute directly for a radial 
compression previously produced or maintained in- 
directly by tangential tensions of (36,000/c) kg/cm’, 
where ¢ is Poisson’s ratio which is near 0.16 for silica 
glass at room temperatures. Such tangential ten- 
sions of 225,000 kg/cm? or 3.2 10° lb/in.’ are in very 
reasonable agreement with estimates of theoretical 
cohesive strength. Shand [39], for example, gives 
2.5 to 510° lb/in.2 for common commercial glasses, 
and thus this abnormal compressibility observed by 


Bridgman seems quantitatively to confirm’ the 
writer’s inference (section 2, above) that vitron 
growth is limited by the strength of the Si—O—Si 


bonds. 

Moreover it will be realized that the compression 
produced by a given pressure should be lower at 
increased temperatures because the process of fissure 
closing by applied pressures is then applicable only 
to the lowered level of fissuring that can exist at the 
higher temperatures. This accords with observa- 
tions and measurements by Birch and Dow [26] 
from room temperatures to 400° C on fused silica 
and a Pyrex glass. 
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9. Internal Friction and Elastic Constants 
at Low Temperatures 


A study of thermal effects on elastic moduli at low 
temperatures invites consideration of the internal 
friction evidenced by the attenuation of ultrasonic 
waves observable from 2° to 190° K and having very 
broad maxima near 30 to 50° K depending on the 
frequencies used (60 ke to 20 Me/sec) 

Such friction effects are very small in quartz and 
observable only with care. In fused silica they were 
observed by Marx and Sivertsen [38], by MeSkimin 
[36], and by Fine et al. [37] before their more detailed 
observation and discussion by Anderson and Bommel 
[24]. Fine’s data show minima in Young’s modulus, 
torsion modulus, and Poisson’s ratio near 65° K 
and a very broad maximum in compressibility near 
65° K. The Anderson and Bémmel data show a 
minimum in rigidity at 65° K for 66 ke and near 

° K for the 20-Me curve. The last named authors 
find that the temperature and frequency behavior 
suggest a structural relaxation with a distribution of 
activation energy averaging 1,030 cal/mole. Further 
they find that the absorption is largely due to a 
shear elastic wave and therefore to a mechanism 
that does not primarily tend to cause fluctuations 
in density; and that the relaxation strength is small, 
only a fraction of all particles contributing, as shown 
by the small increase (about 3%) in shear modulus 
from the minimum value to that at 4° K. 

For discussing these results in relation to vitron 


theory it is convenient to start near room tempera- 
tures where there is a noticeable decrease in co- 
efficient of thermal expansion. This may mean 


merely that the rate of expansive fissuring has fallen 
definitely below that of normal contraction and con- 
tinues so until minimum volume is reached near 
190° K. Here the above described internal friction 
begins very gradually to be detectable by attenuation 
of acoustic waves, and the compressibility after 
remaining relatively constant over an interval of 
+75° near room temperatures resumes an abnormal 
climb as temperatures are lowered. It is here sug- 
gested that because of continuous and_ gradually 
increasing intravitron contraction some of the bonds 
yield in and near the interstitial channels between 
vitrons. This would provide limited numbers of 
less well bound oxygen for the sort of transverse 


vibrations visualized by Anderson and Bémmel, and 
do so in ever increasing number until the whole 


process is checked by the near contacts and strong 
repulsive forces at very low temperatures. Appar- 
ently lower temperatures are necessary to check the 
absorbing vibrations than those required to start 
closure of the fissures and decrease compressibility. 
At the lowest temperatures of observation it will be 
noticed that the rigidity increases and compressi- 
bility decreases as is normal on approaching absolute 
zero. The greater slopes on the low-temperature 
side of all of these curves may indicate that a normal 
‘balance between repulsive and attractive forces is 
effectively reached at the temperatures of their 
respective Maxima or minima near 65° K. 
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10. Vitron Theory and Heat Treatments 


According to the theory of intrinsically stressed 
vitrons, after fused quartz is repeatedly heated and 
cooled in treatments that reach 900° or 1,000°C 
these repeated heat treatments should slightly reduee 
the localized stresses present within the vitrons and 
result in small volume changes at a given room tem. 
perature as compared with previous measurements 
at the same temperature. Such changes should be 
unidirectional (expansions) and permanent, that js 
not reversible by subsequent heat treatments jp 
comparable temperature ranges. Considerable eyj- 
dence of this has been found in the course of anneal- 
ings and in measurements of the thermal expansion 
[42] of vitreous silica. Possibly some effects of this 
sort should be considered in connection with the 
results published by Saunders and Tool [43] on 
Pyrex glasses. 

Such irreversible heat effects, or effects reversible 
only after more fundamental heat treatments at 
much higher temperatures, should be less and _ less 
evident as the nonsilica content of glasses is increased. 
and they were not of noticeable magnitude in the 
results reported by Tool [44] and associates on glasses 
other than Pyrex. 

At very high temperatures in liquid silica the 
stresses incident to vitron growth produce strains 
and distortions easily, and there is no rigidity of 
tetrahedral structure sufficient to resist the tendeney 
for dodecahedra to unite and fill space with a min- 
mum of interstices. Upon cooling and approach to 
the solidification region, even slight trends toward 
greater symmetry within tetrahedra lessen the plia- 
bilitv of dodecahedra, create a few broken bonds, 
and in general form or widen interstitial tissues or 
fissures. This expansive effect may exceed the 
normal contractive effect and vield a small net nega 
tive coefficient of expansion such as was observed by 
Douglas and Isard [27] between 1,000° and 1,500° C, 
Small spheres of vitreous silica had previously been 
quenched by Salmang and von Stoesser [28] from 
high temperatures into water and they had announced 
the existence of a minimum density for this glass 
between 1,.200° and 1.800° C 

These positive changes density of silica glass 
with increase of temperature in and above its anneal- 
ing or “transformation range’’ are very small, say 

1<10~°/deg C, and contrast with negative values 
of from 16 to —46 10° /deg C that were found 
by Tool [44} for various silicate glasses in their an- 
nealing ranges. These larger and negative changes 
are attributable, according to the writer, almost 
entirely to the silica network but the cause may be 
very different. It is suggested by the writer [45] 
that there is an inherent temperature dependence of 
the silicon-oxygen-silicon bond angle that is locked 
in @ semi-invariant or metastable state in pure fused 
silica, but is free to operate in silicate glasses in pro- 
portion as they contain modifiers. The result is @ 
limited, temperature dependent, and reversible de- 
gree of folding or collapse of the network structure. 
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11. Specific Volumes of Soda Silicate 
Glasses 


A useful feature of vitron theory is the degree in 
which it may be quantitatively applied in studies of 
the properties of simple binary glasses as functions 
of composition. Applications of the vitron concept 
to specific volume, viscosity, volatilization, and elec- 
trical resistivity of glasses are being made and data 
on infrared spectra are being studied. It is found 
that the proposed ideas accord with many such data 
for fused silica and for simple alkali-silicate glasses. 
These results will be given in another paper with 
emphasis on property-composition curves. 

However, one such application to binary glasses 
ean be readily summarized here. Since the number 
of cages in R,O silicate glasses is proportional to S/5 
(where constituent proportions R+S=1), it follows 
that the proportions of RO corresponding to integral 
numbers of cations per cage, 2), are readily computed 
as R=R,/(R,+10). At such integral steps or de- 
erees of saturation (and supersaturation) of the 
cages it seems probable that some evidences of 
changes in rate might be found on curves of property 
versus composition. Huggins [35] has already stud- 
ied existing data of volume-composition curves for 
soda-silicate glasses as published by Glaze [41] et al., 
and found evidence of special points at Si to O ratios 
of 0.437, 0.40, 0.38, 0.36, 0.345, and 0.333 which are 
now found to correspond precisely to values of 3, 5, 
6, 8, 9, and 10 for R, according to vitron theory. 
This good agreement between vitron theory and the 
results by Huggins from his investigation is shown 
in further detail in table 3. 


TABLE 3. I pplication of vitron theory to spec ific volumes of 
soda-silicate Glasses 
Special ratios of silicon to oxygen 
Mola 
compo Hol Cat 
Ss Computed Found by Huggins 
Py irom vitron 
nodel 
Nao SiO 1940 1954 
ooo iH) 0. 20 { it 4) a 4) 
vol », G09 x2 476 
7 S34 67 2 $54 
2 769 ‘ 43° 44 44 
2M ; } } 417 
333 7 ; 10) 1) 0) 
375 62 2 “t 7 
1] sv 8 7 
444 t ‘ 7 s4 $a7 
mn { ; { 44 
iL f " {3 333 
+) { oT 
R be vel ‘ } by 
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Study of this model for silica glass emphasizes 
two pertinent characteristics. First, it has a type of 
order with many axes and numerous and varied 
spacings of quasi lattice planes that cannot grow 
into a sizable crystal, and second, it has distortions 


Concluding Remarks 
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and stresses and consequent curvature of the lattice 
planes. Because it has been shown that this model 
is consistent with known data on density, differential 
diffusion of gases, interatomic distances, low tensile 
strength, unusual thermal and pressure effects on 
volume, internal friction at low temperatures, and 
specific volumes of soda-silicate glasses, a new con- 
cept should be considered for general study of the 
glassy state: A small stressed cluster of definitely 
arranged but slightly distorted elemental groups or 
cages (such as pentagonal dodecahedra) that in 
turn are composed of smaller and more nearly regular 
stable structural units (such as tetrahedra), all of 
which as a system tends to form at high temperatures 
in a melt and have an approximate but not exact 
noncrystal symmetry (such as fivefold) and thus be 
capable of limited but not extensive continuous 
growth. Such a cluster may appropriately be called 
& VITRON in order to distinguish it definitely from a 
crystallite, which can increase in size without limit. 
Among the member ions and in connective tissues 
between the vitrons there will exist a wide distribu- 
tion of bond strengths and some variation in angles 
between valence directions that can account not 
only for a range of melting points but for similar 
phenomena indicative of ranges in activation ener- 
gies. The known high energy of dissociation of 
SiO, linkages (106 keal/mole, an energy which is 
many times the mean energy of agitation at most 
industrial temperatures), indicates that vitrons are 
not necessarily completely broken up into tetrahedra 
even at temperatures far above fusion. 

It would seem that the vitron concept should be 
welcomed by those who still regard a modified erys- 
tallite theory with favor. It is not difficult to vis- 
ualize these very small submicroscopic vitrons (say 1 
to 54-cages, approximately 9 to 30 A in total diam- 
eter) partially united to each other through common 
faces and also connected by other Si—O—Si bonds 
or “oxygen bridges,”’ similar to those that unite the 
tetrahedra within a vitron, and thus forming a 
structure approximately filling space with intervening 
connective tissues that vary in width with tempera- 
ture and degree of existing localized stresses. Con- 
sidering their various multiple axes, the vitrons 
could, almost by chance, be positioned to present 
more or less equivalent orientations regardless of 
intervening tissues or fissures. This simple but 
fundamental concept of connected vitrons intro- 
duces more than short-range ordered arrangement 
and thus something that proponents of the crystal- 
lite theory have demanded. A quantitative inter- 
distribution of stresses between elements of a vitron, 
and perhaps between vitrons, helps to explain the 
remarkable continuity of properties as a function of 
temperature through, during, and even above the 
melting range—something that, according to Peychés 
[29], a mere distribution of sizes and shapes does 
not satisfactorily explain. 

Also, the intrinsic properties of a vitron, including 
numerous lattice planes, a distribution of stresses 
among its members, and curvature of the lattice 
planes, can account for much of the observed 
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broadening of diffraction rings without large per- 
centages of voids such as the. crystallite theory 
seemed to require. 

The network theory is modified and strengthened 
by the vitron concept because the existing trouble- 
some requirement for first-order randomness in the 
network is thereby eliminated—a randomness that 
is not entirely compatible with some of the well- 
established properties of glasses and is being chal- 
lenged by electron and neutron diffraction data 
which evidence greater order than did X-ray data. 
Further, the idea of unbroken continuity is replaced 
by the much more useful concept of a quasi-discrete 
or modulated continuity. 

In addition, the vitron concept readily diverts 
excessive attention from the old idea that glass is 
merely supercooled liquid. Thus it accords well with 
certain thermodynamical and other considerations 
that seem to have led Parks and Huffman [30], and 
later Berger [31], to suggest that glass is a fourth 
state of matter resembling but definitely differing 
from the erystalline state rather than a_ twilight 
zone intermediate between liquids and crystals. 
Because in some respects the glassy state differs so 
very little from the ervystalline state, it has always 
been difficult to accept as fully satisfactory the 
undercooled liquid ideas or to account for the sta- 
bility and reproducibility of definite degrees of 
major disorder in networks in the absence of coercive 
fields. An important difference between the present 
vitron concept and Berger’s ideas, as usually under- 
stood, is that the vitron or glassy state forms at 
very high temperatures in the melt rather than at 
or near annealing or softening temperatures. 

The vitron concept of stressed clusters of cages 
having a dodecahedral arrangement of pentagonal 
rings of slightly deformed tetrahedra, may be ap- 
plicable to other than silicate glasses. For boric 
oxide glass with boron in four-coordination state, 
it can be shown that a reasonable density is given by 
such an arrangement starting with nearest-neighbor 
distances to accord with X-ray evidence. 

In the broader field of glasses in general, it is 
possible that the vitron concept (with other non- 
crystal symmetries) may open new vistas for study 
and investigation. Perhaps it will com- 
parable concepts relating to the constitution of 
liquids and solutions. Particularly in the field of 
rheology, it may prove useful and possibly it may 
relate to crystallinity in the field of high polymers 
It would seem that some writers have felt the need 
for and referred to somewhat similar concepts. 
Machatschki [32] speaks of ‘‘precrvstalline”’ arrange- 
ments; Debye and Menke [33] deseribe a ‘‘quasi- 
crystalline’ arrangement instead of a disordered 
atomic distribution in liquid mereury; and Matossi 
[34] infers that in fused silica there are regions of 
distorted order separated by places of still greater 
lattice distortion. Aside from an investigation of 
its possible applications, the idea seems to call for a 
mathematical study of the various classes of sym- 
metry that are adapted to noncrystal-forming ar- 
rangements of matter in space, with a degree of 
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| periodic order and regularity roughly comparable 
to but much more limited than those found jp 
crystals. 


13. Appendix 
13.1. Union and Attachment Between Vitrons 


Some possibilities concerning the growth and unit- 
ing of clusters, or vitrons, may warrant examination 
in more detail. Consider a sphere inscribed in g 
regular pentagonal dodecahedron. It can be sur. 
rounded by 12 equal spheres tangent to the first on 
each of the 12 pentagonal faces to form a cluster of 13 
spheres of which only the central sphere touches more 
than one other. In the next shell one may add 12 
spheres along diameters through the 12 of the first 
shell and also 20 on diameters through the inner 
sphere and through the 20 apices of the original 
polygon. The cluster is now composed of 45 spheres, 
This is a packing that is regular but uniquely refer- 
able to only one definite center. 

In this type of packing of rigid incompressible 
spheres or cages the density decreases toward the 
periphery because of voids between cages as the 
cluster grows. If the cages are deformable, however, 
there can be a small but limited growth of such a 
system with sufficient distortion of cages in the outer 
shells to the interstices and form a cluster 
approximating uniform density. In the case of dodee- 
ahedrons of silicon and oxygen, the chemical bond 
between these elements is very strong and these are 
the forces that tend to maintain vitron structures 
after their formation at high temperatures. The ex- 
tent to which they can be maintained is, however, 
definitely limited. Important possibilities in this 
tvpe of vitron are clusters of 1, or 13, or 45 dodee- 
abedra. 

In grouping dodecahedral elements to form vitrons 
one can arrange them, as just described, around (1) 
a central elemental volume, (2) a common face of 2 
elements, (3) a common edge of 3 elements, or (4) a 
common apex of 4 elements. The first type permits 
the nearest approach to spherical clusters. The 
second is perhaps the least spherical; it can lead to 
chain structure. The third can manifest a evlindrical 
tendeney ; also it seems adapted to sheet or laminated 
structure. The fourth is somewhat tetrahedral but 
the surfaces and edges of the whole cluster are so 
decidedly curved that sphericity is approximated as 
the vitron grows. 

In studying these tvpes of arrangement, the char- 
acteristics of attachable and unitable vitrons were 
computed for several more or less spherical clusters 


close 


consisting of from 1 to 54 elements of volume. The 
results are summarized in table 1. In addition to the 
tabulated characteristics of vitrons in table 1, there 


are other factors that affect their degree of union or 
attachabilitv. The single element as tabulated can 
unite any two clusters each of which has one free 
penta face. There are other useful reduced single 
elements such as (SigQ).)3) which can be joined 
by union to 2 adjacent penta faces, as presented by a 
2-cluster, to form a 3-cluster: and = (Si,O,,)8F°> 
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which is unitable to 3 adjacent faces, as presented by 
a 3-cluster, to form a 4-cluster vitron, ete. 

A dual element as tabulated can unite with one 
single element as reduced for adjacent-2-face attach- 
ability and to one single element reduced to adjacent- 
3-face attachability, to form a 4-element cluster as 
tabulated; that is 1(SixpOw) + 1(SigOis) + 1 (SigOj,) 
1(SigsO7,)"*-. =A reduced dual element that 
unitable to two adjacent faces is (SijgOQxs) 3f7"'® 

A triple element cluster as tabulated presents 3 
convex 4-faces, or 2 reentrant 3-faces plus 3 convex 
9-faces, and can unite clusters presenting corre- 
sponding] opposite combinations of penta faces. 

The tabulated 4-element vitron can, among other 
possibilities, be intermediate (that is form or grow) 
between 4 clusters each presenting a reentrant 3-face. 

The 5-element vitron can unite or have attached to 
it 2 small clusters each with reentrant 3-face plus 3 
small clusters each with reentrant 2-face. 

The 7-element cluster might unite or have attached 
to it 5 small clusters each with reentrant 2-face. 

The 13-cluster vitron is very nearly spherical. It is 
directly and fully unitable through its 12 free faces; 
or it can accept unions simultaneously at 8 of the 20 
places where it has reentrant 3-face combinations. 
A few oxvgen connections are also available. 

The 14-element vitron presents 12 places of re- 
entrant 4-face and a few opportunities for connecti- 
bility through oxygens; or it presents convex 2-face 
possibilities of union in 6 places, together with added 
attachability through oxygens. 

The 26-element vitron is directly unitable on each 
of 12 faces and has added attachability through 
oxygen. 

The 45-element vitron is the largest of the spherical 
or element-center tvpe that has been studied. It is 
directly unitable through each of 12 penta faces and 
also attachable through 20 oxygens at places equally 
distributed over the vitron surface. 

The 54-element tetrahedral vitron is unique among 
those in table 1 in that it presents no penta faces for 
union with other vitrons. All connections must be 
through oxygen at 4 triple edges, 12 double edges, 
and 12 single edges. 


IS 


13.2. “Unit Cell’’ of Vitreous Silica 


The symmetry within vitrons of limited size has 
different aspects for the different arrangements. 
For type 1 clusters with a central element there is a 
center of symmetry and the whole vitron has, 
basically, all of the axes that were mentioned in 
section 2 in describing a dodecahedral element (see 
fig. 7); also there are quasi-reflecting planes similar to 
those for single elements. It is useful, in attempting 
to discuss the concept of a “unit cell’ for vitreous 
silica, to note that it is possible to select 3 mutually 
perpendicular axes running through oxygen atoms 
and intersecting at the center of an element. As one 
travels along any of these orthogonal axes of a type 
| cluster, one proceeds 4.19 A from origin to first 
oxygen, and 4.91 A further to first silicon, then 
orients 90° around the axis and proceeds 4.91 A to 
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second oxygen, and then 4.19 A to center of second 
elemental cage. The second element, however, is at 
90° with the element of origin and one must repeat 
the four above-mentioned steps to arrive at an 
element center apparently identical with the origin 
(but in reality at 180° therewith). In other words 
the unit cell in the sense used in crystallography 
might seem to be 36.4 A on each edge, if the structure 
could grow to such an extent. Such a cell would 
contain 960 SiO,. Apparently, one-eighth of this 
cell with edge of 18.2 A and containing 120 SiO, 
might be representative of the structure if suitable 
orientations and inversions are made. 

The estimate of a unit cell just given does not seem 
satisfactory, however, because along axes parallel to, 
but somewhat displaced from, the principal axes the 
repeat distances are not always 36.4 A. Other 
repeat distances of 14.26 and 24.34 A are easily 
found. 

If vitrons of type 4 with silicon as a center are 
considered, it is obvious that there are four axes 
equally distributed in space and making angles of 
109° 28’ with each other as in a tetrabedron. Along 
any of these axes of a cluster of type 4, one proceeds 
1.60 A from origin to first oxygen, 1.60 A to first 
silicon, 8.97 to second silicon, then orients 60° 
around the axis and proceeds to repeat these 3 steps 
in the same order to arrive at a silicon apparently 
identical with the origin (but in reality at 120° 
therewith). The apparent repeat distance is 24.34 A. 

The vitron of type 3 with a common edge and 
oxygen at its center, has a principal axis with a 
repeat distance of 24.34 A, numerically equal to that 
just found for the tetrahedral axes of type 4 vitrons. 
Perpendicular to the principal axis there are three 
axes at 120° with each other. Along each of these 
the repeat distance is 36.4 A, numerically equal to 
that found for the orthogonal axes of type 1 vitrons. 

More complicated is the vitron of type 2 with 
center in a face. Here one has a principal axis and 2 
groups of 5 axes each in the plane of the central or 
common face. Along the principal axis, perpendic- 
ular to the central face, one proceeds 7.13 A to the 
center of the opposite face, then orients 36° and 
proceeds 7.13 A to a face center that appears as did 
the origin (but in reality at 72° therewith). The 
apparent repeat distance is 14.26 A. For one group 
of 5 axes, through the 5 oxygens of the central face, 
the repeat distance is 36.4 A; and for the other 
group, through the 5 silicons, the repeat distance is 
also 36.4 A. 

The repeat distance of 36.40 A along a two-fold 
axis through an element center (type 1) must be 
multiplied by 2 to reach an element center identical 
with the origin in true orientation; the repeat dis- 
tance of 24.34 A along a three-fold axis through a 
silicon center (type 4) must be multiplied by 3 to 
reach a silicon identical with the origin in true 
orientation; and the repeat distance of 14.26 A along 
a five-fold axis through a face center (type 2) must 
be multiplied by 5 to reach a face center identical 
with the origin in true orientation. The average of 
these 3 estimates gives the true repeat distance for 





the proposed structure of vitreous silica along any of 


its multiple axes as 72.4 A and the corresponding cell 
would contain 7680 SiO,. This can be interpreted 
as the size of the unit cell of vitreous silica that 


would be found if the structure were capable of 


growing to such an extent. 
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